Tropical dung beetle morphological traits predict functional traits and show intraspecific differences across land uses.

Functional traits and functional diversity measures are increasingly being used to examine land use effects on biodiversity and community assembly rules. Morphological traits are often used directly as functional traits. However, behavioral characteristics are more difficult to measure. Establishing methods to derive behavioral traits from morphological measurements is necessary to facilitate their inclusion in functional diversity analyses. We collected morphometric data from over 1,700 individuals of 12 species of dung beetle to establish whether morphological measurements can be used as predictors of behavioral traits. We also compared morphology among individuals collected from different land uses (primary forest, logged forest, and oil palm plantation) to identify whether intraspecific differences in morphology vary among land use types. We show that leg and eye measurements can be used to predict dung beetle nesting behavior and period of activity and we used this information to confirm the previously unresolved nesting behavior for Synapsis ritsemae. We found intraspecific differences in morphological traits across different land use types. Phenotypic plasticity was found for traits associated with dispersal (wing aspect ratio and wing loading) and reproductive capacity (abdomen size). The ability to predict behavioral functional traits from morphology is useful where the behavior of individuals cannot be directly observed, especially in tropical environments where the ecology of many species is poorly understood. In addition, we provide evidence that land use change can cause phenotypic plasticity in tropical dung beetle species. Our results reinforce recent calls for intraspecific variation in traits to receive more attention within community ecology.

INTRODUCTION

Trait‐based approaches are considered essential for establishing general principles in community ecology (McGill, Enquist, Weiher, & Westoby, 2006). Functional traits describe morpho‐physio‐phenological traits which indirectly impact an organism's fitness through their effects on growth, reproduction, and survival (Violle et al., 2007). The distribution and frequency of functional traits can be used to calculate a range of indices for a sample, for example functional diversity, evenness, and dispersion (Albert et al., 2012; Mouillot, Graham, Villéger, Mason, & Bellwood, 2013). Both functional traits and associated diversity measures are increasingly being used to examine human impacts on biodiversity, and the consequences for key ecosystem functions and services (Mouillot et al., 2013; Flynn & Gogol‐Prokurat, 2009; Hamer et al. 2015). The study of functional traits is also providing insights into how communities are structured (Messier, McGill, & Lechowicz, 2010; Violle et al., 2012) and the vulnerability of communities to future disturbance (Laliberté et al., 2010).

Morphological measurements (e.g., body mass, beak size, egg size) are frequently used directly as functional traits, as morphology typically reflects an organism's interaction with the environment and the functions it performs (Bertossa, 2011; Wainwright, 1996). These morphological traits are measured on individuals caught in the field (e.g., Edwards et al., 2014; Yates & Andrew, 2011), housed in museum collections (Bregman, Sekercioglu, & Tobias, 2014; Matthews, Sheard, & Cottee‐Jones, 2015), or taken from existing literature (e.g., Edwards, Edwards, Hamer, & Davies, 2013; Flynn & Gogol‐Prokurat, 2009). In contrast, behavioral characteristics that relate to an animal's function generally need to be observed directly from living individuals and can be more difficult to obtain. Establishing methods to derive behavioral traits from morphological measurements would therefore facilitate their inclusion in functional diversity analyses (Bertossa, 2011). In turn, the inclusion of behavioral traits would then increase the utility of functional trait analyses for ecological research.

The use of mean trait values in calculations of functional diversity indices is necessary where the resources are not available to measure all individuals, but is vulnerable to biases where sample sizes are small and could overlook important within‐species variation across different treatments or habitat types (Albert, Thuiller, Yoccoz, Douzet et al., 2010; Griffiths, Louzada, Bardgett, & Barlow, 2016; Hulshof & Swenson, 2010). The latter is a particular risk for species that show high morphological plasticity, as morphological traits could vary greatly among individuals of the same species in different habitats (Albert, Thuiller, Yoccoz, Soudant et al., 2010; Petchey & Gaston, 2006). The extent and magnitude of intraspecific variation in traits needs to be understood in order to interpret the effects of land use change, the processes structuring communities, such as species coexistence (Violle et al., 2012), and the interactions between populations (Bolnick et al., 2011). Despite an increasing interest in intraspecific trait variation in ecological research (Albert et al., 2012; Griffiths et al., 2016; Violle et al., 2012), how traits vary across different land uses has not been widely examined outside of producer systems, and variability in animal traits has been largely overlooked (Albert, Thuiller, Yoccoz, Soudant et al., 2010; Albert, Grassein, Schurr, Vieilledent, & Violle, 2011 but see Moretti et al., 2017).

Dung beetles (Coleoptera: Scarabaeidae) are an ideal group in which to study variation in functional traits within and among species, as well as across different land uses. Dung beetles can be separated into distinct functional groups based on body size, period of activity (diurnal/nocturnal), and nesting behavior. They can nest within a dung pat (endocoprids or dwellers), beneath the dung (paracoprids or tunnellers), or by rolling a dung ball horizontally away from the source (telocoprids or rollers) (Hanski & Cambefort, 1991). Dung beetle species can also vary in temperature tolerance (Addo‐Bediako, Chown, & Gaston, 2000) and food resource use (Filgueiras, Liberal, Aguiar, Hernández, & Iannuzzi, 2009; Hanski & Cambefort, 1991; Holter, Scholtz, & Wardhaugh, 2002). All of these functional traits directly affect the ecosystem functions performed by dung beetles, including dung removal, bioturbation, and secondary seed dispersal (Griffiths et al., 2015; Nichols et al., 2008; Slade, Mann, Villanueva, & Lewis, 2007).

To date, dung beetle behavioral traits have been derived from direct observations of dung beetle activity. However, the natural history of many species remains incomplete particularly in the tropics (Nichols, Uriarte, Bunker, & Favila, 2013). Phylogenetic relatedness is often used to predict the behavioral ecology of a species or genus where information is missing (Barton, Gibb, Manning, Lindenmayer, & Cunningham, 2011; Bregman et al., 2014; Hunt, Bergsten, & Levkanicova, 2007; Matthews et al., 2015). However, in dung beetles, some behaviors (e.g., rolling behavior) may have evolved many times from tunneling clades (Monaghan, Inward, Hunt, & Vogler, 2007), indicating that generalizations based on phylogeny alone may be an inaccurate indication of dung beetle behavioral ecology. Dung beetle morphology is known to vary between nesting behaviors (Gregory, Gómez, Oliveira, & Nichols, 2015; Hanski & Cambefort, 1991; Inward, Davies, Pergande, Denham, & Vogler, 2011; Shafiei, Moczek, & Nijhout, 2001; Simmons & Ridsdill‐Smith, 2011) and activity period (Baird, Byrne, & Scholtz, 2010; Caveney, Scholtz, & McIntyre, 1995; Hernández, Monteiro, & Favila, 2011; McIntyre & Caveney, 1998). Therefore, morphological data could be used to predict dung beetle behavioral traits more reliably than using the behavioral traits of congeneric species.

Dung beetles have been shown to display high phenotypic plasticity in certain traits in response to resource availability (Moczek & Nijhout, 2004) and environmental change (Alves & Hernández, 2017; Emlen, 1994; Moczek, 2002; Pomfret & Knell, 2006). In particular, the morphology of males of the same species has been found to vary greatly as a result of resource quality (Moczek & Emlen, 2000; Moczek & Nijhout, 2004), and differences in horn size and shape, both within and between species, are well documented (e.g., Emlen, Marangelo, Ball, & Cunningham, 2005). Intraspecific differences in morphological traits are therefore likely to be easily detected where they occur, making dung beetles the ideal study group to assess whether land use change and associated changes in resource availability may be causing intraspecific differences in functional traits.

Here, we test whether morphological traits can be used to reliably identify the behavioral traits of dung beetle species, and whether there is intraspecific variation in these traits across different land uses. Our focus is on tropical forested landscapes, where logging and conversion to agriculture continue to cause a reduction in species richness and the abundance of many taxa, including dung beetles (Gray, Slade, Mann, & Lewis, 2014; Nichols et al., 2007; Slade, Mann, & Lewis, 2011). These changes correspond to a loss of functional diversity in dung beetles (Edwards et al., 2013, 2014; Gray et al., 2014; Senior et al., 2013) and impact ecosystem functions (Slade et al., 2011). We used measurements from dung beetle species collected in old‐growth forest, logged forest, and oil palm plantations in Sabah, Borneo, to ask the following questions:

Can morphological characteristics be used to separate known behavioral functional traits in dung beetles?

Is there intraspecific variation in dung beetle morphological traits across different tropical land uses?

METHODS

Data collection

Dung beetles (Coleoptera: Scarabaeinae) were sampled in Sabah, Malaysian Borneo during February and March 2011. We sampled nine sampling points at three sites in each of three land uses (n = 27): old‐growth forest (OG) which had never been logged; logged forest (LF) which had experienced 2–3 periods of timber extraction, and oil palm plantations (OP) in which palms were planted between 2000 and 2006. All sampling points were located within the Stability of Altered Forest Ecosystems (S.A.F.E.) project, encompassing the old‐growth forest at Maliau Basin, (see Reynolds, Payne, Sinun, Mosigil, & Walsh, 2011; Ewers et al., 2011 for further details). The nine sampling points in each site were 178 m apart, and the three sites in each land use were 1.5 km apart. At each sampling point, we placed a pitfall trap baited with 25 g of human feces and filled with a solution of water, salt, and detergent (Gray et al., 2014). Trap contents were collected after 48 hours and stored in 75% ethanol in a freezer until beetles were identified and measurements taken.

Specimens were identified using Balthasar (1963), Boucomont (1914), the works on Bornean Scarabaeinae by Ochi and Kon (e.g., Ochi, 1996), and the reference collections housed in the Oxford University Museum of Natural History (OUMNH). Voucher collections are deposited at the OUMNH, Universiti Malaysia Sabah, and the Forest Research Centre in Sandakan, Sabah.

To ensure a minimum of 25 individuals from each species from each site were measured, we focused on 12 of the most abundant and widespread species, including at least one from each of the major functional groups. We combined the samples from all nine traps from each site. We measured as many randomly selected individuals as possible in the available time (n = 25–36 individuals per species per site). Species names, functional groups (as identified by Slade et al., 2007), and the number of individuals measured per species are given in Supporting information Table S1.

We were specifically interested in the species Synapsis ritsemae (Lansberge) as its nesting behavior is unresolved. Synapsis ritsemae is the only member of its genus to occur in Borneo. This genus has traditionally been classified as a nocturnal tunneller based on its taxonomic position within the Coprini (Davis, Scholtz, & Philips, 2002; Philips, Pretorius, & Scholtz, 2004), and field observations of Masumoto (1973), yet species within this genus have been recorded to perform dung ball formation for rolling as well as tunneling activity (Kon, Ochi, Kusakabe, & Araya, 2004; Zidek & Pokorny, 2010). If morphological traits alone can be reliably used to separate species according to their behaviors, we may be able to resolve the behavior of Synapsis ritsemae from morphological measurements. As no individuals were collected from pitfall trap contents in this study, 29 mounted specimens from earlier work at nearby locations in Sabah, Borneo (Slade et al., 2011), were measured to address this question. As samples were already dry, wing measurements could not be taken from the mounted specimens.

For each individual beetle, we measured 13 linear dimensions and one area (details given in Supporting information Table S2 and Figure S1) using a LEICA M165 stereo microscope and morphometric software Leica Application Suite (version 3.0), from which we calculated a set of six morphological traits of potential importance for ecological functioning: relative body size, abdomen size, wing aspect ratio, wing loading, eye size, and hind leg size (description of each measurement given in Table 1). Measurements were always taken on the left appendages to account for asymmetry between sides. We took measurements of tibial width across all species at approximately one‐quarter of the leg length before the base to avoid the issue of tibial wear from digging which can vary based on beetle age. For wing measurements, wing cases were opened with forceps and the left wing removed. The wing was flattened on a glass slide using a brush and water to remove creases. All morphological measurements were accurate to ±0.005 mm. The two main forms of male dung beetles (major and minor) frequently differ in size (Simmons & Ridsdill‐Smith, 2011). We therefore recorded the sex and phenotype of the beetles (female, major male, minor male).

Table 1

Behavioral traits and associated morphological traits, with explanation of the morphometric measurements taken

Behavioral trait	Morphological trait	Morphometric measurement	Values near unity indicate:
Activity period	Eye size (Caveney et al., 1995)	Eye Length/Body Length	A large eye relative to body length
Reproductive capacity	Abdomen size (Srygley & Chai, 1990)	Thorax Length/Abdomen Length	A large abdomen relative to thorax
Dispersal	Wing aspect ratio (Berwaerts, Van Dyck, & Aerts, 2002)	Wing Width/Wing Length	A thin, long wing
	Wing loading (Berwaerts et al., 2002)	Wing Area/Body Area (Body Area = Body Length × Thorax Width)	A large wing area relative to body area
Relative body size	Body size (Marden, 2000)	Body Length/Body Length of that species	A large body size
Nesting behavior	Hind leg size	Hind Leg Length/Body Length	A long leg relative to body length
	Hind leg robustness	Hind Leg width/Hind Leg length	A wide, short leg

ANALYSIS

Can morphological characteristics be used to determine behavioral functional traits?

We ran a permutational analysis of variance (PERMANOVA) to test whether morphology differed between functional groups. Data on hind leg size, hind leg robustness, eye roundness, and eye size were used as the response variables, with functional group (diurnal tunneller, nocturnal tunneller, diurnal roller, nocturnal roller, and nocturnal roller/tunneller) as the predictor. Site was used as a grouping variable to account for spatial autocorrelation. To visualize the differences between functional groups we used a principal component analysis (PCA).

Is there intraspecific variation in morphological traits across different tropical land uses?

We tested whether dung beetle morphological traits varied across land uses with general linear mixed models (GLMMs). We specified land use (old‐growth, logged, oil palm), species, sex (female, major male), and the two‐way interactions Land use:Species and Land use:Sex as fixed effects, and set site as a random factor. Model selection was based on p values, starting with the full model and dropping terms for which model comparison using chi‐squared tests gave p > 0.05. A Gaussian error distribution was appropriate for all response variables. We excluded data on minor males, as there were insufficient measurements for all species across all land uses to estimate the Land use:Sex interaction with sex as a three‐level factor (female, major male, minor male). We selected focal species that occurred in more than one land use type. However, because dung beetles show strong affinities to particular land uses, many species were not found in all three land uses. Therefore, for each morphological trait we ran separate models for the species occurring in:

Set A) Old‐growth and logged forest (n = 1,315 individuals measured from 11 species across six sites: Catharsius dayacus, Copris sinicus, Microcopris doriae, Onthophagus obscurior, Onthophagus rugicollis, Onthophagus vulpes, Paragymnopleurus maurus, Paragymnopleurus sparsus, Paragymnopleurus striatus, Proagaderus watanabei, and Sisyphus thoracicus).

Set B) Logged forest and oil palm (n = 408 individuals measured from four species across six sites: Catharsius renaudpauliana, Onthophagus obscurior, Onthophagus rugicollis, and Proagaderus watanabei).

As we were unable to incorporate data on minor males into the models described above, we ran a separate analysis on the ratio of major to minor males to assess whether there were differences in the relative abundances of the two phenotypes across land use types. We used only data from tunnellers, as roller species do not have male morphs that can be clearly distinguished based on morphological traits. For each site, we calculated the number of major and minor males in the random sample and used this as a two‐column response variable in a binomial generalized linear model (GLM) with a quasibinomial error structure to account for overdispersion. Land use, species, and the two‐way interaction were fixed effects. As above, we separated data into comparisons of logged forest and old‐growth forest (39 comparisons of seven species across six sites) and logged forest versus oil palm (21 comparisons of four species across six sites), and used model selection as described above.

All analyses were run using the R software (version 3.0.2, R Core Team, 2017) and the packages vegan (version 2.0‐10, (Oksanen et al., 2017) and lme4 (version 1.1‐7, Bates, Maechler, Bolker, & Walker, 2013).

RESULTS

We found statistically significant differences among the morphological traits (hind leg size, hind leg robustness, eye size, and eye roundness) of different functional groups (PERMANOVA, F 4,1715 = 1182.7, p < 0.001). Differences between functional groups were visualized using a principal components analysis (Figure 1); PC1 explained 56.3% of the variation in the data and PC2 explained 38.6% of the total variation. More positive values of PC1 and more negative values of PC2 represent larger leg size and lower leg robustness. More positive values of both PC1 and PC2 represent rounder, smaller eyes. The tunneller and roller species separated into clusters with the rollers having greater hind leg size (relative to body size) and the tunnellers greater leg robustness (i.e., greater width for any given leg length). The nocturnal and diurnal species also separated out, with the diurnal species having rounder (i.e., greater width for any given length), smaller eyes than the nocturnal species, although this distinction was more pronounced for the tunnellers. Nocturnal rollers were represented by only one species, Paragymnopleurus striatus. Individuals of this species had larger, less round eyes than many of the diurnal rollers, but still fell within the overall range of values for diurnal rollers. Synapsis ritsemae fell clearly between the tunneller and roller clusters, but was slightly closer to the tunneller cluster. Synapsis ritsemae also clustered nearer the nocturnal functional groups.

Figure 1

Principal component analysis of data on hind leg size, hind leg robustness, eye size, and eye roundness (see Table 2). Plot shows values for individual beetles and standard deviation of the data around the mean. The tunneller and facultative roller Synapsis ritsemae (■) falls between the tunneller and roller functional groups

Table 2

Results of model selection based on Chi‐squared tests and p‐values. A separate model selection process was run for each morphological trait in each of the data sets. χ2, df and p‐values are given for comparisons between models with and without the term specified – following term removal. NB – in set A, C. sinicus has only 5 major males in old growth, and P. striatus only 5 individuals in total in logged forest (4 female, 1 male)

		Set A: Old growth and logged forest (11 species)			Set B: Logged forest and oil palm (4 species)
Response	Terms	Df	χ2	p	Df	χ2	p
Relative body size	Land Use	1	0.9	0.344	1	0.02	0.898
	Sex	1	65.2	<0.001	1	35.3	<0.001
	Species	10	4450	<0.001	3	1520	<0.001
	Land Use:Species	10	29.5	<0.001	3	0.31	0.957
	Land Use:Sex	1	0.01	0.907	1	0.57	0.451
Abdomen size	Land Use	1	9.08	0.003	1	2.42	0.12
	Sex	1	120	<0.001	1	38.9	<0.001
	Species	10	2440	<0.001	3	516	<0.001
	Land Use:Species	10	98.3	<0.001	3	10.3	0.017
	Land Use:Sex	1	0.04	0.838	1	0.57	0.452
Wing aspect ratio	Land Use	1	12.1	<0.001	1	3.36	0.067
	Sex	1	10.5	0.001	1	4.62	0.032
	Species	10	2220	<0.001	3	378	<0.001
	Land Use:Species	10	15.4	0.118	3	1.32	0.724
	Land Use:Sex	1	0.58	0.445	1	0.44	0.505
Wing loading	Land Use	1	10.3	0.001	1	0.82	0.365
	Sex	1	4.79	0.029	1	7.53	0.006
	Species	10	3720	<0.001	3	849	<0.001
	Land Use:Species	10	29.9	0.001	3	7.55	0.056
	Land Use:Sex	1	0	0.987	1	0.57	0.451
Eye size	Land Use	1	0.26	0.612	1	0.13	0.719
	Sex	1	0.2	0.652	1	5.43	0.02
	Species	10	3170	<0.001	3	1130	<0.001
	Land Use:Species	10	35.1	<0.001	3	1.84	0.606
	Land Use:Sex	1	2.12	0.146	1	0.18	0.672
Hind leg size	Land Use	1	7.26	0.007	1	4.59	0.032
	Sex	1	5.78	0.016	1	1.83	0.177
	Species	10	4670	<0.001	3	325	<0.001
	Land Use:Species	10	31	0.001	3	12.1	0.007
	Land Use:Sex	1	2.26	0.132	1	0.29	0.59

There were significant intraspecific morphological differences between dung beetles found in old‐growth and logged forest for all but one of the traits we analyzed (Table 2, Figure 2). The Land use:Species term was significant for all morphological traits apart from wing aspect ratio, suggesting that there are intraspecific differences in body size, abdomen size, wing loading, hind leg size, and eye size among land uses. For example, four species (Catharsius dayacus and three Onthophagus species) showed significantly higher wing loading in old‐growth relative to logged forest sites, and one species (Proagaderus watanabei) showed lower wing loading (Figure 2). The modeled mean increase in wing loading corresponds to a 6% increase in the ratio of wing area to abdomen area overall. Similarly, diurnal tunnellers that showed a change in abdomen size were found to be larger in old‐growth than that of individuals from logged forest, although the magnitude of the difference varied between species. For example, the greatest difference in abdomen size (the ratio of abdomen to thorax: higher values indicate larger abdomen, see Table 1) was shown by Onthophagus obscurior; the mean abdomen size of individuals was 16% larger in old‐growth compared to logged forest. Intraspecific differences in body size between old‐growth and logged forest were driven by an increase in body size of Catharsius dayacus in logged forest.

Figure 2

Standardized regression coefficients for morphological trait response variables for individuals in old‐growth forest relative to individuals from the same species in logged forest for each species in Set A where a significant Land use:Species interaction was observed. Responses less than zero indicate a decline the response variable moving from logged forest to old‐growth; values above zero indicate an increase. Figures show means and 95% CI obtained from model coefficients, note variation in y‐axes. Each response variable is a ratio (see Table 1) [Colour figure can be viewed at http://wileyonlinelibrary.com]

For wing aspect ratio (the only response variable for which no significant Land use:Species interaction was found when comparing old‐growth and logged forest), Land Use alone was significant (Table 2). Wing aspect ratio was found to be consistently smaller for individuals of all species found in logged forest in comparison with individuals found in old‐growth. The model coefficients corresponded to a 1% decrease in the ratio of wing width to wing length from old‐growth to logged forest. The Land use:Sex interaction was not significant for any traits compared between old‐growth and logged forest, indicating that both females and major males showed similar differences across land use types. For the species found in old growth and logged forest, all morphological measurements apart from eye size differed between sexes. Males showed larger values for all measurements apart from wing loading, which was larger in females.

Intraspecific trait differences between oil palm and logged forest only occurred in two of the six morphological traits measured: hind leg size and abdomen size (Table 2, Figure 3). For both of these traits, these results were driven by one species, the diurnal tunneller, Onthophagus rugicollis, that had larger hind leg size and abdomen size in oil palm compared to logged forest. The other species showed no response in relation to land use (oil palm versus logged forest). There were no significant Land use:Sex interactions; females and major males showed similar differences between logged forest and oil palm. Across all species found in logged forest and oil palm, morphological measurements differed significantly between males and females apart from hind leg size. In each case, males were larger apart from wing loading, which was larger in females.

Figure 3

Standardized regression coefficients for morphological trait response variables for individuals in oil palm relative to individuals from the same species in logged forest for each species in Set B (logged forest:oil palm) where a significant Land use:Species interaction was observed. Responses less than zero indicate a decline the response variable moving from logged forest to oil palm; values above zero indicate an increase. Figures show means and 95% CI obtained from model coefficients, note variation in y‐axes. Each response variable is a ratio (see Table 1) [Colour figure can be viewed at http://wileyonlinelibrary.com]

There were three species that were found in all three land uses (old‐growth forest, logged forest, and oil palm: Proagaderus watanabei, Onthophagus obscurior, and Onthophagus rugicollis). By comparing Figures 2 and 3, it is evident that there were no consistent patterns of traits increasing or decreasing in magnitude across the land use gradient.

There were no statistically significant differences between the proportion of major and minor males in species found in logged forest compared to old‐growth forest (F 1,37 = 1.39, p = 0.248), nor in species found in logged forest compared to oil palm (F 1,19 = 0.16, p = 0.693).

DISCUSSION

Given the increasing importance of functional trait measures in community ecology, we sought to address two issues to help inform the use of morphological trait measurements. First, whether behavioral traits can be predicted from morphological measurements and second, whether there are intraspecific differences in morphological traits across land use types. Our analyses demonstrate that the functional groups of activity period and nesting behavior in dung beetles can be resolved from specific morphological measurements. We also found that for several morphological traits there are statistically significant differences among individuals of the same species from different land uses. Phenotypic plasticity was found for traits associated with dispersal (wing aspect ratio) and reproductive capacity (abdomen size).

Dung beetle functional groups are most often classified based on direct observation of behavior (e.g., Doube, 1990). This can be time‐consuming and a lack of data can lead to inaccuracies. It also may not be possible to observe all behavioral traits. Our results show that once behaviors and morphological measurements are known for some species, is it possible to resolve both nesting behavior and activity period (two of the most common traits used in ecological studies) for new species using morphological measurements.

Rollers and tunnellers separated along PC1, which corresponded largely to hind leg morphology. The low variance within tunneller and roller clusters suggests hind leg morphology is strongly conserved within functional groups. We were particularly interested in the position of Synapsis ritsemae along these axes, as this species has been classified as a tunneller (Masumoto, 1973), despite being seen to perform ball production and rolling behavior (Kon, Ochi, Kusakabe, & Araya, 2004; Zidek & Pokorny, 2010). We found that individuals of this species were clustered between the two nesting behavior guilds, albeit slightly closer to the tunneller guilds, strongly suggesting that this species is indeed capable of both behaviors. This suggests that either Synapsis does not belong within the Coprini or is a basal or atypical coprine in terms of its phylogeny (Tarasov & Dimitrov, 2016). This highlights how relying on phylogeny to classify functional groups may be misleading, and morphological measurements may be better predictors. Interestingly, this flexible behavior, may have evolved as a response to competition; individuals were observed to roll dung balls only when competition for dung was high, suggesting that they are facultative rollers (Slade & Mann, unpublished obs.).

Nocturnal and diurnal tunneling species separated out clearly along the PC2 axis, corresponding to eye size and roundness metrics. This confirms that for tunnellers, eye morphological characteristics can be used to determine activity period. Since P. striatus has been classified as a nocturnal roller, we expected that this species would form a cluster distinct from the diurnal rollers, but this was not the case. Currently, there is an incomplete understanding of the natural history of Bornean dung beetles, and species that could potentially be crepuscular such as P. striatus may have been misclassified as nocturnal (Slade et al., 2007). Overall, this result indicates the added value of using morphological measurements to examine and confirm behavioral traits, particularly if such behaviors are difficult to observe or are unrecorded.

The approach of classifying function through morphology could be extended to distinguish more precise, refined behavioral functional group classifications for dung beetle species (such as Synapsis ritsemae in this study), or to move beyond discrete classifications of behavior to metrics on a continuous scale (Fountain‐Jones 2015; Petchey & Gaston, 2006). Based on the high intraspecific variability detected in some morphological traits and the overlap between morphology and functional groups, we suggest that morphometric data be used as an additional tool to behavioral classifications of functional traits in dung beetles. Due to the lack of a resolved phylogeny for dung beetles (Monaghan et al., 2007; Tarasov & Dimitrov, 2016), it was not possible to control for the influence that phylogenetic relatedness could have on the morphology of the species studied. Despite this potential bias, the separation between functional groups is strong enough to suggest these patterns would remain if relatedness was accounted for.

Changes in dung beetle functional diversity across different land uses have been previously recorded (Andresen, 2003; Barragán, Moreno, & Escobar, 2011; Edwards et al., 2014). However, these studies have used mean trait values which do not take into account variation within species. We found significant differences in trait measurements among land use types, both within species and across all species. Griffiths et al. (2016) have recently shown intraspecific differences in neotropical dung beetle species from distinct communities in lowland tropical rainforest. Here we show that intraspecific differences are also seen with increasing human modification. It therefore seems likely that for taxa demonstrating high morphological plasticity, such as dung beetles, functional diversity measures may be incorrect if mean species trait values are used. More studies on other taxonomic groups are needed to establish the extent to which intraspecific differences in morphology occur across different land use types and the effect that including intraspecific variation has on functional diversity indices.

We found evidence of intraspecific differences between old‐growth and logged forests for several of the traits measured. Wing aspect ratio was the only trait that did not show intraspecific variation; where relatively shorter, wider wings were found for individuals of all species in logged forest in comparison with old‐growth. Abdomen size and wing loading were found to be on average larger across all species in old‐growth in comparison with logged forest. Relatively wider and shorter but smaller wings and smaller abdomens in logged forest could be indicative of a trade‐off between dispersal traits (Gibb et al., 2006) and reproductive traits (Thomas, Hill, & Lewis, 1998). This suggests traits associated to flight ability are of particular importance for dung beetle species in disturbed habitats and changes in investment between flight ability and reproductive potential could favor survival in modified land use types.

Larger bodied dung beetle species are thought to be more vulnerable to the effects of habitat disturbance, showing both reduced abundance and species richness in degraded forest (Larsen, Williams, & Kremen, 2005; Senior et al., 2013; Slade et al., 2011). However, trends seen in body size at the interspecific level are not found at the intraspecific level, as body size did not differ among land uses for most species. There was also no difference found in the relative abundance of male major/minor phenotypes across land uses. Both body size and male phenotype vary with parental investment (Moczek, 1998), which indicates that any changes in resource availability between land use types in this landscape are small, or are not being translated into reduced brood ball quality.

Hind leg size and eye size varied between species, and land use comparisons. This suggests that changes in resource availability, temperature exposure, and habitat structure do not have the same effects on all species or all traits. For example, eye size was found to be larger for some species in logged forest in comparison with old‐growth forest, and smaller for others. Thus, any possible explanation for changes in trait values affecting visibility and eyesight, such as changes in light levels through canopy structure, cannot be acting on all species in the same way. Similarly, other processes, such as competition between species, and the trade‐offs of size in any morphological dimension, appear to be acting differently on different species.

IMPLICATIONS

We have demonstrated that behavioral traits can be resolved from morphology in dung beetles, using seven morphological traits to successfully confirm or identify five behavioral traits, associated with nesting behavior, dispersal, and activity period in 12 dung beetle species. This indicates the potential to further the use of morphological characteristics in defining functional traits in dung beetles and other groups, and their use to support and supplement the classification of behavioral traits. Our findings also support recent calls for increased reporting and use of intraspecific variation in studies of morphological and functional traits (Bolnick et al., 2011; Laughlin, Joshi, Bodegom, & Bastow, 2012; Violle et al., 2012), particularly in terrestrial animal groups (see Griffiths et al., 2016). Current studies linking dung beetle communities and ecosystem functions tend to focus on community‐level attributes such as species richness, abundance, and community composition (e.g., Barnes, Emberson, Chapman, Krell, & Didham, 2014; Filgueiras, Iannuzzi, & Leal, 2011). Given that the morphological and functional traits of dung beetles affect their contribution to ecosystem function (Slade et al., 2007), biologically meaningful intraspecific differences in functionally important traits could be driving some of the changes in ecosystem function observed with habitat conversion. These subtle intraspecific effects are not yet well understood and should be incorporated into studies of how forest degradation alters ecosystem functioning.

CONFLICT OF INTEREST

None declared.

AUTHOR CONTRIBUTIONS

CG, DJM, and EMS conceived the study. EMS collected the samples , EHR carried out morphometric data collection and DJM assisted in species identification. CG, EHR, and EMS carried out data analysis. EHR lead the writing of the manuscript with input from all authors.

DATA ACCESSIBILITY

Dung beetle morphometric data can be accessed the Dryad Digital Repository: https://doi.org/10.5061/dryad.1484dn3.

Click here for additional data file.

Click here for additional data file.