# Abundance, not diversity, of host beetle communities determines abundance and diversity of parasitoids in deadwood

Sebastian Vogel, Andreas Prinzing, Heinz Bussler, Jorg Mueller, Stefan Schmidt, Simon Thorn

► To cite this version:

Sebastian Vogel, Andreas Prinzing, Heinz Bussler, Jorg Mueller, Stefan Schmidt, et al.. Abundance, not diversity, of host beetle communities determines abundance and diversity of parasitoids in deadwood. Ecology and Evolution, 2021, 11 (11), pp.6881-6888. 10.1002/ece3.7535. hal-03216899

HAL Id: hal-03216899

https://hal.science/hal-03216899

Submitted on 26 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Abundance, not diversity, of host beetle communities determines abundance and diversity of parasitoids in deadwood

Sebastian Vogel<fnref n="1" /> ID

Andreas Prinzing<fnref n="2" />

Heinz Bußler<fnref n="1" />

Jörg Müller<fnref n="1" /> ID

Stefan Schmidt<fnref n="3" />

Simon Thorn<fnref n="1" /> ID

### Correspondence

Simon Thorn, Department of Animal Ecology and Tropical Biology (Zoology III), Field Station Fabrikschleichach, Julius Maximilians University Würzburg, Glashüttenstraße 5, 96181 Rauhenebrach, Germany.
Email: simon@thornonline.de

### Funding information

Scholarship Program of the German Federal Environmental Foundation, Grant/Award Number: 20016/466; Federal Ministry of Education and Research; Genome Canada

### Abstract

Most parasites and parasitoids are adapted to overcome defense mechanisms of their specific hosts and hence colonize a narrow range of host species. Accordingly, an increase in host functional or phylogenetic dissimilarity is expected to increase the species diversity of parasitoids. However, the local diversity of parasitoids may be driven by the accessibility and detectability of hosts, both increasing with increasing host abundance. Yet, the relative importance of these two mechanisms remains unclear. We parallelly reared communities of saproxylic beetle as potential hosts and associated parasitoid Hymenoptera from experimentally felled trees. The dissimilarity of beetle communities was inferred from distances in seven functional traits and from their evolutionary ancestry. We tested the effect of host abundance, species richness, functional, and phylogenetic dissimilarities on the abundance, species richness, and Shannon diversity of parasitoids. Our results showed an increase of abundance, species richness, and Shannon diversity of parasitoids with increasing beetle abundance. Additionally, abundance of parasitoids increased with increasing species richness of beetles. However, functional and phylogenetic dissimilarity showed no effect on the diversity of parasitoids. Our results suggest that the local diversity of parasitoids, of ephemeral and hidden resources like saproxylic beetles, is highest when resources are abundant and thereby detectable and accessible. Hence, in some cases, resources do not need to be diverse to promote parasitoid diversity.

### KEYWORDS

barcoding, deadwood, experiment, host-parasitoid interaction, natural enemy, specialization

## 1 | INTRODUCTION

Most species worldwide are natural enemies of other species, as exemplified by parasitoids and parasites. Those species tend to be neither specialized on a single host species nor completely generalist,

colonizing a set of similar and often closely related hosts (Strong et al., 1984). Natural enemies have hence adapted to overcome the chemical, morphological, physiological, and immunological defense characteristics of their host species, which can be phylogenetically conserved (Gross, 1993; Strand & Pech, 1995). As a result, hosts

<footnote n="1">¹Department of Animal Ecology and Tropical Biology (Zoology III), Field Station Fabrikschleichach, Julius Maximilians University Würzburg, Rauhenebrach, Germany</footnote>

<footnote n="2">²Research Unit “Ecosystèmes Biodiversité, Evolution” («UMR 6553»), Centre National de la Recherche Scientifique, University Rennes 1, Rennes, France</footnote>

<footnote n="3">³SNSB-Zoologische Staatssammlung München, Munich, Germany</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.</footnote>

<footnote>© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.</footnote>
with different functional characteristics or evolutionary ancestry
tend to be exploited by different enemies (Frank, 1993). Hence, in-
creasing phylogenetic or functional dissimilarity of host communi-
ties is thought to increase the species diversity of their enemies in
a community (Schuler et al., 2015; Vialatte et al., 2010). This might
be the case in systems in which functionally different hosts require
fundamentally different adaptations from their enemies and the
performance of an enemy to colonize one host species is traded off
against its performance on colonizing another (Egas et al., 2005;
Straub et al., 2011).

Enemy diversity may also depend on the detectability and
accessibility of hosts: If hosts are hidden or sheltered, only the
most specialized enemy species will be able to detect or use them
(Price, 1973). Host detectability and accessibility, in turn, may in-
crease with increasing host abundance. Host abundance increases
host detectability by increasing olfactory, acoustical, or optical sig-
nals available for enemies (Yguel et al., 2014). Host abundance may
also increase among-host competition forcing to leave the enemy-
free space, that is, microsites or periods where or when hosts are
inaccessible to enemies (Aukema & Raffa, 2000). Also, increased
competition may require increased investment into competitive-
ness, binding energy that could otherwise be invested into defense
against enemies (Bashey, 2015). In sum, low host abundance may
reduce the diversity of natural enemies to a narrow subset of enemy
species.

Parasitoids constitute several different taxonomic groups and
are a key component of terrestrial ecosystems (Heraty, 2009). It
is estimated that 10%-20% of all insect species belong to parasit-
oid Hymenoptera, representing the insect group with the highest
species diversity worldwide (Forbes et al., 2018; Gaston, 1991).
Comparable to parasites, the majority of parasitoid Hymenoptera is
specialized on a certain set of hosts (Forbes et al., 2018), although
many host species are so far unknown. As a result, a given parasitoid
Hymenoptera species is typically most abundant only on one or a
few focal host species.

Parasitoid Hymenoptera are one of the most important nat-
ural enemies of saproxylic beetle communities in forests, as
mainly demonstrated for bark beetles (Curculionidae, Scolytinae)
(Wegensteiner et al., 2015). Saproxylic beetles differ fundamen-
tally in functional properties, such as in their preference for certain
host tree species, diameters of deadwood, or their dependence on
specific microclimatic conditions (Seibold et al., 2015). Moreover,
communities of saproxylic beetles change with succession
during deadwood decomposition (Parisi et al., 2018). Parasitoid
Hymenoptera mainly attack the eggs, larvae, and pupae of sap-
roxylic beetles and exhibit a clear preference for one develop-
mental stage over the others (Strand & Pech, 1995; Wegensteiner
et al., 2015). Many of these and other functional properties are
phylogenetically conserved (Seibold et al., 2015). The functional
or phylogenetic dissimilarity of saproxylic species might hence
increase species diversity of parasitoid enemies. However, given
that most saproxylic species are extremely ephemeral and hidden,

the diversity of their parasitoids might also increase with saprox-
ylic abundance.

We investigated communities of saproxylic beetle hosts and their
associated parasitoid Hymenoptera. Both groups of insects were
reared from experimentally felled trees during the early succession
of deadwood. For each tree, the species diversity of saproxylic beet-
tles and parasitoid Hymenoptera emerging from those beetles were
quantified. In addition, the saproxylic beetle communities of each
tree erized according to their similarities across seven functional
traits and their phylogenetic similarity, to measure functional and
phylogenetic dissimilarities. Our aim was to test whether abundance
or functional and phylogenetic dissimilarity of host communities de-
termines the abundance and species richness of associated parasit-
oid Hymenoptera.

2 | MATERIALS AND METHODS

## 2.1 | Study area and experimental design

Our study was conducted in the Bavarian Forest National Park, located in southeastern Germany. The park consists of approximately 24,850 ha of mountainous forests at elevations between 600 and 1,460 m a.s.l. Depending on the altitude, the annual average temperature ranges from 3.8 to 5.8°C. Yearly precipitation varies between 1,200 and 1,800 mm (Bässler et al., 2010). Forest stands in higher elevations of the park area are naturally dominated by Norway spruce (*Picea abies*). Within the last two decades, extensive waves of natural disturbances have generated highly diverse deadwood structures, contributing substantially to partially more than 300 m³ of deadwood per hectare (Thorn et al., 2017).

Observational studies do naturally depend on the occurrence
and spatial distribution of natural disturbances, which typically
do not match the requirements of a standardized scientific study
design (Lindenmayer et al., 2010). Also, often only a single host
species is dominant, and results might be idiosyncratic to that
species. Hence, we created artificial wind throws and applied bark
treatments to create variance in the community composition of
saproxylic beetles, that is, to avoid the dominance of the European
bark beetle (*Ips typographus*). *Ips typographus* is the most import-
ant pest species of mature spruce stands throughout Eurasia with
a preference for weakened or freshly dead Norway spruce trees
(Wermelinger, 2004).

Our field experiment was established in April 2013. The design
consisted of 12 artificial wind throws (plots), each composed of three
pulled down mature Norway spruce trees with similar physical attri-
butes (Figure 1). Two trees per plot were uprooted and debranched,
and their root plates were cut off. One tree was bark-scratched (dis-
ruption of the phloem every 3 cm), and a second tree was debarked
/removal of all phloem). The remaining tree served as control. The
minimum distance between plots was 200 m (for more details of the
experimental design, see Thorn, Bässler, et al. (2016), Thorn, Bußler,
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

**FIGURE 1** (a) Different mechanical bark treatments used to alter abundance and species composition of saproxylic beetle communities. On each plot, one tree was fully debarked, one tree was bark-scratched (disruption of the phloem every 3 cm), and the third tree served as the control. (b) Stem emergence trap used for sampling of saproxylic beetles and parasitoid Hymenoptera. (c) The European spruce bark beetle (Ips typographus) and a parasitoid Hymenoptera (Superfamily Chalcidoidea) as representatives for the trapped species and analyzed communities (pictures were used after the Creative Commons license (CC BY-SA-2.0 and CC BY-3.0); copyright by U. Schmidt and M.A. Broussard)

et al. (2016). Within the conditions provided in our experiment, communities of beetles and parasitoids assembled naturally.

## 2.2 | Arthropod sampling

Communities of saproxylic beetles and parasitoid Hymenoptera were sampled by using one stem emergence trap per tree, which were established from April to September 2014. Stem emergence traps are covering a given tree section and are only capturing individuals directly emerging from this section, representing the local community more complete the, for example, flight interception traps (Sverdrup-Thygeson & Birkemoe, 2009). To preserve species material for further steps, we used 90% ethanol and emptied traps monthly. All trapped beetles were identified to the species level according to Freude et al. (1963–1984) and classified as saproxylic according to Schmidl & Bußler (2004). To identify species of parasitoid Hymenoptera, one leg of each specimen was removed and submitted to the Canadian Centre for DNA Barcoding for DNA sequence analysis. DNA was extracted according to Ivanova et al. (2006), and standardized primer sets were used to amplify and sequence the 658-bp barcode region (Folmer et al., 1994; Hebert et al., 2003). The sequences were aligned with those in the Barcode of Life data system (Ratnasingham & Hebert, 2007). If the identification to the species level was not possible, barcoding data were used as molecular operational taxonomic units (MOTUs), which are represented by barcode index numbers and closely approximate species-level identifications (Ratnasingham & Hebert, 2013). Hereafter, MOTUs are included in

"species" for simplicity. Single MOTUs were defined by a 97% sequence similarity. Species or MOTUs with sequences that matched those of parasitoids certainly not using coleopteran hosts according to Yu et al. (2016) and Noyes (2019) were excluded. The beetle data used in our paper are available online in Chao et al. (2019).

## 2.3 | Host phylogeny and traits

We used the phylogenetic tree of European saproxylic beetles provided by Seibold et al. (2015) according to the approach of Kuhn et al. (2011). The latter allows the input of partially resolved trees with known topology and node ages as constraints and applies a Markov chain Monte Carlo algorithm to permute polytomies using a constant-rate birth–death model. The phylogenetic tree was calibrated using 25 calibration points obtained from fossil records.

Functional dissimilarity was analyzed using the resource-related and morphological traits of saproxylic beetles described by Gossner et al. (2013), Thorn et al. (2014), and Seibold et al. (2015). These databases represent the most comprehensive information about saproxylic beetle functional traits in Central Europe. Information on the mean diameter niche position (<15, 15–35, 35–70, >70 cm), the decay niche position (alive, freshly dead, initiated, advanced decomposition, extremely decomposed), and the canopy cover niche (open, semi-open, closed) was included, together with information on the preference of saproxylic beetle species for host trees (coniferous, broadleaved, both types of trees), the microhabitat guilds of larvae (wood and bark, cavities,
fungi), the mean body size of single species, and the larval feeding type (detritivorous, mycetophagous, xylophagous, zoophagous). For the degree of phylogenetic correlation in functional traits, see Doerfler et al. (2020).

## 2.4 | Dissimilarity of host communities

All statistical analyses were carried out in R version 4.0.3 (www.r-project.org). The phylogenetic tree of European saproxylic beetles was used to quantify phylogenetic dissimilarity based on the mean nearest taxon distance (MNTD), reflecting the mean distance separating each species from its closest relative in the same community. MNTD was additionally used to quantify the similarities among co-occurring species that may influence colonization by parasitoids. A value of MNTD less than zero thereby represents a clustered (i.e., similar) community of hosts, while a value greater than zero represents an overdispersed (i.e., dissimilar) community. To calculate MNTD, we first calculated the cophenetic distances between co-occurring species using the function "cophenetic" from the R package "stats." Because MNTD values can decrease as the number of species in a community increases, an abundance-weighted null model was applied to compare the observed MNTD values with those of 999 randomly generated communities containing the same number of species by using the function "ses.mntd" from the R package "picante" (Cavender-Bares et al., 2006; Kembel et al., 2010). This procedure resulted in a standardized effect size of MNTD that reflected the mean MNTD, that is, the phylogenetic dissimilarity of a beetle community standardized by species numbers. Standardized effect sizes of functional and phylogenetic diversities were only weakly correlated (Pearson's *r*: 0.29).

The functional dissimilarity of saproxylic beetle communities was calculated using the same procedure as described above, except that the distance matrices of traits were calculated as Euclidean (numeric variables) and Gower (non-numeric or in combination with numeric variables) distances using the function "daisy" from the R package "cluster" (Gower, 1971).

## 2.5 | Statistical analyses

The effects of abundance, species richness, phylogenetic dissimilarity, and functional dissimilarity of saproxylic beetle communities were tested on the abundance, species richness, and Shannon diversity of parasitoids as response variables. We standardized the number of sampled saproxylic beetle species to 0.95 sample completeness using a rarefaction/extrapolation framework according to Hsieh et al. (2016), implemented in the R package "iNext." We estimated the sample coverage of saproxylic beetles and parasitoid Hymenoptera by the same approach.

The effect of saproxylic beetle communities on the abundance of parasitoids was tested by applying a generalized linear mixed model with poisson error distribution by using the R package "lme4"

(Bates et al., 2015). This model included the abundance of parasitoids as response variable and the log-transformed abundance, species richness, functional dissimilarity, and phylogenetic dissimilarity of the respective host beetle communities as predictor variables. Furthermore, we included the plot as random effect to account for the nested study design (Bolker et al., 2009) and the sample id to account for possible Poisson overdispersion (Elston et al., 2001). We used the same model formula to model species numbers and additionally included the log-transformed abundance of parasitoid Hymenoptera to the model. Shannon diversity of parasitoids was modeled using the same model formula as for species numbers but a Gaussian error distribution.

## 3 | RESULTS

Our dataset included 15,516 individuals of saproxylic beetles from 106 species and 146 individuals from 44 species of parasitoids (Table 1), corresponding to a 0.93% rate of successful parasitization. The mean species number per log was 23.61 ± 7.5 species for saproxylic beetles and 2.64 ± 2.6 for parasitoids. Species accumulation curves revealed a high sample coverage of both, saproxylic beetles and parasitoid Hymenoptera (Figure 2).

Saproxylic beetles were dominated by bark beetles (Curculionidae, Scolytinae), encompassing 10,869 individuals from 18 species. The most abundant beetle species, *Trypodendron lineatum*, was represented by 7,496 trapped individuals. *Ips typographus* was trapped by 824 (control trees), 180 (bark-scratched trees), and 19 (debarked trees) individuals, respectively. Parasitoids were recorded from 11 families, dominated by Pteromalidae (40 individuals, 10 species), Platygastridae (36 individuals, 10 species), and Braconidae (33 individuals, 6 species). Many of the parasitoids recorded in our study were rare: 36 species of parasitoids with ≤3 individuals contributed 66 of the 146 sampled individuals (Table 1).

Abundance of parasitoids increased with increasing abundance and species richness of saproxylic beetles (Figure 3a). Species numbers and Shannon diversity of parasitoid Hymenoptera strongly increased with increasing abundance of parasitoid Hymenoptera (Figure 3b,c). Shannon diversity of parasitoid Hymenoptera increased with increasing host abundance (Figure 3c). We did not find any significant effect of functional or phylogenetic dissimilarity on any of our response variables (Figure 3).

## 4 | DISCUSSION

Our results indicate that abundance and species richness in communities of parasitoid Hymenoptera increase with increasing abundance of hosts, whereas functional and phylogenetic dissimilarities of host communities were of minor importance. Furthermore, species numbers of parasitoids were strongly promoted by their abundance, which in turn increases with species richness of hosts. Thereby, our results do not support that functional or phylogenetic
<table><caption>**TABLE 1** Species of parasitoid Hymenoptera recorded in our study with their respective abundances</caption><thead><tr><th>ID</th><th>Scientific name</th><th>Abundance</th></tr></thead><tbody><tr><td colspan="3">Braconidae (6 species/33 individuals)</td></tr><tr><td>BOLD:AAU9839</td><td>Chelonus sp.</td><td>1</td></tr><tr><td>BOLD:ACM7216</td><td>Ropalophorus clavicornis</td><td>1</td></tr><tr><td>BOLD:ACM7563</td><td>Dendrosoter middendorffii</td><td>5</td></tr><tr><td>BOLD:ACQ8673</td><td>Brachistinae</td><td>1</td></tr><tr><td>BOLD:ACQ9535</td><td>Cosmophorus regius</td><td>2</td></tr><tr><td>BOLD:ACQ9771</td><td>Blacus sp.</td><td>23</td></tr><tr><td colspan="3">Ceraphronidae (3 species/5 individuals)</td></tr><tr><td>BOLD:ACF6025</td><td>Aphanogmus sp.</td><td>3</td></tr><tr><td>BOLD:ACG4508</td><td>Aphanogmus sp.</td><td>1</td></tr><tr><td>Ceraph_01</td><td>Aphanogmus sp.</td><td>1</td></tr><tr><td colspan="3">Diapriidae (2 species/2 individuals)</td></tr><tr><td>BOLD:ACG4055</td><td>Pantoclis sp.</td><td>1</td></tr><tr><td>Diapri_01</td><td>Trichopria sp.</td><td>1</td></tr><tr><td colspan="3">Eulophidae (4 species/7 individuals)</td></tr><tr><td>BOLD:ACQ8898</td><td>Necremnus croton</td><td>2</td></tr><tr><td>BOLD:ACQ9006</td><td>Necremnus leucarthros</td><td>3</td></tr><tr><td>Euloph_01</td><td></td><td>1</td></tr><tr><td>Euloph_02</td><td></td><td>1</td></tr><tr><td colspan="3">Eupelmidae (1 species/2 individuals)</td></tr><tr><td>Eupelm_01</td><td></td><td>2</td></tr><tr><td colspan="3">Eurytomidae (1 species/3 individuals)</td></tr><tr><td>BOLD:ACM7745</td><td>Eurotoma arctica/afra</td><td>3</td></tr><tr><td colspan="3">Figitidae (1 species/1 individual)</td></tr><tr><td>BOLD:ACQ9714</td><td>Eucoilinae</td><td>1</td></tr><tr><td colspan="3">Ichneumonidae (4 species/13 individuals)</td></tr><tr><td>BOLD:ABU6543</td><td>Enclisis vindex</td><td>2</td></tr><tr><td>BOLD:ACQ9146</td><td>Phrudus monilicornis</td><td>3</td></tr><tr><td>BOLD:ACR0681</td><td>Phradis sp.</td><td>1</td></tr><tr><td>BOLD:ACR0964</td><td>Rhimphoctona teredo</td><td>7</td></tr><tr><td colspan="3">Platygastridae (10 species/36 individuals)</td></tr><tr><td>BOLD:AAN8098</td><td>Telenomus sp.</td><td>9</td></tr><tr><td>BOLD:ACC2809</td><td>Sclexioninae</td><td>2</td></tr><tr><td>BOLD:ACF7380</td><td></td><td>1</td></tr><tr><td>BOLD:ACF9487</td><td></td><td>2</td></tr><tr><td>BOLD:ACI4334</td><td>Telenomus sp.</td><td>10</td></tr><tr><td>BOLD:ACI4527</td><td>Telenomus sp.</td><td>1</td></tr><tr><td>BOLD:ACI9091</td><td>Telenomus sp.</td><td>4</td></tr><tr><td>BOLD:ACI9128</td><td>Telenomus sp.</td><td>3</td></tr><tr><td>BOLD:ACR1479</td><td></td><td>2</td></tr><tr><td>BOLD:ACR1888</td><td></td><td>2</td></tr><tr><td colspan="3">Proctotrupidae (2 species/4 individuals)</td></tr><tr><td>BOLD:ACR1488</td><td></td><td>2</td></tr></tbody></table>

<table><caption>**TABLE 1** (Continued)</caption><thead><tr><th>ID</th><th>Scientific name</th><th>Abundance</th></tr></thead><tbody><tr><td colspan="3">BOLD:ACR1985</td></tr><tr><td colspan="3">Pteromalidae (10 species/40 individuals)</td></tr><tr><td>BOLD:AAN8215</td><td>Pteromalus sp.</td><td>2</td></tr><tr><td>BOLD:AAZ7417</td><td>Mesopolobus gemellus</td><td>3</td></tr><tr><td>BOLD:ACA9177</td><td></td><td>1</td></tr><tr><td>BOLD:ACM7334</td><td>Roptrocerus mirus</td><td>15</td></tr><tr><td>BOLD:ACM7652</td><td></td><td>2</td></tr><tr><td>BOLD:ACQ8466</td><td>Dinotiscus eupterus</td><td>3</td></tr><tr><td>BOLD:ACQ8826</td><td>Holcaeus compressus</td><td>7</td></tr><tr><td>BOLD:ACQ9876</td><td>Perniphora robusta</td><td>3</td></tr><tr><td>Pterom_01</td><td></td><td>3</td></tr><tr><td>Pterom_03</td><td></td><td>1</td></tr></tbody></table>

dissimilarity of hosts begets the diversity of parasitoids within dead-wood objects, but rather point toward an abundance-driven system.

The increasing abundance and species richness of parasitoids with increasing abundance of their saproxylic beetle hosts are in line with the »more individuals hypothesis«, an extension of the »species-energy theory« (Clarke & Gaston, 2006; Schuler et al., 2015). The »more individuals hypothesis« predicts an increase in the abundance and consequently in species richness of a given species group (i.e., parasitoids) in response to an increase in the availability of resources (i.e., hosts). The »more individuals hypothesis« is valid as long as species use different niches and do not rely on the same resources (MacArthur & MacArthur, 1961; Tews et al., 2004). In our study system, an increasing abundance could increase the detectability (Vinson, 1998) and accessibility (Price, 1973) of saproxylic beetle hosts.

In our study system, host accessibility may be limited, given the short spatiotemporal window during which a specific developmental stage of a host exists. Furthermore, parasitoids may be limited to a specific set of host species and do not utilize patches with extremely low host abundance (Hassell, 2000; Murdoch, 1977).

High host abundances may force some hosts into enemy-exposed microenvironments. In our study system, this could apply to bark beetles, which reached highest abundances of potential hosts and then might be forced to occupy parts of the tree trunk which are more easier to access to parasitoids (Aukema & Raffa, 2000).

Higher abundances of hosts might also increase the detectability, that is the olfactory signal to parasitoids, since saproxylic beetles, especially bark beetles, emit sex pheromones (Vega & Hofstetter, 2015). Furthermore, at extremely high abundances, bark beetles can emit anti-aggregation pheromones (Sun et al., 2006). Such pheromones might additionally increase host detectability at high abundances. Overall, abundance of bark beetles might render them distinctly more detectable or accessible to a larger number of parasitoid species and thereby increase their species diversity.

We did not find any effect of functional or phylogenetic dissimilarity of hosts on parasitoid communities (Figure 3). Functional
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Sample coverage of saproxylic beetles and parasitoid Hymenoptera based on rarefaction/extrapolation up to twice the actual sampling effort. Note that the number of sampled individuals has been transformed to percentage for better comparability</figcaption></figure>

dissimilarity might to some degree reflect functional identity. If there is one dominant functional group, such as in our case bark beetles, then a greater functional similarity may reflect the dominance of that group. Specifically, our samples were composed of species colonizing mainly the phloem of recently killed Norway spruce. Specialization on particular host groups may reflect constraints on, for example, ovipositor length. For example, bark thickness limits parasitization by *Caenopachys hartigii* (Braconidae, Doryctinae) (Mancini et al., 2003). Thus, host species that do not feed directly under the bark or other superficial parts of deadwood might not be reachable by this enemy. A highly abundant bark beetle species in our study was *T. lineatum*, which create galleries into sapwood, where it might be difficult to access for parasitoids. Indeed, only few species of parasitoids enter bark beetle galleries for oviposition (Vega & Hofstetter, 2015).

It is possible that the effect of functional identity overlays and partly conceals the effect of functional dissimilarity, explaining why we did not find any effect of functional dissimilarity. Despite we used the most comprehensive database of saproxylic beetles, we are still lacking information about, for example, temporal occurrence of saproxylic beetles. An extension of our trait database by such information may hence change the relative importance of functional dissimilarity. Phylogenetic dissimilarity did not show a significant effect on parasitoid abundance (Figure 3). One of the reasons for the lacking effect of phylogenetic dissimilarity might be that closely related beetle hosts occur in high abundances, such as bark beetles in our study system. This is ultimately caused by the fact that host beetle communities cannot be established in a fully experimental design, that is, covering a range of phylogenetic dissimilarity, due to their host tree specializations. Moreover, for a set of locally coexisting species, phylogenetic dissimilarity may not necessarily reflect the dissimilarity of functional traits (Prinzing et al., 2008).

Overall, our study suggests that resource diversity does not mandatory begets the diversity of local enemy communities. This intuitive hypothesis implicitly assumes that any given type of resources is detectable and accessible to the consumers. In that case, a diverse set of resources can be used by a diverse set of consumers.

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**FIGURE 3** Relative effect strengths (z-values) of predictors on (a) abundance of parasitoid Hymenoptera, (b) species numbers of parasitoid Hymenoptera, and (c) Shannon diversity of parasitoid Hymenoptera. Gray shading indicates the range of nonsignificant z-values, significant effects depicted by black dots</figcaption></figure>

However, if resources are hosts that are ephemeral and hidden, they may easily remain undetected by or inaccessible to many potential parasitoids. Only at high host abundance, detection and access may become possible for many parasitoid species, leading to an increase of parasitoid diversity with host abundance.

## ACKNOWLEDGMENTS

SV was supported by the Scholarship Program of the German Federal Environmental Foundation (fellowship number 20016/466). DNA barcode sequences were obtained as part of the German Barcode of Life Project, funded by the Federal Ministry of Education and Research, with financial support by Genome Canada (Ontario Genomics Institute) in the framework of the iBOL program, WG 1.9.

## CONFLICT OF INTEREST

The authors declare no conflict of interest.

## AUTHOR CONTRIBUTIONS

**Sebastian Vogel:** Conceptualization (equal); Formal analysis (equal); Investigation (equal); Writing—original draft (lead); Writing—review & editing (equal). **Andreas Prinzing:** Conceptualization (equal);
Writing-review & editing (equal). Heinz Bußler: Data curation (equal); Resources (equal). Jörg Müller: Resources (equal); Writing-review & editing (equal). Stefan Schmidt: Data curation (equal); Investigation (equal); Resources (equal). Simon Thorn: Conceptualization (equal); Formal analysis (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing-original draft (equal); Writing-review & editing (equal).

## DATA AVAILABILITY STATEMENT

All specimen data are accessible in the Barcode of Life Database (BOLD, http://www.boldsystems.org) as a single citable dataset (https://dx.doi.org/10.5883/DS-GBBWA). The data include collecting locality, geographic coordinates, elevation, collector, one or more digital images, identifier, and voucher depository.

## ORCID

Sebastian Vogel

<figure><img src="image_8.png" /></figure>

https://orcid.org/0000-0002-3500-3960

Jörg Müller

<figure><img src="image_9.png" /></figure>

https://orcid.org/0000-0002-1409-1586

Simon Thorn

<figure><img src="image_10.png" /></figure>

https://orcid.org/0000-0002-3062-3060

## REFERENCES

Aukema, B. H., & Raffa, K. F. (2000). Chemically mediated predator-free space: Herbivores can synergize intraspecific communication without increasing risk of predation. *Journal of Chemical Ecology*, 26(8), 1923-1939. https://doi.org/10.1023/A:1005505009408

Bashey, F. (2015). Within-host competitive interactions as a mechanism for the maintenance of parasite diversity. *Philosophical Transactions of the Royal Society B: Biological Sciences*, 370(1675), 20140301. https://doi.org/10.1098/rstb.2014.0301

Bässler, C., Müller, J., & Dzioc, F. (2010). Detection of climate-sensitive zones and identification of climate change indicators: A case study from the Bavarian Forest National Park. *Folia Geobotanica*, 45(2), 163-182. https://doi.org/10.1007/s12224-010-9059-4

Bates, D., Mächler, M., Bolker, B., & Walker, S. (2015). Fitting Linear Mixed-Effects Models using lme4. *Journal of Statistical Software*, 67(1), 1-48. https://doi.org/10.18637/jss.v067.i01

Bolker, B. M., Brooks, M. E., Clark, C. J., Geange, S. W., Poulsen, J. R., Stevens, M. H. H., & White, J.-S.-S. (2009). Generalized linear mixed models: A practical guide for ecology and evolution. *Trends in Ecology & Evolution*, 24(3), 127-135. https://doi.org/10.1016/j.tree.2008.10.008

Cavender-Bares, J., Keen, A., & Miles, B. (2006). Phylogenetic structure of Floridian plant communities depends on taxonomic and spatial scale. *Ecology*, 87(7), 109-122. https://doi.org/10.1890/0012-9658(2006)87[109:PSOFPC]2.0.CO;2

Chao, A., Chiu, C.-H., Villéger, S., Sun, I.-F., Thorn, S., Lin, Y.-C., Chiang, J.-M., & Sherwin, W. B. (2019). An attribute-diversity approach to functional diversity, functional beta diversity, and related (dis)similarity measures. *Ecological Monographs*, 89(2), e01343. https://doi.org/10.1002/ecm.1343

Clarke, A., & Gaston, K. J. (2006). Climate, energy and diversity. *Proceedings of the Royal Society B: Biological Sciences*, 273(1599), 2257-2266. https://doi.org/10.1098/rspb.2006.3545

Doerfler, I., Cadotte, M. W., Weisser, W. W., Müller, J., Gossner, M. M., Heibl, C., Bässler, C., Thorn, S., & Seibold, S. (2020). Restoration-oriented forest management affects community assembly patterns of deadwood-dependent organisms. *Journal of Applied Ecology*, 57(12), 2429-2440. https://doi.org/10.1111/1365-2664.13741

Egas, M., Sabelis, M. W., & Dieckmann, U. (2005). Evolution of specialization and ecological character displacement of herbivores along

a gradient of plant quality. *Evolution*, 59(3), 507-520. https://doi.org/10.1111/j.0014-3820.2005.tb01011.x

Elston, D. A., Moss, R., Boulinier, T., Arrowsmith, C., & Lambin, X. (2001). Analysis of aggregation, a worked example: Numbers of ticks on red grouse chicks. *Parasitology*, 122, 563-569. https://doi.org/10.1017/S0031182001007740

Folmer, O., Black, M., Hoeh, W., Lutz, R., & Vrijenhoek, R. (1994). DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. *Molecular Marine Biology and Biotechnology*, 3(5), 294-299. https://doi.org/10.1371/journal.pone.0013102

Forbes, A. A., Bagley, R. K., Beer, M. A., Hippee, A. C., & Widmayer, H. A. (2018). Quantifying the unquantifiable: Why Hymenoptera, not Coleoptera, is the most speciose animal order. *BMC Ecology*, 18(1), 21. https://doi.org/10.1186/s12898-018-0176-x

Frank, S. A. (1993). Evolution of host-parasite diversity. *Evolution*, 47(6), 1721-1732.

Freude, H., Harde, K., & Lohse, G. A. (1963-1984). *Die Käfer Mitteleuropas*. Goecke und Evers.

Gaston, K. J. (1991). The magnitude of global insect species richness. *Conservation Biology*, 5(3), 283-296. https://doi.org/10.1111/j.1523-1739.1991.tb00140.x

Gossner, M. M., Lachat, T., Brunet, J., Isacsson, G., Bouget, C., Brustel, H., Brandl, R., Weisser, W. W., & Müller, J. (2013). Current near-tot nature forest management effects on functional trait composition of saproxylic beetles in beech forests. *Conservation Biology*, 27(3), 605-614. https://doi.org/10.1111/cobi.12023

Gower, J. C. (1971). A general coefficient of similarity and some of its properties. *Biometrics*, 27(4), 857-871. https://doi.org/10.2307/2528823

Gross, P. (1993). Insect behavioral and morphological defenses against parasitoids. *Annual Review of Entomology*, 38(1), 251-273. https://doi.org/10.1146/annurev.ento.38.1.251

Hassell, M. P. (2000). *The spatial and temporal dynamics of host-parasitoid interactions*. Oxford University Press.

Hebert, P. D. N., Cywinska, A., Ball, S. L., & DeWaard, J. R. (2003). Biological identifications through DNA barcodes. *Proceedings of the Royal Society of London. Series B: Biological Sciences*, 270(1512), 313-321. https://doi.org/10.1098/rspb.2002.2218

Heraty, J. (2009). Parasitoid biodiversity and insect pest management. In R. G. Footit, & P. H. Adler (Eds.), *Insect biodiversity: Science and society, second*. Wiley-Blackwell. https://doi.org/10.1002/9781444308211.ch19

Hsieh, T. C., Ma, K. H., & Chao, A. (2016). iNEXT: An R package for rarefaction and extrapolation of species diversity (Hill numbers). *Methods in Ecology and Evolution*, 7(12), 1451-1456. https://doi.org/10.1111/2041-210X.12613

Ivanova, N. V., Dewaard, J. R., & Hebert, P. D. N. (2006). An inexpensive, automation-friendly protocol for recovering high-quality DNA. *Molecular Ecology Notes*, 6(4), 998-1002. https://doi.org/10.1111/j.1471-8286.2006.01428.x

Kembel, S. W., Cowan, P. D., Helmus, M. R., Cornwell, W. K., Morlon, H., Ackerly, D. D., Blomberg, S. P., & Webb, C. O. (2010). Picante: R tools for integrating phylogenies and ecology. *Bioinformatics*, 26(11), 1463-1464. https://doi.org/10.1093/bioinformatics/btq166

Kuhn, T. S., Mooers, A., & Thomas, G. H. (2011). A simple polytomy solver for dated phylogenies. *Methods in Ecology and Evolution*, 2(5), 427-436. https://doi.org/10.1111/j.2041-210X.2011.00103.x

Lindenmayer, D. B., Likens, G. E., & Franklin, J. F. (2010). Rapid responses to facilitate ecological discoveries from major disturbances. *Frontiers in Ecology and the Environment*, 8(10), 527-532. https://doi.org/10.1890/090184

MacArthur, R. H., & MacArthur, J. W. (1961). On bird species diversity. *Ecology*, 42(3), 594-598. https://doi.org/10.2307/1932254
Mancini, D., Priore, R., Battaglia, D., & Achterberg, C. V. (2003). *Caenopachys hartigii* (Ratzeburg)(Hymenoptera: Braconidae: Doryctinae) confirmed for Italy, with notes on the status of the genus Caenopachys Foerster. *Zoologische Mededelingen Zool. Med. Leiden*, 77(4), 459-470.

Murdoch, W. W. (1977). Stabilizing effects of spatial heterogeneity in predator-prey systems. *Theoretical Population Biology*, 11(2), 252-273. https://doi.org/10.1016/0040-5809(77)90028-4

Noyes, J. S. (2019). *Universal Chalcidoidea Database*. https://ci.nii.ac.jp/naid/10026455980/

Parisi, F., Pioli, S., Lombardi, F., Fravolini, G., Marchetti, M., & Tognetti, R. (2018). Linking deadwood traits with saproxylic invertebrates and fungi in European forests - A review. *I Forest*, 11(June), 423-436. https://doi.org/10.3832/ifor2670-011

Price, P. W. (1973). Reproductive strategies in parasitoid wasps. *The American Naturalist*, 107(957), 684-693. https://doi.org/10.1086/282867

Prinzing, A., Reiffers, R., Braakhekke, W. G., Hennekens, S. M., Tackenberg, O., Ozinga, W. A., Schamine, J. H. J., & van Groenendael, J. M. (2008). Less lineages more trait variation: Phylogenetically clustered plant communities are functionally more diverse. *Ecology Letters*, 11(8), 809-819. https://doi.org/10.1111/j.1461-0248.2008.01189.x

Ratnasingham, S., & Hebert, P. (2007). The barcode of life data system BOLD: (www.barcodinglife.org). *Molecular Ecology Notes*, 7, 355-364. https://doi.org/10.1111/j.1471-8286.2006.01678.x

Ratnasingham, S., & Hebert, P. D. N. (2013). A DNA-based registry for all animal species: The barcode index number (BIN) system. *PLoS One*, 8(7), 1-16. https://doi.org/10.1371/journal.pone.0066213

Schmidl, J., & Bußler, H. (2004). Ökologische Gilden xylobionter Käfer Deutschlands. *Naturschutz Und Landschaftsplanung*, 36(7), 202-218.

Schuler, M. S., Chase, J. M., & Knight, T. M. (2015). More individuals drive the species energy-area relationship in an experimental zooplankton community. *Oikos*, 124(8), 1065-1070. https://doi.org/10.1111/oik.01931

Seibold, S., Brandl, R., Buse, J., Hothorn, T., Schmidl, J., Thorn, S., & Müller, J. (2015). Association of extinction risk of saproxylic beetles with ecological degradation of forests in Europe. *Conservation Biology*, 29(2), 382-390. https://doi.org/10.1111/cobi.12427

Strand, M. R., & Pech, L. L. (1995). Immunological basis for compatibility in parasitoid-host relationships. *Annual Review of Entomology*, 40(1), 31-56. https://doi.org/10.1146/annurev.en.40.010195.000335

Straub, C. S., Ives, A. R., & Gratton, C. (2011). Evidence for a trade-off between host-range breadth and host-use efficiency in aphid parasitoids. *American Naturalist*, 177(3), 389-395. https://doi.org/10.1086/658177

Strong, D., Lawton, J., & Southwood, R. (1984). *Insects and plants. Community patterns and mechanisms*. Harvard University Press.

Sun, X., Yang, Q., Sweeney, J. D., & Gao, C. (2006). A review: Chemical ecology of *Ips typographus* (Coleoptera, Scolytidae). *Journal of Forestry Research*, 17(1), 65-70. https://doi.org/10.1007/s11676-006-0016-2

Sverdrup-Thygeson, A., & Birkemoe, T. (2009). What window traps can tell us: Effect of placement, forest openness and beetle reproduction in retention trees. *Journal of Insect Conservation*, 13(2), 183-191. https://doi.org/10.1007/s10841-008-9141-x

Tews, J., Brose, U., Grimm, V., Tielbörger, K., Wichmann, M. C., Schwager, M., & Jeltsch, F. (2004). Animal species diversity driven by habitat heterogeneity/diversity: The importance of keystone structures. *Journal of Biogeography*, 31(1), 79-92. https://doi.org/10.1046/j.0305-0270.2003.00994.x

Thorn, S., Bässler, C., Bußler, H., Lindenmayer, D. B., Schmidt, S., Seibold, S., Wende, B., & Müller, J. (2016). Bark-scratching of storm-felled trees preserves biodiversity at lower economic costs compared to debarking. *Forest Ecology and Management*, 364, 10-16. https://doi.org/10.1016/j.foreco.2015.12.044

Thorn, S., Bässler, C., Gottschalk, T., Hothorn, T., Bußler, H., Raffa, K., & Müller, J. (2014). New insights into the consequences of post-windthrow salvage logging revealed by functional structure of saproxylic beetles assemblages. *PLoS One*, 9(7), e101757. https://doi.org/10.1371/journal.pone.0101757

Thorn, S., Bässler, C., Svoboda, M., & Müller, J. (2017). Effects of natural disturbances and salvage logging on biodiversity - Lessons from the Bohemian Forest. *Forest Ecology and Management*, 388(11), 113-119. https://doi.org/10.1016/j.foreco.2016.06.006

Thorn, S., Bußler, H., Fritze, M.-A., Goeder, P., Müller, J., Weiß, I., & Seibold, S. (2016). Canopy closure determines arthropod assemblages in microhabitats created by windstorms and salvage logging. *Forest Ecology and Management*, 381, 188-195. https://doi.org/10.1016/j.foreco.2016.09.029

Vega, F. E., & Hofstetter, R. (2015). *Bark beetles-Biology and ecology of native and invasive species*, 1st ed. Academic Press.

Vialatte, A., Bailey, R. I., Vasseur, C., Matocq, A., Gossner, M. M., Everhart, D., Vitrac, X., Belhadj, A., Ernoult, A., & Prinzing, A. (2010). Phylogenetic isolation of host trees affects assembly of local Heteroptera communities. *Proceedings of the Royal Society B: Biological Sciences*, 277, 2227-2236. https://doi.org/10.1098/rspb.2010.0365

Vinson, S. B. (1998). The general host selection behavior of parasitoid hymenoptera and a comparison of initial strategies utilized by Larvaphagous and Oophagous species. *Biological Control*, 11(2), 79-96. https://doi.org/10.1006/bcon.1997.0601

Wegensteiner, R., Wermelinger, B., & Herrmann, M. (2015). Natural enemies of bark beetles: Predators, parasitoids, pathogens, and nematodes. In F. Vega, & R. Hofstetter (Eds.), *Bark beetles - Biology and ecology of native and invasive species* (1st ed., pp. 247-304). Academic Press.

Wermelinger, B. (2004). Ecology and management of the spruce bark beetle *Ips typographus* - A review of recent research. *Forest Ecology and Management*, 202(1-3), 67-82. https://doi.org/10.1016/j.foreco.2004.07.018

Yguel, B., Bailey, R. I., Villemant, C., Brault, A., Jactel, H., & Prinzing, A. (2014). Insect herbivores should follow plants escaping their relatives. *Oecologia*, 176(2), 521-532. https://doi.org/10.1007/s00442-014-3026-3

Yu, D. S. K., van Achterberg, C., & Horstmann, K. (2016). *Taxapad 2016, Ichneumonoidea 2016, Ottawa, Ontario*. Database on flash-drive. http://www.taxapad.com

**How to cite this article:** Vogel S, Prinzing A, Bußler H, Müller J, Schmidt S, Thorn S. Abundance, not diversity, of host beetle communities determines abundance and diversity of parasitoids in deadwood. *Ecol Evol.* 2021;00:1-8.
https://doi.org/10.1002/ece3.7535