2urn:lsid:arphahub.com:pub:45048D35-BB1D-5CE8-9668-537E44BD4C7Eurn:lsid:zoobank.org:pub:91BD42D4-90F1-4B45-9350-EEF175B1727AZooKeysZK1313-29891313-2970Pensoft Publishers10.3897/zookeys.625.91169116Research ArticleCervidaeColeopteraBiodiversity & ConservationJapanEffects of an increase in population of sika deer on beetle communities in deciduous forestsIidaTaichitaiyakiiida@gmail.comhttps://orcid.org/0000-0003-1868-36821SogaMasashi2KoikeShinsuke1Graduate School of Agriculture, Tokyo University of Agriculture and Technology, 3-5-8, Saiwai, Fuchu 183-8509, JapanGraduate School of Agriculture, Tokyo University of Agriculture and TechnologyFuchuJapanSchool of Engineering, The University of Tokyo, 7-3-1, Hongo, Bunkyo, Tokyo, 113-8656, JapanSchool of Engineering, The University of TokyoTokyoJapan
2016191020166256785FFDEFFAA-FFED-C857-A936-FFA9FFA0A8117AE785C5-2827-4A44-A653-A25BFB5694C84959081005201613092016Taichi Iida, Masashi Soga, Shinsuke KoikeThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/7AE785C5-2827-4A44-A653-A25BFB5694C8
The overabundance of large herbivores is now recognized as a serious ecological problem. However, the resulting ecological consequences remain poorly understood. The ecological effects of an increase in sika deer, Cervusnippon Temminck (Cervidae), on three insect groups of beetles was investigated: ground beetles (Carabidae), carrion beetles (Silphidae), and dung beetles (Scarabaeidae and Geotrupidae) on Nakanoshima Island, Hokkaido, northern Japan. We collected beetles on Nakanoshima Island (experimental site) and lakeshore areas (control site) and compared the species richness, abundance, diversity index, and community composition of beetles between the sites. Results showed that although both species diversity and abundance of carabid beetles were significantly higher at the lakeshore site, those of dung and carrion beetles were higher at the island site. It was additionally observed that abundance of larger carabid beetles was higher at the lakeshore site, whereas that of small-sized carabid beetles did not differ between the lakeshore and island sites. For dung beetles, abundance of smaller species was higher at the island site, whereas that of large species did not differ between the lakeshore and island sites. Abundance of two body sizes (small and large) of carrion beetles were both higher at the island site. Overall, the findings of this study demonstrated that an increase in deer population altered the insect assemblages at an island scale, suggesting further changes in ecosystem functions and services in this region.
Iida T, Soga M, Koike S (2016) Effects of an increase in population of sika deer on beetle communities in deciduous forests. ZooKeys 625: 67–85. doi: 10.3897/zookeys.625.9116
Introduction
The overabundance of large herbivores is now recognized as one of the serious ecological issues worldwide, especially in the northern hemisphere (Gill 1992, Côté et al. 2004, Takatsuki 2009). Indeed, there is mounting evidence demonstrating the serious consequences of large herbivore overabundance on forest ecosystems (Takatsuki 2009, Foster et al. 2014). As large herbivores selectively browse palatable species, browsing has negative effects on understory structure and species composition (Rooney 2009, Martin et al. 2011, Tanentzap et al. 2011). Large herbivores additionally alter forest structure (e.g. Kanda et al. 2005, White 2012) by limiting tree regeneration (Akashi and Nakashizuka 1999, Horsley et al. 2003), thereby inducing a cascading effect on animal species (Kanda et al. 2005, Bressette et al. 2012). Indeed, it has been reported that smaller herbivorous invertebrates are negatively influenced by large herbivore overabundance because of interspecies competition for food resources (Gómez and González-Megías 2002, Wheatall et al. 2013). In addition, browsing results in the degradation of habitat for invertebrates and birds (Stewart 2001, Allombert et al. 2005, Chollet and Martin 2013).
To evaluate the ecological impacts of large herbivore overabundance on ecosystems, previous studies have commonly used manipulations with inclusion/exclusion treatments of large herbivores (Goméz and González Megías 2002, Mysterud et al. 2010, Holt et al. 2011, Bush et al. 2012). However, these previous studies have two major limitations. First, as past studies investigated the effects of large herbivore overabundance up to 10 years (e.g. Suominen et al. 2003, Beguin et al. 2011), the effects on ecosystems were unlikely to be completely detected. Indeed, it is established that the effects of browsing by large herbivores on vegetation structure can persist for an extended period (Tanentzap et al. 2011, Nuttle et al. 2014). Second, to understand the responses of ecosystems to perturbation, it is necessary to conduct a large-scale experiment (Carpenter 1998), whereas previous studies have employed relatively small-scale experiments (e.g. Dennis et al. 1998, Gómez and González-Megías 2002; Iida et al. 2016). Nevertheless, field studies that employed a relatively large-scale experiment remain scarce, with most having been conducted in Canadian coniferous forest (e.g. Allombert et al. 2005, Martin et al. 2010).
Lake Toya located in western Hokkaido, northern Japan, provides an ideal study site to investigate the long-term impacts of deer overabundance on ecosystems (Fig. 1). On Nakanoshima Island (476.7 ha, hereafter, island), which is one of the islands situated within Lake Toya, three sika deer individuals [Cervusnippon Temminck (Cervidae)] were introduced approximately 50 years ago and deer density at the island has now reached over 50 deer/km2 (Ikeda et al. 2013). Takahashi and Kaji (2001) reported a significant impoverishment of understory vegetation at the island because of the rapid growth in the deer population (Fig. 1B), suggesting a further consequence for other animal species. In the areas around Lake Toya (hereafter, lakeshore), and at the time of the current study, deer density was approximately 1.1 × 10−2 deer/km2 (Akaba et al. 2014), which is considerably lower than the density at the island site (Fig. 1C). As the island is geographically isolated from lakeshore areas, the island ecosystem could be considered as an experimental site for examining the effects of an increase in deer population on natural ecosystems.
The location and forest floor of the study area. A The location of Lake Toya B Forest floor of Nakanoshima Island C Forest floor of the lakeshore of Lake Toya. The map of Lake Toya was modified from the Geospatial Information Authority of Japan 2015.
https://binary.pensoft.net/fig/107248
Here, we evaluated the effects of an increase in sika deer population on four taxonomic groups of beetles: carabid (Carabidae), carrion (Silphidae) and dung beetles (Scarabaeidae and Geotrupidae). These beetle groups were selected for several reasons. First, as these beetles inhabit the forest floor and are known to be sensitive to microclimatic changes (Rainio and Niemelä 2003, Arellano et al. 2005), they are likely to sharply respond to the increase in the large herbivores (Rooney 2009, Yamada and Takatsuki 2015). Second, many beetles display species-specific ecological traits, which facilitates the investigation of the association between species life-history traits and their responses to large herbivore overabundance (Sumways 2007, Bachand et al. 2014). Third, the clarification of the responses of beetles to environmental changes is crucial for understanding the overall changes in forest ecosystems, as beetles constitute a large proportion of biomass in forest ecosystems and play important roles that maintain ecosystem functions, such as nutrient cycling and pollination (Kevan and Baker 1983, Speight et al. 2008, Barton et al. 2013).
Materials and methodsStudy area
Our study area is located in Western Hokkaido, northern Japan. The deer density at the island site (42°36'N, 140°51'E (DDM)) is dramatically higher than that at the lakeshore site because of artificial introduction. The mean annual temperature of this area is 7.3 °C and the mean monthly temperature ranges from −5.1 °C to 20.2 °C. The mean annual precipitation is 984.8 mm and the mean annual snow depth is 30 cm. The study area is situated in a deciduous forest dominated by Quercuscrispula Blume (Fagaceae); Kalopanaxseptemlobus (Thunb.) Koidz. (Araliaceae); Magnoliaobovata Thunb. (Magnoliaceae); Acerpictum Thunb. (Sapindaceae); Maackiaamurensis Rupr. et Maxim. (Fabaceae); Tiliajaponica (Miq.) Simonk. (Tiliaceae) and Ostryajaponica Sarg. (Betulaceae) (Kaji et al. 1991). Because of overgrazing at the island site, the structure of understory vegetation and the forest floor differ between the island and lakeshore sites (Fig. 1), with the island site dominated by Seneciocannabifolius Less. (Asteraceae); Saginajaponica (Sw.) Ohwi (Caryophyllaceae); Pachysandraterminalis Siebold et Zucc. (Buxaceae) and Chloranthusserratus (Thunb.) Roem. et Schult. (Chloranthaceae), which sika deer find unpalatable (Takahashi and Kaji 2001, Miyaki and Kaji 2009). The forest floor of the lakeshore site is dominated by Dryopteriscrassirhizoma Nakai (Dryopteridaceae); Cardamineleucantha (Tausch) O.E. Schulz (Brassicaceae); Cephalotaxusharringtonia (Knight ex Forbes) K.Koch var. nana (Nakai) Rehder (Cephalotaxaceae) and Sasasenanensis (Franch. et Sav.) Rehder (Poaceae) (observed by authors). Although other large animals including Vulpesvulpes Linnaeus (Canidae); Nyctereutesprocyonoides Gray (Canidae) and Procyonlotor (Linnaeus) (Canidae) occur at the lakeshore site, their abundance is negligible compared to sika deer (Akaba et al. 2014). Besides the sika deer, the island site contains no medium- or large-sized mammals (Akaba et al. 2014). We established 30 sampling plots which are 100 m apart from each other (along a 3 km sampling transect) at each of the island and lakeshore sites (42°34'N, 141°54'E (DDM)).
Beetles sampling
At both sites, we sampled carabid, carrion and dung beetles using pitfall traps baited with cattle dung and fermented milk. Fermented milk is one of the major bates in collecting ground-dwelling beetles in Japan (Suttiprapan and Nakamura 2007). Pitfall traps were constructed using plastic containers (22.5 cm diameter and 26.6 cm deep) and plastic cups (8.3 cm diameter and 11.5 cm deep). Each container of fermented milk was buried to the rim in the ground. Plastic cups containing cattle dung were then hung inside the container using wires. We set plastic roofs on the traps to prevent interference from rain and fallen leaves. These traps were set at a density of one trap per plot. Surveys were conducted during early September 2012 and 2013.
Species classification
For each beetles group, the insect species collected in the field were divided into different size groups according to body length (see Fujita et al. 2008, Koike et al. 2014). Information on body size was collected from Ueno et al. (1985). Carabid and carrion beetles were divided into groups for small, medium or large species (small, <10 mm; medium, ≥10 and <20 mm and large, ≥20 mm), and dung beetles were divided into groups for small or large species (small, <10 mm; large, ≥10 mm). However, because carrion beetle species which are divided into small size group were not sampled, carrion beetles were divided into the medium or large size groups.
Data analyses
All analyses were performed using R ver. 3.2.1 (R core team 2015). In order to test the differences in abundance and species richness of each taxonomic group or body length between the island and lakeshore sites, we used
generalized linear models
(GLMs), using the ‘glm’ function. For GLMs, the abundance and species richness of each taxonomic or body size group were used by a response variable with a Poisson distribution and a log link function. Sampling sites were used as categorical explanatory variables (we used the lakeshore site as a reference). To test the difference in species diversity between two sites, Shannon-Wiener diversity index of each taxonomic group which were calculated by using package ‘vegan’ (Oksanen et al. 2015) were also used as a response variable with a Gaussian (normal) distribution and an identity function. Species diversity indices of each functional group were not calculated and were excluded from the analysis because of the small sample size. Because we could not identify Synuchus spp. to the species level, these species were excluded from the analysis of species richness and diversity. We did not find large carabid species in the island sites, which hinders the fitting of our data using GLMs. Hence, we excluded large carabid species in this analysis. Because of small sampling sizes, we used ‘zeroinfl’ function in the package ‘pscal’ for abundance and species richness of small carabid species, and medium and large carrion beetles. We calculated estimated abundance, species richness and Shannon-Wiener indices using the estimates of parameters in GLMs.
To investigate the difference in species composition between the island and lakeshore sites, we performed a
non-metric multidimensional scaling ordination
(NMDS) using the metaMDS function with the Bray–Curtis measure within the package ‘labdsv’ (Roberts 2015) as well as cluster analysis. The validity of clustering was evaluated using the Calinski–Harabasz criterion. The indicator species value (IndVal) was calculated to identify species, which are indicators of assemblages of each site as illustrated by NMDS and cluster analysis.
Results
3,876 individuals in total were collected, comprising 824 carabid beetles (18 species), 148 carrion beetles (four species) and 2,902 dung beetles (five species) (see Suppl. material 1).
Comparisons between species richness, abundance, and diversity index of insect species between the island and lakeshore sites were summarized in Figs 2–4 and Tables 1–3 (see also Supplementary Tables 2–4 for the differences between estimated and observed values). For carabid beetles, the abundance of individuals at the lakeshore site was higher than that at the island site (p < 0.001; Fig. 2A; Table 1). Although the abundance of medium species were higher at the lakeshore site than at the island site (p < 0.001), that of small species did not significantly differ between the two sites (p = 0.24; Fig. 2A; Table 1). The species richness of all carabid species (p = 0.49), small (p = 0.42) and medium species (p = 0.38) did not significantly differ between the island and lakeshore sites (Fig. 3A; Table 2). The species diversity index was higher at the lakeshore site than at the island site (p = 0.05; Fig. 4A; Table 3).
Abundance of each insect taxonomic and functional group per plot at the island and lakeshore sites estimated by using GLMs. A Carabid B Carrion C Dung beetles. Carabid and carrion beetles were classified into small, medium and large species and dung beetles were classified into small and large species (see the main text). Asterisks indicate a significant difference (p < 0.05).
Species richness of each insect taxonomic and functional group per plot at the island and lakeshore sites estimated by using GLMs. A Carabid B Carrion and C Dung beetles. Carabid and carrion beetles were classified into small, medium and large species and dung beetles were classified into small and large species (see the main text). Asterisks indicate a significant difference (p < 0.05).
Shannon–Wiener index of each insect taxonomic group per plot at the island and lakeshore sites estimated by using GLMs. A Carabid B Carrion and C Dung beetles. Asterisks indicate a significant difference (p < 0.05).
https://binary.pensoft.net/fig/107251
Estimates of parameter differences between sites (reference = lakeshore site, i.e., coefficients of lakeshore site are zero),
standard errors
(SEs), z values, and p values in generalized linear models of each of taxonomic group and body size group abundance. Bold letters represent p < 0.05.
Parameters
Estimates
SEs
z values
p values
Carabid beetles
Whole species
(Intercept)
2.892
0.043
67.271
< 0.001
Island
-0.596
0.072
-8.266
< 0.001
Small species
(Intercept)
0.170
0.484
0.351
0.73
Island
0.799
0.682
1.172
0.24
Medium species
(Intercept)
2.526
0.052
48.911
< 0.001
Island
-0.707
0.090
-7.864
< 0.001
Carrion beetles
Whole species
(Intercept)
-0.511
0.236
-2.168
0.03
Island
1.985
0.251
7.897
< 0.001
Medium species
(Intercept)
-0.051
0.316
-0.160
0.88
Island
1.280
0.340
3.767
< 0.001
Large species
(Intercept)
-2.226
0.736
-3.023
0.003
Island
3.536
0.745
4.758
< 0.001
Dung beetles
Whole species
(Intercept)
2.686
0.048
56.330
< 0.001
Island
1.776
0.052
34.450
< 0.001
Small species
(Intercept)
2.230
0.060
37.250
< 0.001
Island
2.165
0.063
34.260
< 0.001
Large species
(Intercept)
1.680
0.079
21.319
< 0.001
Island
0.031
0.111
0.276
0.78
Estimates of parameter differences between sites (reference = lakeshore site, i.e., coefficients of lakeshore site are zero),
standard errors
(SEs), z values, and p values in generalized linear models of each of taxonomic group and body size group species richness. Bold letters represent p < 0.05.
Parameters
Estimates
SEs
z values
p values
Carabid beetles
Whole species
(Intercept)
0.588
0.136
4.319
< 0.001
Island
-0.139
0.200
-0.696
0.49
Small species
(Intercept)
-2.015
0.500
-4.028
< 0.001
Island
-0.693
0.866
-0.800
0.42
Medium species
(Intercept)
0.210
0.164
1.276
0.20
Island
0.196
0.222
0.882
0.38
Carrion beetles
Whole species
(Intercept)
-0.836
0.277
-3.015
0.003
Island
1.149
0.318
3.609
< 0.001
Medium species
(Intercept)
-1.003
0.302
-3.327
< 0.001
Island
0.738
0.367
2.012
0.04
Large species
(Intercept)
-2.708
0.707
-3.830
< 0.001
Island
2.197
0.745
2.948
0.003
Dung beetles
Whole species
(Intercept)
0.758
0.125
6.061
< 0.001
Island
0.426
0.161
2.651
0.008
Small species
(Intercept)
0.154
0.169
0.912
0.36
Island
0.693
0.207
3.348
< 0.001
Large species
(Intercept)
-0.034
0.186
-0.183
0.86
Island
-0.035
0.265
-0.132
0.90
Estimates of parameter differences between sites (reference = lakeshore site, i.e., coefficients of lakeshore site are zero),
standard errors
(SEs), t values, and p values in generalized linear models of each of taxonomic group and body size group Shannon-Wiener diversity index. Bold letters represent p < 0.05.
Parameters
Estimates
SEs
t values
p values
Carabid beetles
(Intercept)
0.841
0.105
7.995
< 0.001
Island
-0.302
0.149
-2.030
0.05
Carrion beetles
(Intercept)
0.067
0.069
0.971
0.34
Island
0.285
0.097
2.935
0.005
Dung beetles
(Intercept)
0.828
0.051
16.300
< 0.001
Island
0.347
0.072
4.828
< 0.001
The indicator species values of each species. Bold letters represent p < 0.05.
Taxa
Community
IndVals
p values
Caccobiusjessoensis Harold
Scarabaeidae
island
0.988
0.001
Onthophagusater Waterhouse
Scarabaeidae
island
0.822
0.001
Eusilphajaponica (Motschulsky)
Silphidae
island
0.604
0.001
Silphaperforata Gebler
Silphidae
island
0.405
0.01
Pterostichusleptis Bates
Carabidae
island
0.375
0.001
Liatongusphanaeoides (Westwood)
island
0.345
0.001
Nicrophorusquadripunctatus Kraatz
Silphidae
island
0.177
0.23
Coprisochus Motschulsky
Scarabaeidae
island
0.172
0.02
Chlaeniuspallipes Gebler
island
0.115
0.39
Pterostichusplanicollis (Motschulsky)
Carabidae
island
0.103
0.11
Pterostichusyoritomus Bates
Carabidae
island
0.069
0.21
Pterostichushaptoderoides (Tschitscherin)
Carabidae
island
0.034
0.49
Chlaeniusvariicornis Morawitz
island
0.034
0.49
Lithochlaeniusnoguchii (Bates)
island
0.034
0.49
Pterostichussamurai (Lutshnik)
Carabidae
island
0.034
0.48
Hemicarabustuberculosus (Dejean et Boisduval)
island
0.034
0.48
Oiceoptomathoracicum (Linnaeus)
island
0.034
0.47
Pterostichusprolongatus Morawitz
Carabidae
island
0.031
0.74
Synuchus spp. Gyllenhal
Carabidae
lakeshore
0.827
0.001
Pterostichusthunbergi Morawitz
Carabidae
lakeshore
0.806
0.001
Geotrupeslaevistriatus Motschulsky
Geotrupidae
lakeshore
0.477
0.56
Leptocarabusarboreus (Lewis)
Carabidae
lakeshore
0.355
0.001
Trichotichnuslongitarsis Morawitz
Carabidae
lakeshore
0.075
0.58
Pterostichusorientalis (Motschulsky)
Carabidae
lakeshore
0.065
0.50
Damasterblaptoides Kollar
Carabidae
lakeshore
0.065
0.50
Cychrus morawitzi Gehin
lakeshore
0.065
0.49
Leptocarabusopaculus (Putzeys)
Carabidae
lakeshore
0.032
1.00
The abundance and species richness of carrion beetles at the island site were higher than those at the lakeshore site (p < 0.001; Figs 2B and 3B; Tables 1 and 2). The abundance of the two body size groups (medium and large species) were higher at the island site than at the lakeshore site (p < 0.001; Fig. 2B; Table 1). The species richness of the two body size groups were also higher at the island site than at the lakeshore site (medium: p = 0.04; large: p = 0.003). The mean species diversity index was higher at the island site than at the lakeshore site (p = 0.005; Fig. 4B; Table 3).
For dung beetles, the abundance (p < 0.001) and species richness (p = 0.008) at the island site were higher than those at the lakeshore site (Figs 2C, 3C; Tables 1 and 2). Although the abundance and species richness of small species were similarly higher at the island site than at the lakeshore site (p < 0.001; Figs 2C, 3C; Tables 1and 2), those were not the case for large species (abundance: p = 0.78, species: p = 0.90; Fig. 2C, 3C; Tables 1 and 2). The mean species diversity index was higher at the island site than at the lakeshore site (p < 0.001; Fig. 4C; Table 3).
NMDS and cluster analysis identified a significant difference in species composition between the island and lakeshore sites (Fig. 5). In addition, according to indicator species analysis, indicator species were different between the island and lakeshore sites. Indicator species analysis showed that Caccobiusjessoensis Harold; Onthophagusater Waterhouse; Liatongusphanaeoides (Westwood); Coprisochus Motschulsky; Eusilphajaponica (Motschulsky); Silphaperforata Gebler and Pterostichusleptis Bates had significant IndVals for the island site (Table 4). Conversely, Synuchus spp. Gyllenhal; Pterostichusthunbergi Morawitz and Leptocarabusarboreus (Lewis) had significant IndVals for the lakeshore site.
Non-metric multidimensional scaling ordination of each plot of the island and lakeshore sites. From I1 to I30: plots at the island site; from L1 to L30: plots at the lakeshore site. The differences of color and symbol represent the result of clustering.
https://binary.pensoft.net/fig/107252Discussion
The current study demonstrated that an increase in sika deer at Lake Toya significantly changed both the abundance and species richness of beetle species belonging to three different taxonomic groups. The abundance and diversity index of carabid beetles at the island site were significantly lower than that at the lakeshore site (Figs 2, 4). In addition, we observed lower species richness and diversity indices at the island site; however they were not statistically significant. As carabid beetles utilize the forest floor as a habitat, the decline of understory vegetation because of overbrowsing is likely to result in the alteration of micro-climate conditions of their habitat, such as a decrease in humidity and an increase in light availability, and a higher mortality by predation (Rooney and Waller 2003, Melis et al. 2007, Cerda et al. 2015). For example, Melis et al. (2007) reported that shade tolerant and hygrophilous carabid species were negatively affected by moose browsing as browsing resulted in a decrease of bilberry coverage and changed the humidity and light intensity. Although we did not measure vegetation and other environmental variables, a significant decline in understory cover has been reported at the island site (Takahashi and Kaji 2001; Fig. 1B).
Conversely, carrion and dung beetles responded positively to deer overabundance, with abundance, species richness and the diversity index higher at the island site (Figs 2–4). Adult and larvae dung beetles utilize the faces of mammals as their main food resource (Andresen and Laurance 2007, Nichols et al. 2009). Hence, deer overabundance at the island site is likely to produce a greater quantity of deer faces, thereby facilitating a more abundant food source for dung beetles as compared to the lakeshore site. Indeed, Kanda et al. (2005) reported a positive relationship between dung beetle abundance and deer density. In addition, as adult and larvae carrion beetles utilize animal carcasses as a food resource (Scott 1998, Dekeirsschieter et al. 2011), deer overabundance may provide greater abundance of food resources for carrion beetles. In support of these assertions, it has similarly been reported that a deer carcass increases the activity of necrophagous beetles (Melis et al. 2004).
Although small carabid species were not affected by deer overabundance, medium and large carabid species were negatively affected by deer overabundance (note: large carabid species were not sampled at the island site). This result suggests that larger carabid species are more sensitive to changing habitat condition than smaller species (Magura et al. 2006, Jelaska and Durbes˘ić 2009). Because the understory vegetation protects these species against extreme microclimate conditions (Ikeda et al. 2005), understory decline because of large herbivore overbrowsing results in a harsh ground floor environment for larger carabid species (Melis et al. 2007). Indeed, a higher susceptibility of larger species to other environmental changes, including habitat fragmentation, deforestation and urbanization, has been reported (Magura et al. 2006, Fujita et al. 2008).
Although the abundance of small dung beetle species was significantly higher at the island site than at the lakeshore site, that of large species did not differ between the two sites. This result suggests that small species are likely to favor a deer abundant environment. In contrast to our results, in Japan, Koike et al. (2014) observed that small species did not prefer environments with higher deer density. One possible reason for the difference would be the difference of dung usage among study species. Koike et al. (2014) found that small species tended to be dwellers which simply lay their eggs into the dung on the ground (Camberfort and Hanski 1991). Therefore, drying because of sunlight exposure, which results from deer overbrowsing, negatively affects these species. Conversely, in our study, the majority of small species were tunnelers, which are species burying dung under the ground before oviposition (Camberfort and Hanski 1991). Because of dung burying, tunnelers are not influenced by the drying of the ground surface.
In our study, NMDS showed a difference in species composition between the island and lakeshore sites (Fig. 5). In addition, indicator species analysis suggested that whereas beetle assemblage at the lakeshore site was characterized by carabid beetle species, that at the island site was characterized by carrion and dung beetle species (Table 4). These results, along with other reported findings indicate a biological homogenization of insect communities at the island site (McKinney and Lockwood 1999). In our study, large herbivore overabundance is likely to cause the homogenization of the beetle assemblage through the increased number of carrion and dung beetles and the decreased number of carabid beetles (Figs 2–4). Moreover, our findings suggest a shift in body size in insect communities; as the number of deer increase, insect communities are more likely to be dominated by smaller individuals.
Although our study brings a valuable contribution, the survey data is limited in its ability for generalizing the current results. Indeed, since the data used in this study was obtained from short term surveys with a limited number of samples (see also the species accumulation curves in Suppl. material 5: Figure 1), our results would provide only partial support for our hypotheses. Thus, further longitudinal studies must be needed to understand the ecological impacts of deer overabundance on insect communities more comprehensively.
Conclusions
In the present study, evidence is provided that an increase in deer population altered species richness, abundance and diversity of beetles within three different taxonomical groups. Whether such changes affect the ecosystem functions provided by these beetles is unknown. Nevertheless, the observed change in compositions of the three taxonomic beetles groups raises the potential that ecosystem functions may be altered through cascading effects (Larsen et al. 2005, Rouabah et al. 2014). To conserve and maintain overall forest ecosystems, an investigation into the responses of other taxonomic species, including birds or amphibians, to the rapid increase in deer population and its impacts on ecosystem functions is required.
Acknowledgements
We thank K. Kaji, T. Yoshida, T. Ikeda, T. Hino and Toyako Town for providing valuable information of the study areas. We thank T Ikeda for providing the picture of study sites (Fig. 1). This study was partly funded by JSPS Grant-in-Aid for Scientific Research (No. 25292085, No. 16H02555) and 26th Pro-Natura Fund, Japan in 2015.
ReferencesAkabaSHinoTYoshidaT (2014) Impacts of high density of sika deer on dung beetle community.58: 269–274. doi: 10.1303/jjaez.2014.269AkashiNNakashizukaT (1999) Effects of bark-stripping by sika deer (Cervusnippon) on population dynamics of a mixed forest in Japan.113: 75–82. doi: 10.1016/S0378-1127(98)00415-0AllombertSStocktonSMartinJ-L (2005) A natural experiment on the impact of overabundant deer on forest invertebrates.19: 1917–1929. doi: 10.1111/j.1523-1739.2005.00280.xAndresenELauranceSGW (2007) Possible indirect effects of mammal hunting on dung beetle assemblages in Panama.39: 141–146. doi: 10.1111/j.1744-7429.2006.00239.xArellanoLFavilaMEHuertaC (2005) Diversity of dung and carrion beetles in a disturbed Mexican tropical montane cloud forest and on shade coffee plantations.14: 601–615. doi: 10.1007/s10531-004-3918-3BachandMPellerinSMorettiMAubinITremblayJ-PCôtéSDPoulinM (2014) Functional responses and resilience of boreal forest ecosystem after reduction of deer density.9: e90437. doi: 10.1371/journal.pone.0090437BartonPSCunninghamS a.LindenmayerDBManningAD (2013) The role of carrion in maintaining biodiversity and ecological processes in terrestrial ecosystems.171: 761–772. doi: 10.1007/s00442-012-2460-3BeguinJPothierDCôtéSD (2011) Deer browsing and soil disturbance induce cascading effects on plant communities: A multilevel path analysis.21: 439–451. doi: 10.1890/09-2100.1BressetteJWBeckHBeauchampVB (2012) Beyond the browse line: Complex cascade effects mediated by white-tailed deer.121: 1749–1760. doi: 10.1111/j.1600-0706.2011.20305.xBushERBueschingCDSladeEMMacdonaldDW (2012) Woodland recovery after suppression of deer: Cascade effects for small mammals, wood mice (Apodemus sylvaticus) and bank voles (Myodes glareolus).7: e31404. doi: 10.1371/journal.pone.0031404CamberfortYHanskiI (1991) Dung beetle population biology. In: HanskiICamberfortY (Eds) . Princeton University Press, Princeton, 36–50.CarpenterSR (1998) The need for large-scale experiments to assess and predict the response of ecosystems to perturbation. In: . Springer New York, New York, 287–312. doi: 10.1007/978-1-4612-1724-4_12CerdaYGrezAASimonettiJA (2015) The role of the understory on the abundance, movement and survival of Ceroglossus chilensis in pine plantations: An experimental test.19: 119–127. doi: 10.1007/s10841-015-9752-yCholletSMartinJ-L (2013) Declining woodland birds in North America: Should we blame Bambi?19: 481–483. doi: 10.1111/ddi.12003CôtéSDRooneyTPTremblayJ-PDussaultCWallerDM (2004) Ecological impacts of deer overabundance.35: 113–147. doi: 10.1146/annurev.ecolsys.35.021103.105725DekeirsschieterJVerheggenFLognayGHaubrugeE (2011) Large carrion beetles (Coleoptera, Silphidae) in Western Europe: A review.15: 435–447. http://www.scopus.com/inward/record.url?eid=2-s2.0-80052580984&partnerID=40&md5=c91f47becba8a56ffd60beacd96c1384DennisPYoungMRGordonIJ (1998) Distribution and abundance of small insects and arachnids in relation to structural heterogeneity of grazed, indigenous grasslands.23: 253–264. doi: 10.1046/j.1365-2311.1998.00135.xFosterCNBartonPSLindenmayerDB (2014) Effects of large native herbivores on other animals.51: 929–938. doi: 10.1111/1365-2664.12268FujitaAMaetoKKagawaYItoN (2008) Effects of forest fragmentation on species richness and composition of ground beetles (Coleoptera: Carabidae and Brachinidae) in urban landscapes.11: 39–48. doi: 10.1111/j.1479-8298.2007.00243.xGeospatial Information Authority of Japan (2015) . http://www.gsi.go.jp/kiban/ [accessed 9 December 2015]GillRMA (1992) A review of damage by mammals in north temperate forests: 3. Impact on trees and forests.65: 363–388. doi: 10.1093/forestry/65.4.363-aGómezJMGonzález-MegíasA (2002) Asymmetrical interactions between ungulates and phytophagous insects: Being different matters.83: 203–211. doi: 10.1890/0012-9658(2002)083[0203:AIBUAP]2.0.CO;2HoltCAFullerRJDolmanPM (2011) Breeding and post-breeding responses of woodland birds to modification of habitat structure by deer.144: 2151–2162. doi: 10.1016/j.biocon.2011.05.004HorsleySStoutSDeCalestaD (2003) White-tailed deer impact on the vegetation dynamics of a northern hardwood forest.13: 98–118. doi: 10.1890/1051-0761IidaTSogaMHiuraTKoikeS (2016) Life history traits predict insect species responses to large herbivore overabundance: a multitaxonomic approach.20: 295–304. doi: 10.1007/s10841-016-9866-xIkedaHHommaKKubotaK (2005) Biotic and abiotic factors affecting the structures of ground invertebrate communities in Japanese cedar dominant forests.8: 1–13. http://133.87.26.249/dspace/handle/2115/22186IkedaTTakahashiHYoshidaTIgotaHKajiK (2013) Evaluation of camera trap surveys for estimation of sika deer herd composition.38: 29–33. doi: 10.3106/041.038.0103JelaskaLSDurbes˘ićP (2009) Comparison of the body size and wing form of carabid species (Coleoptera: Carabidae) between isolated and continuous forest habitats.45: 327–338. doi: 10.1080/00379271.2009.10697618KajiKYajimaTIgarashiT (1991) Forage selection by deer introduced on Nakanoshima Island and its effect on the forest vegetation. In: MaruyamaNBobekBOnoYRegelinWBartosLRatcliffePR (Eds) . Japan Wildlife Research Center, Tokyo, 52–55.KandaNYokotaTShibataESatoH (2005) Diversity of dung-beetle community in declining Japanese subalpine forest caused by an increasing sika deer population.20: 135–141. doi: 10.1007/s11284-004-0033-6KevanPGBakerHG (1983) Insects as flower visitors and pollinators.28: 407–453. doi: 10.1146/annurev.en.28.010183.002203KoikeSSogaMNemotoYKozakaiC (2014) How are dung beetle species affected by deer population increases in a cool temperate forest ecosystem?293: 227–233. doi: 10.1111/jzo.12138LarsenTHWilliamsNMKremenC (2005) Extinction order and altered community structure rapidly disrupt ecosystem functioning.8: 538–547. doi: 10.1111/j.1461-0248.2005.00749.xMaguraTTóthmérészBLöveiGL (2006) Body size inequality of carabids along an urbanisation gradient.7: 472–482. doi: 10.1016/j.baae.2005.08.005MartinJ-LStocktonSAAllombertSGastonAJ (2010) Top-down and bottom-up consequences of unchecked ungulate browsing on plant and animal diversity in temperate forests: Lessons from a deer introduction.12: 353–371. doi: 10.1007/s10530-009-9628-8MartinTGArcesePScheerderN (2011) Browsing down our natural heritage: Deer impacts on vegetation structure and songbird populations across an island archipelago.144: 459–469. doi: 10.1016/j.biocon.2010.09.033McKinneyMLLockwoodJL (1999) Biotic homogenization: A few winners replacing many losers in the next mass extinction.14: 450–453. doi: 10.1016/S0169-5347(99)01679-1MelisCSundbyMAndersenRMoksnesAPedersenBRøskaftE (2007) The role of moose Alcesalces L. in boreal forest – the effect on ground beetles (Coleoptera, Carabidae) abundance and diversity.16: 1321–1335. doi: 10.1007/s10531-005-6230-yMelisCTeurlingsILinnellJDCAndersenRBordoniA (2004) Influence of a deer carcass on Coleopteran diversity in a Scandinavian boreal forest: A preliminary study.50: 146–149. doi: 10.1007/s10344-004-0051-2MiyakiMKajiK (2009) Shift to litterfall as year-round forage for sika deer after a population crash. In: McCulloughDRTakatsukiSKajiK (Eds) . Springer Japan, Tokyo, 171–180. doi: 10.1007/978-4-431-09429-6_12MysterudAAaserudRHansenLOÅkraKOlbergSAustrheimG (2010) Large herbivore grazing and invertebrates in an alpine ecosystem.11: 320–328. doi: 10.1016/j.baae.2010.02.009NicholsEGardnerTaPeresCaSpectorS (2009) Co-declining mammals and dung beetles: An impending ecological cascade.118: 481–487. doi: 10.1111/j.1600-0706.2009.17268.xNuttleTRistauTERoyoAA (2014) Long-term biological legacies of herbivore density in a landscape-scale experiment: Forest understoreys reflect past deer density treatments for at least 20 years.102: 221–228. doi: 10.1111/1365-2745.12175OksanenJBlanchetFGKindtRLegendrePMinchinPRO’HaraRBSimpsonGLSolymosPStevensMHHWagnerH (2015) . http://CRAN.R-project.org/package=vegan [accessed 5 November 2015]RainioJNiemeläJ (2003) Ground beetles (Coleoptera: Carabidae) as bioindicators.12: 487–506. doi: 10.1023/A:1022412617568R Core Team (2015) . R Foundation for Statistical Computing, Vienna.RobertsDW (2015) . http://CRAN.R-project.org/package=labdsv [accessed 19 November 2015]RooneyTP (2009) High white-tailed deer densities benefit graminoids and contribute to biotic homogenization of forest ground-layer vegetation.202: 103–111. doi: 10.1007/s11258-008-9489-8RooneyTPWallerDM (2003) Direct and indirect effects of white-tailed deer in forest ecosystems.181: 165–176. doi: 10.1016/S0378-1127(03)00130-0RouabahALasserre-JoulinFAmiaudBPlantureuxS (2014) Emergent effects of ground beetles size diversity on the strength of prey suppression.39: 47–57. doi: 10.1111/een.12064ScottMP (1998) The ecology and behavior of burying beetles.43: 595–618. doi: 10.1146/annurev.ento.43.1.595SpeightMRHunterMDWattAD (2008) . 2nd ed. Wiley-Blackwell, Oxford, 640 pp.StewartAJA (2001) The impact of deer on lowland woodland invertebrates: A review of the evidence and priorities for future research.74: 259–270. doi: 10.1093/forestry/74.3.259SuttiprapanPNakamuraH (2007) Species composition and seasonal abundance of carabid beetles by three sampling methods on the campus of the faculty of agriculture, Shinshu University.18: 83–90.SumwaysMJ (2007) (1st edn). Springer Netherlands, Heidelberg, 358 pp.SuominenONiemeläJMartikainenPNiemeläPKojolaI (2003) Impact of reindeer grazing on ground-dwelling Carabidae and Curculionidae assemblages in Lapland.26: 503–513. doi: 10.1034/j.1600-0587.2003.03445.xTakahashiHKajiK (2001) Fallen leaves and unpalatable plants as alternative foods for sika deer under food limitation.16: 257–262. doi: 10.1046/j.1440-1703.2001.00391.xTakatsukiS (2009) Effects of sika deer on vegetation in Japan: A review.142: 1922–1929. doi: 10.1016/j.biocon.2009.02.011TanentzapAJBazelyDRKohSTimciskaMHaggithEGCarletonTJCoomesDA (2011) Seeing the forest for the deer: Do reductions in deer-disturbance lead to forest recovery?144: 376–382. doi: 10.1016/j.biocon.2010.09.015UenoSKurosawaYSatoM (1985) Hoikusha publishing, Osaka, 514 pp.WheatallLNuttleTYergerE (2013) Indirect effects of pandemic deer overabundance inferred from caterpillar-host relations.27: 1107–1116. doi: 10.1111/cobi.12077WhiteMA (2012) Long-term effects of deer browsing: Composition, structure and productivity in a northeastern Minnesota old-growth forest.269: 222–228. doi: 10.1016/j.foreco.2011.12.043YamadaHTakatsukiS (2015) Effects of deer grazing on vegetation and ground-dwelling insects in a larch forest in Okutama, Western Tokyo.2015: 1–9. doi: 10.1155/2015/687506Supplementary materials10.3897/zookeys.625.9116.suppl1903706383F3883-43F6-51D1-9FE6-193E606D357E
Table 1
Data table
A list of beetle species observed in our study.
https://binary.pensoft.net/file/107253This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Taichi Iida, Masashi Soga, Shinsuke Koike10.3897/zookeys.625.9116.suppl29037085C2E7CEB-AFF4-5B6B-9450-B147339A532F
Table 2
Data table
Estimated and observed values of beetle abundance at the island and lakeshore sites.
https://binary.pensoft.net/file/107254This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Taichi Iida, Masashi Soga, Shinsuke Koike10.3897/zookeys.625.9116.suppl3903710D40EA483-44FD-573D-BD6A-1616D13A97C6
Table 3
Data table
Estimated and observed values of beetle species richness at the island and lakeshore sites.
https://binary.pensoft.net/file/107255This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Taichi Iida, Masashi Soga, Shinsuke Koike10.3897/zookeys.625.9116.suppl4903712841E242C-F8FC-5205-A484-1A37FF57F807
Table 4
Data table
Estimated and observed values of beetle diversity (Shannon-Wiener diversity index) at the island and lakeshore sites.
https://binary.pensoft.net/file/107256This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Taichi Iida, Masashi Soga, Shinsuke Koike10.3897/zookeys.625.9116.suppl5903714BB63E89C-DE8E-5AA4-B8AB-E4632EC87C76
Figure 1
Figure
Species accumulation curves for carabid, carrion and dung beetles at the island (blue lines) and lakeshore sites (red lines).
https://binary.pensoft.net/file/107257This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Taichi Iida, Masashi Soga, Shinsuke Koike