Abstract
Background:
Nitrogen (NOx, NHy) and acidifying (NOx, NHy, SOx) deposition has reduced the biodiversity of European dry heathlands. Restoration efforts such as sod-cutting (removal of vegetation, litter and humus layer) often shifted these systems from N to P limitation and have had limited success in restoring the invertebrate community. Possible reasons for this include the unresolved acidification and a change in food plant stoichiometry. Here, we investigate how liming and P addition change food nutritional quality and their consequences for invertebrate performance.
Methods:
We performed feeding experiments with field crickets (Gryllus campestris), using plant material collected from a full factorial field experiment with liming and P addition. We related female reproduction as measure of individual fitness to elemental ratios of plants fed to the crickets.
Results:
P addition stimulated cricket daily reproduction and shortened their reproductive period, resulting in no difference in total reproduction. Liming greatly reduced both daily and total reproduction and resulted in more females cannibalizing on their male mates. Females that did so could partly offset the liming induced reduction in reproduction, suggesting dietary deficiency. P-addition improved food quality (lower N:P ratios) while liming led to skewed Mn:Mg and Fe:Mg ratios that compare unfavorably to ratios found in terrestrial invertebrates.
Conclusion:
Increased plant N:P ratio following sod-cutting constrains the reproductive potential in Gryllus campestris in a non-linear way. Liming reduced nutritional quality, likely by inducing deficiencies in Fe or Mn.
Management Implications:
High-impact restoration management practices such as sod cutting and liming cause new problems for invertebrates rooted in ecological stoichiometry. Since P-addition only partially offsets these negative effects, we instead advocate the use of less intensive N removal management and weaker buffering agents to reduce soil acidification. Furthermore, a reduction in N emission is paramount as it will remove the need for disruptive interventions.
Introduction
Increased deposition of nitrogen (NOx, NHy) alter soil biochemical cycles (), result in increased acidification rates (; ) and impact plant and invertebrate diversity (; ; ). Western European dry heathlands are naturally oligotrophic and poorly buffered systems and therefore highly vulnerable to increased N (NOx, NHy) and acidifying (NOx, NHy, SOx) deposition (). In North Western Europe, current annual N deposition rates often exceed critical loads, and intensified management is needed to prevent N accumulation (). In dry heathland areas subjected to high N deposition, restoration management focuses on maximizing removal of biomass and consequently N, using high-impact measures such as sod cutting (removal of vegetation, litter and humus layer) across large scales (; ; ). Such high-impact, large scaled management practices have been criticized however, as they can harm invertebrate populations due to disturbance or habitat destruction and homogenization (). Also, current management aimed at maximizing N removal often indiscriminately removes other nutrients as well, including P and base cations. This has led to undesirable shifts in soil nutrient stoichiometry, shifting these systems from N to P limitation (; ; ). Increased availability of inorganic N greatly increases aboveground plant tissue N:P ratios () and decreased plant available P has been found to further increase plant tissue N:P ratios in dry heathlands (). Such shifts in plant tissue N:P ratio has been implicated in the decline of dry heathland invertebrate diversity (). Compared to plants, invertebrates have in general a lower tissue N:P ratio and also a lower plasticity in this and other macronutrient ratios (). Consequently, shifts in plant tissue N:P ratio can greatly impact invertebrate population composition and dynamics (). In European dry heathlands that underwent sod-cutting management, there is a need to address how stoichiometric mismatches created by sod-cutting can be mitigated, and to test whether such mitigation is effective in restoring invertebrate diversity.
In addition to creating stoichiometric imbalances, management aimed at N removal does not restore soil buffer capacity. Increased soil acidity is known to hamper the restoration of characteristic heathland flora and fauna as well (; ; ; ; ), and additional measures to increase soil buffer capacity are often necessary in order to achieve restoration goals (, ; ). Increased acidification results in changed elemental content of plant tissue, reducing cations (Ca, Mg, Na, and K), and increasing Al and heavy metal (Cd and Pb) content of plants (). Depending on soil type and species, concentrations and ratios of trace elementals (Mn, Cu, and Zn) can also shift considerably. Such shifts in elemental content or ratios can ultimately result in a significant reduction in uptake of essential trace elements by invertebrate consumers. For example, reduced Ca intake resulting from soil acidification has been linked to reduced hatchling success in insectivorous birds (; ; ). Whether such changes also negatively impact fitness in invertebrates is however largely unexplored (). Variation in plant (trace) elemental content has however been found to influence invertebrate herbivore density and their species composition (), the degree of herbivory () and the fitness of honeybees (), suggesting that the availability of these elements play a role in the ecological functioning of invertebrate communities. Deficiencies that arise may however likely relate to trace elements other than Ca, as demands are much lower in invertebrates compared to birds and other vertebrates ().
Here we test the effectiveness of additional management following sod-cutting from an invertebrate’s perspective. Specifically we examine the effects of liming (to ameliorate acidification), P addition (to resolve P limitation created by sod-cutting) and their interaction, using the field cricket (Gryllus campestris) as model species. This is a large insect and mostly confined to xero-thermophilic habitats in its northern distribution range. Dry heathland is the typical habitat of this species in the area of study (). It overwinters as half grown nymph, adults emerge in early to mid-spring. Female crickets are able to produce a large amount of eggs, which develop without diapause and hatch approximately 2 weeks after oviposition (). Earlier feeding experiments with plant tissue differing in N:P ratio indicated that this species can exhibit P limited growth (). We expected that P replenishment as well as liming would improve the nutritional quality of plant tissue quality for female Gryllus campestris, by reducing plant N:P ratio to more optimal conditions (P-replenishment), and by increasing Ca, Mg and trace elemental concentration in food plants via reduced root Al stress, increasing root elemental uptake efficiency. We tested this by performing a climate controlled feeding experiment using field collected plant material from a full factorial management experiment in recently sod cut dry heathland vegetation (see ; for soil mesofauna and vegetation effects) as food source and used female reproductive success as measure of fitness response.
Materials and Methods
Field Experiment
Experimental plots were established in Hoge Veluwe National Park, in the center of the Netherlands (52°.2’0.30”N; 5°.49’0.50”) at a site that underwent sod cutting management in the first quarter of 2012. Before sod cutting the vegetation was highly dominated by Molinia caerulea. At the start of the experiment in autumn 2012, the soil was still bare, no vegetation had yet established. For more detailed information on the experimental site and set-up see and .
Experimental plots measured 15 m × 15 m in size and were placed in a randomized block design each consisting of four plots. In total five blocks were installed (n = 5). On the plots, the reduction of soil acidity treatment (further referred as Ca+, but note that Mg was also added in this treatment) was achieved by addition of dolomite (200 g.m–2 granular Dolokal; 80% CaCO3 and 20% MgCO3). P-replenishment treatment (further referred as P+) was achieved by adding 5 g P⋅m–2 by spraying an aqueous solution of NaH2PO4⋅2H2O, adding in total 5.67 kg NaH2PO4⋅2H2O per plot. P+ treatment (early April 2013) was carried out substantially later than Ca+ treatment (December 2012), in order to avoid high PO43– sorption resulting from liming acid soils, as reported by . Each plot received either P+ or no addition (P–), Ca+ or no addition (Ca–), resulting in four treatment combinations P+Ca+, P+ Ca–, P–Ca+ and P–Ca–. Vegetation was then allowed to develop naturally for two growing seasons in order to provide sufficient amount of food for the feeding experiment.
Feeding Experiment
Food Acquirement and Preparation
Vegetation for food preparation was collected at the experimental sites in April and May 2014. In all plots, the dominant graminoid (Molinia caerulea) and ericoid (Calluna vulgaris) was collected separately by clipping aboveground living biomass. The obtained plant tissue was pooled by treatment combination, using an equal volume of plant material per individual plot. In addition to bulk collected aboveground biomass for food preparation, smaller samples of Calluna and Molinia shoots were collected as well and stored for all individual plots and for both species separately for later analysis of elemental content. At the same day of collection, all obtained bulk vegetation samples were dried for 48 h at 60°C. Dried samples were then finely ground and homogenized at equal amounts based on dry mass. In order to produce food pellets that were easier to handle, the following procedure was done for each experimentally manipulated food. A 5 g l–1 agar solution was prepared in distilled water at 100°C. Each food stock was suspended in the agar solution after cooling of the solution to 60°C in batches of 15 g of stock in 45 ml solution. Each batch was transferred on a 10 × 10 cm patch of filter paper, spread out to a maximum thickness of 0.5 cm, cut in small squares of approximately 1 cm2 and subsequently dried in an oven for 48 h at 60°C. After drying, each food batch was separated from the filter paper and broken into pellets using the cut lines. Pellets from each treatment combination were stored at 5°C in closed jars.
Breeding Experiment
A climate controlled breeding chamber was used in which temperature was kept constant at 20°C and relative humidity at 55% (±5%). In order to reflect the light/dark regime in the field during the peak reproductive period of field crickets, artificial light regime was held constant at 16/8 h L/D. Field crickets were obtained in April 2014 by collecting penultimate and ultimate nymphal stages in the field in a nearby site within the National Park The Hoge Veluwe. Crickets were then transported to a climate controlled breeding chamber and transferred to separate plastic containers, which were randomly assigned to receive one of the four foods obtained from the field experiment. Containers (Figure 1) measured 17 × 10 × 7 cm (L × W × H). In each container, ca. 0.9 g of food pellets were offered on a 3 cm diameter petri dish. To allow for drinking, a water filled petri dish covered with gauze in which a piece of filter paper was stuck into the water was placed in each container. To allow for hiding, a tube of 4.7 cm length and 1.4 cm diameter covered in black plastic tape was placed in the container. To allow females to deposit eggs, each container was filled with 3 cm of sand which was kept moist by regular spraying with tap water. At intervals of (2)–3–4 days, water was refreshed and food pellets were replaced in each container. The remaining food pellets were often too fragmented to allow for weighing of residual food biomass. In general, food pellets were eaten to a high degree (90% or more of food consumed), unless pellets showed clear signs of fungal decay. Occurrence of fungal decay did not differ between treatments (32.6; 34.4; 31.0; and 35.5% for P+Ca+; P+Ca–; P–Ca+ and P– Ca–, respectively), and was often limited to one of the several food pellets offered. At inset of the experiment, and at each refreshment period, nymphs and adults were weighed on a microbalance. At the start of the experiment, 14 female nymphs were assigned to each food type, resulting in a total of 56 individuals. Following molt into the adult stage, females were kept from mating with males in a 2-week abstinence period in order to minimize prior field effects on nutrient and energy reserves and to allow for eggs to fully develop. After this period, males were placed with females in the same container in order to allow for mating to take place. Initially, males were kept with females constantly. However, in the first week of the experiment, a high degree of cannibalism by the females became apparent. Subsequently, males were kept with females for 1 day each week, which proved effective in eliminating the occurrence of cannibalism, while still allowing females to successfully mate (confirmed by the presence of spermatophores). During the reproductive period, containers were replaced at each of the 3–4 days refreshment intervals by new containers with a layer of fresh sand. Each rearing container was labeled by time period of female present, female identity and food treatment. Containers were subsequently incubated under the same climate regime and regularly sprayed with tap water to prevent dehydration of eggs. Containers were subsequently checked for hatched first instar nymphs at least once every 2 days. First hatchlings typically emerged two and a half weeks after replacement. Hatched nymphs were removed from each container using an insect exhauster and counted regularly until no nymphs were found for 17 consecutive days in the respective container. For each individual, date of mortality was recorded in order to obtain the length of the reproductive period for each individual. The feeding experiment ended when the last female died (88 days after adult eclosion), but counting of hatched nymphs continued for a longer period, by following the aforementioned procedure.
FIGURE 1
Food Quality
Collected Molinia and Calluna shoots were analyzed separately for elemental content. 200 mg of dried plant material was digested using 4 ml HNO3 (65%) and 1 ml H2O2 (35%) in teflon vessels using an Ethos D microwave (1200 MLS, Milestone, Sorisole, Italy), after which digests were analyzed for phosphorus (P), calcium (Ca), aluminum (Al), iron (Fe), potassium (K), magnesium, (Mg), manganese (Mn), sulfur (S), silicon (Si) and zinc (Zn) using an inductively coupled plasma emission spectrophotometer (ICP-OES; model IRIS Intrepid II XDL, Thermo Fisher Scientific, Franklin, United States). C and N contents of 3 mg homogenized dried plant material were determined by an elemental CNS analyzer (model NA 1500, Carlo Erba; Thermo Fisher Scientific, Franklin, United States). Food elemental content was calculated using the mean of Calluna + Molinia nutrient content, thus reflecting the food quality of the food mixture offered in the experiment.
Data Analysis
Breeding Experiment
Prior to analyzing the effect of treatments, we assessed whether other life-history parameters could have influenced reproduction, either dependent or independent of treatment. We analyzed whether initial weight, weight at start of mating, weight change and female reproductive lifespan (time between start of mating and death of the female in days) differed between food treatments, by fitting P×Ca, P+Ca, P, Ca and intercept only predictor models and compared these using the model AIC scores and AIC-derived model fit parameters (). We subsequently analyzed the model output of the minimum adequate model or second ranked model, in case the intercept only model was the minimum adequate model. Only for female reproductive lifespan, a significant negative treatment effect of P addition was apparent (Supplementary Figure 1 and Supplementary Tables 1, 2). Apart from life history parameters, we analyzed whether the observed incidences of cannibalism were influenced by the experimental food treatments. Based on the preliminary analyses, we used total number of hatched nymphs as well as daily reproduction, calculated as total number of hatched nymphs/female reproductive lifespan as measure of reproductive success. Since occurrence of cannibalism was significantly influenced by Ca+ treatment (see Results section “Occurrence of cannibalism”), we also included the occurrence of cannibalism as a covariate in the models, and present the results separately for all females and all non- cannibalistic females. Responses in reproductive parameters from cannibalistic females only are provided in Supplementary Figure 3. In assessing the effects of treatments on reproduction, we used an information theoretic approach by fitting generalized linear models varying between maximum (P × Ca × cannibalism) and minimum (intercept only) number of parameters and comparing the degree of fit of these models using the model AIC scores and AIC-derived model fit parameters (). In order to compare the effect of additional parameters on model outcome with the minimum adequate model, second ranking model outcomes are reported as well. A total overview of model ranking results is provided in Supplementary Table 3. In order to further assess the influence of cannibalism on model outcome, the same model procedures, excluding cannibalism as covariate, was performed on reproduction data from non-cannibalistic females. The results of this secondary analysis procedure are summarized in Supplementary Tables 4, 5. Total hatched eggs models were fitted using a negative binomial GLM with log link, daily reproduction models were fitted using a Gamma GLM with identity link.
Food Quality
Treatment effects on food elemental content and elemental ratios were analyzed using linear regression using an initial full factorial (P × Ca) model formulation. When no significant effect of the interaction term was apparent from this model, this parameter was dropped from the final model. An exception was made for plant tissue P, where a near significant effect for the P × Ca interaction was kept in the model in order to better compare model results with tissue N:P ratio (containing a significant P × Ca interaction term). In order to meet normality assumptions, elemental ratios were log-transformed prior to the analysis.
Comparisons of Food and Invertebrate Elemental Stoichiometry
We explored possible causal mechanisms underlying the observed responses in the feeding experiment by comparing data of invertebrate macro and micro-elemental content from literature (), a previously obtained dataset of elemental content of invertebrate taxa (Supplementary Table 6), and by comparing results with a heathland vegetation feeding experiment with Gryllus campestris performed earlier, which differed markedly in plant tissue N:P stoichiometry (for results and methods, see ).
Results
Breeding Experiment
Occurrence of Cannibalism
The frequency of occurrence of cannibalism in female crickets during the first week of allowed male mating was significantly higher in the lime addition (Ca+ treatment) groups compared to the Ca– groups (Figure 2, Ca+: 46.4%, Ca–: 14.3%, χ2 = 5.41, p = 0.02), but not in the P addition (P+ treatment) groups compared to the P– groups (Figure 2, P+: 35.7%, P–: 25%, χ2 = 0.34, p = 0.56).
FIGURE 2
Reproduction
Total number of hatched eggs was best explained by the model incorporating lime addition and occurrence of cannibalism as explaining variables (Table 1, AIC = 698.3). Based on these model parameters, female crickets that consumed their male in the first week of the experiment produced 177 eggs more than females that did not consume their male counterpart (Figures 3A,B and Table 1, intercept = 351 eggs; male eaten = 528 eggs, z = 3.12, p < 0.01). Liming significantly reduced the total number of hatched eggs by 107 eggs (Table 1 and Figures 3A,B, intercept = 351 eggs, Ca+ = 242, z = –3.1, p < 0.05). A model that also included P addition (P+ treatment) received less support (Table 1, AIC = 700.2, relative likelihood = 0.37), and the estimated effect of P addition in this model was minor and not significant, while the modeled effects of liming and cannibalism were highly comparable. Models excluding cannibalistic females show the similar results, with again the best performing model incorporating lime addition as explaining variable (Supplementary Table 4, AIC = 494.9), followed by a model also incorporating P addition (Supplementary Table 4, AIC = 496.9, relative likelihood = 0.37). Liming reduced total number of hatched eggs in non-cannibalistic females by 102 (Supplementary Table 5, intercept = 347, Ca+ = 245, z = –2.3, p < 0.05). P treatment only had a minor, non-significant effect on total hatched eggs in non-cannibalistic females (Supplementary Table 5).
TABLE 1
| Response parameter | Model | Parameter | Estimate | St. Error | z/t value |
| Total reproductive success (hatched eggs)a | Ca+ Male eaten | Intercept | 5.86 | 0.08 | 73.07*** |
| AIC = 698.3, params = 4, | Ca+ | –0.37 | 0.12 | –3.1** | |
| wi = 0.51 | Male eaten | 0.41 | 0.13 | 3.12** | |
| P+Ca+ Male eaten | Intercept | 5.86 | 0.1 | 60.81*** | |
| AIC = 700.2, params = 5, | P+ | –0.01 | 0.11 | –0.08NS | |
| ΔAIC = 1.99, wi = 0.19, | Ca+ | –0.37 | 0.12 | –3.11** | |
| rl = 0.37 | Male eaten | 0.41 | 0.13 | 3.1** | |
| Daily reproductive success (hatched eggs)/ | P+Ca | Intercept | 9.03 | 1.12 | 8.07*** |
| reproduction lifespan)b | AIC = 324.8, params = 4, | P+ | 3.27 | 1.25 | 2.62* |
| wi = 0.32 | Ca+ | –3.26 | 1.25 | –2.61* | |
| Ca+P+ Male eaten | Intercept | 8.88 | 1.15 | 7.71*** | |
| AIC = 325.9, params = 5, | P+ | 3.14 | 1.3 | 2.42* | |
| ΔAIC = 1.16, wi = 0.18, | Ca+ | –3.56 | 1.33 | –2.68** | |
| rl = 0.56 | Male eaten | 1.23 | 1.42 | 0.86NS | |
Model summaries for treatment effects on reproductive success parameters.
Minimum adequate models are presented first. For both response parameters, the second best performing model (selected using AIC and Akaike weight calculation) results are also presented under the minimum adequate model. P+: P addition. Ca+: lime treatment. Male eaten: female cannibalized on male counterpart. n = 55 for all models. NS, no significant effect.*, **, ***Significance p levels of 0.05; 0.01; and 0.001, respectively. aNegative binomial GLM with log link, estimates not back-transformed. bGamma GLM with identity link, estimates are number of hatched eggs/day. AIC, Akaike information criterion; wi, Akaike weight; ΔAIC, AIC difference between minimum adequate and second ranking model; rl, relative likelihood of second ranking model, calculated as 1/(wi(minimum adequate model)/wi(secondranking model)).
FIGURE 3
Daily reproduction was best explained by a model incorporating both P and lime addition as explanatory variables (Table 1, AIC = 324.8). Daily reproductive success was significantly increased by P addition (P+: Table 1 and Figure 3Cz = 2.62, p < 0.05) and decreased by lime addition (Ca+: Table 1 and Figure 3C, z = 2.61, p < 0.05) by about the same amount (3.3 eggs/day). A model that also included cannibalism received slightly less support (Table 1, AIC = 325.9, rl = 0.56). Cannibalism did not significantly influence daily reproductive success (compare Figures 3C,D), but improved the fit for lime addition as well as the estimated effect of this treatment.). Models excluding cannibalistic females show similar results with the best performing model incorporating P and lime addition as explaining variable (Supplementary Table 4, AIC = 237.4), followed by a model incorporating only P addition as explanatory variable (Supplementary Table 4, AIC = 239.1, relative likelihood = 0.42). Liming reduced daily reproduction in non-cannibalistic females by 3.2 eggs/day, while P addition increased daily reproduction by 3.4 eggs/day in both models (Supplementary Table 5). In these models, both treatment effects were only significant at the 0.1 threshold, probably due to the lower number of replicates and imbalanced grouping (i.e., more individuals in the Ca– groups).
Food Quality
As expected, P+ treatment significantly increased food P and significantly reduced food N:P ratio (Figure 4A and Table 2, +32.9 P, t = 9.22, p < 0.001, –21.9 N:P, i = –10.94, p < 0.001). This reduction in N:P ratio appeared to be less pronounced when lime was added, resulting in a significant interaction between P and lime addition treatments in N:P ratio (Figure 4A and Table 2, –19.4 N:P, t = 2.43, p < 0.05). Food N content was not significantly altered by either of the treatments (Table 2). Food N:P ratio was substantially higher than those reported for terrestrial invertebrates in the treatments without P addition (), whereas in treatments with P addition, plant tissue N:P closely matched this range (Figure 4A).
FIGURE 4
TABLE 2
| Estimate (SE) | t | Estimate (SE) | t | |
| Parameter | N | Fe | ||
| Intercept (P-Ca–) | 1,983.4 (26.8) | 74.0*** | 6.13 (0.78) | 7.86*** |
| P+ | –50.7 (32.3) | –1.57NS | –0.25 (0.90) | –0.28NS |
| Ca+ | 13.4 (32.3) | 0.42NS | –0.87 (0.90) | –0.97NS |
| P+ × Ca+ | – | – | – | – |
| P | N:Pa | |||
| Intercept (P-Ca–) | 41.43 (2.52) | 16.44*** | 3.87 (0.04) | 93.57*** |
| P+ | 32.87 (3.56) | 9.22*** | –0.61 (0.06) | –10.94*** |
| Ca+ | 4.42 (3.56) | 1.24NS | –0.13 (0.07) | –1.99NS |
| P+ × Ca+ | –10.12 (5.04) | –2.01⋅ | 0.22 (0.09) | 2.434* |
| Mg | Mn:Mga | |||
| Intercept (P-Ca–) | 56.00 (2.71) | 20.64*** | –3.30 (0.10) | –33.41*** |
| P+ | 10.06 (3.13) | 3.21** | -0.60 (0.11) | –5.29*** |
| Ca+ | 30.31 (3.13) | 3.13*** | –0.97 (0.11) | –8.55*** |
| P+ × Ca+ | – | – | – | – |
| Mn | Fe:Mga | |||
| Intercept (P-Ca–) | 2.01 (0.11) | 17.54*** | –2.27 (0.12) | –18.59*** |
| P+ | –0.60 (0.13) | –4.53*** | –0.19 (0.14) | –1.32NS |
| Ca+ | –0.72 (0.13) | –5.39*** | –0.56 (0.14) | –3.95** |
| P+ × Ca+ | – | – | – | – |
Model summaries for plant chemistry treatment effects.
P+: P addition. Ca+: lime treatment. P+ × Ca+: interaction between P and lime addition. n = 5 for all models. Model parameters and t-scores are estimates from the simplified model formulation: response ∼P × Ca.–: Term removed from model formulation during model simplification steps. NS, no significant effect.⋅, *, **, ***Significance p levels of 0.1, 0.05, 0.01, and 0.001, respectively. amodel performed on log transformed data in order to meet normal distribution assumption; estimates not back-transformed.
Both treatments influenced plant tissue elemental content of the food (Supplementary Figure 2). As expected, food Ca and Mg increased in lime treatments. P addition also increased food Ca and Mg content. Mn decreased strongly in lime treatments and to a lesser extent also in P addition treatment, while for Fe this was also reduced in lime treatments, although this was not significant (Supplementary Figure 2 and Table 2). Together, changes in food Mn, Fe and Mg concentrations resulted in pronounced changes in food Mn:Mg ratio (Figure 4B) and Fe:Mg ratio (Figure 4C), with lime treatment having a stronger or consistent negative effect (Figure 4B and Table 2: –0.023 Mn:Mg:, t = –8.08, p < 0.001; –0.04 Fe:Mg, t = p < 0.01) compared to P treatment (Figure 4B and Table 2: –0.017 Mn:Mg, t = –4.96, p < 0.001, –0.02 Fe:Mg, t = –1.32, p = 0.20). As a result, food Mn:Mg and Fe:Mg ratio declined to levels substantially lower than those observed for Ensifera and other invertebrates (Figures 4B,C).
Comparison of Breeding Experiments
When comparing the results of our hatching experiment with results obtained in a previous feeding experiment, there appears to be a non-linear response in the number of hatched eggs to tissue N:P ratio of the plants fed to the crickets (Figure 5). At N:P molar ratio below 50, the number of hatched eggs varied around 300 for a large range of N:P ratios, but at ratios above 50 (i.e., the Strabrecht “heath” treatment), there was a substantial reduction in the mean number of hatched eggs (107.3). Furthermore, at the lowest treatment N:P ratio of 26.2 (current experiment; P+Ca– treatment), mean number of hatched eggs after 30 days was substantially higher (396.6) compared to treatments with intermediate food N:P ratios (282.2–333).
FIGURE 5

Comparison of reproductive potential of female Gryllus campestris with mean plant tissue N:P ratio under two feeding experiments. Blues: current experiment (by experimental treatment). Greens: experiment performed with plant material collected from Strabrechtse Heide heathland under contrasting nutrient availability (
Discussion
P Amelioration Effects
P amelioration successfully achieved the aim of removing stoichiometric imbalances, with food N:P ratio decreasing to values that align with the mean ratio of terrestrial invertebrates. Daily reproductive success also improved when crickets were fed food from P addition plots, with female crickets being able to reach their maximal reproductive potential in a shorter timespan. However, higher daily reproduction did not translate in a higher total reproductive success. In allowing crickets to produce eggs until death, those who were fed higher (i.e., less favorable) N:P food were able to eventually catch up with low (i.e., more favorable) N:P food fed crickets, having access to more food over this protracted period. The faster rate of egg production in individuals fed with food from P addition plots may be linked to an increased potential for protein synthesis, as ribosomal RNA is the major pool of cellular P (
Liming Effects
Contrary to our expectations, the liming treatment did not improve but rather resulted in a significant reduction in reproductive performance. We also observed increased cannibalistic behavior of female crickets fed with food from limed plots, suggesting that the lower reproductive success was related to a nutrient deficiency induced by liming. For total reproduction, a significant increase in eggs produced under cannibalism partly compensated for the decrease in eggs produced when fed food from limed plots, but even when considering only cannibalistic females (Supplementary Figure 3), there was a clear reduction in total reproduction with liming. Since cannibalism affected total reproduction, but not daily reproduction (for which P addition and subsequent N:P ratio was found to be important), our results suggest that cannibalism was induced due to a deficiency in other (trace) elements rather than a deficiency in CNP macronutrients.
Liming resulted in substantial shifts in food Mn:Mg and Fe:Mg ratios. The skewed ratios arose because plant tissue Mg increased, while Mn and Fe decreased, with Fe decreasing only in the combined P and lime addition treatment. Decreased plant tissue Mn as a result of liming has been reported earlier (
Both Mn and Fe deficiencies can be linked to reduced insect performance. In invertebrates, uptake of Mn2+ and Fe2+ is partially regulated through the divalent metal transporter (DMT1) ortholog malvolio (
Cannibalism and Feeding Limitations
The occurrence of cannibalism in the first week of the experiment was an unexpected outcome of the experiment. In nature, cannibalism is rare (T. Tregenza, pers. comm.1), but when males and females were kept within a confined space, females were apparently able to predate on males. Although cannibalism complicated the analysis and interpretation of our results, the higher prevalence of cannibalism in the lime addition treatment provided us with an additional clue to the effect of liming on reproduction. Including cannibalism in our models improved the model fit and revealed a more negative effect of liming on reproduction (as cannibalism enabled females to partly compensate the loss in fecundity due to liming). The results are not confounded by cannibalism, since models excluding cannibalistic females showed the same outcome, albeit with a lower statistical power due to large reductions in the replicates for some treatments.
In our study female crickets were not allowed to actively select food of differing quality. This deprived the females the option to regulate nutrient intake. It is well established that insects are able to complement low intake of specific nutrients by selecting food sources rich in this nutrient. Such regulation in insects and crickets has been extensively reported by several authors for protein and carbohydrates (
In the experimental plots used as food source for this experiment, P and lime addition treatments resulted in a significant increase in plant species richness compared to control (
Management Implications
Our results demonstrate that high-impact, large scale management practices such as sod cutting to reduce the effects of chronic increased nitrogen and acidifying deposition include a high risk of unwanted side effects for invertebrate populations. The nature of these unwanted side effects lies in disturbance, habitat homogenization and destruction and in increasing stoichiometric imbalances as shown in this and earlier studies (
Statements
Data availability statement
The raw data supporting the conclusions of this article has been made publicly available via the DANS easy repository: https://doi.org/10.17026/dans-zsa-f3y9.
Author contributions
JV and HS conceived the ideas and designed methodology. JV, MW, and JK collected the data. JV and JK performed the breeding experiment. JV and WV analyzed the data and drafted the manuscript. All authors contributed to the drafts and gave final approval for publication.
Funding
This research was funded by the O+BN program, financed by the Dutch Ministry of Agriculture, Nature and Food Quality and BIJ12 (Grant No. OBN207-DZ).
Acknowledgments
We thank Evi Verbaarschot for assistance with field work, Jan-Willem Wolters for laboratory analyses of invertebrate elemental content and Jakob Leidekker from the National Park De Hoge Veluwe for providing access to the park and allowing us to perform experimental research.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2021.659363/full#supplementary-material
Footnotes
References
1
AuC.BenedettoA.AschnerM. (2008). Manganese transport in eukaryotes: the role of DMT1.NeuroToxicology29569–576. 10.1016/j.neuro.2008.04.022
2
BehmerS. T. (2009). Insect herbivore nutrient regulation.Annu. Rev. Entomol.54165–187. 10.1146/annurev.ento.54.110807.090537
3
Ben-ShaharY. (2018). The impact of environmental Mn exposure on insect biology.Front. Genet.9:70. 10.3389/fgene.2018.00070
4
BernerD.BlanckenhornW. U.KornerC. (2005). Grasshoppers cope with low host plant quality by compensatory feeding and food selection: N limitation challenged.Oikos111525–533.
5
BobbinkR.HeilG. W. (1993). “Atmospheric deposition of sulphur and nitrogen in heathland ecosystems,” in Heathlands: Patterns and Processes in a Changing Environment, edsAertsR.HeilG. W. (Dordrecht: Kluwer Academic Publishers), 25–50.
6
BobbinkR.HettelinghJ. P. (2011). “Review and revision of empirical critical loads and dose-response relationships,” in Proceedings of an Expert Workshop, Noordwijkerhout, 23-25 June 2010, (Bilthoven: Coordination Centre for Effects).
7
BobbinkR.HicksK.GallowayJ.SprangerT.AlkemadeR.AshmoreM.et al (2010). Global assessment of nitrogen deposition effects on terrestrial plant diversity: a synthesis.Ecol. Appl.2030–59. 10.1890/08-1140.1
8
BobbinkR.HornungM.RoelofsJ. G. M. (1998). The effects of air-borne nitrogen pollutants on species diversity in natural and semi-natural European vegetation.J. Ecol.86717–738.
9
BoersmaM.ElserJ. J. (2006). Too much of a good thing: on stoichiometrically balanced diets and maximal growth.Ecology871325–1330.
10
BorerE. T.LindE. M.FirnJ.SeabloomE. W.AndersonT. M.BakkerE. S.et al (2019). More salt, please: global patterns, responses and impacts of foliar sodium in grasslands.Ecol. Lett.221136–1144. 10.1111/ele.13270
11
BurnhamK. P.AndersonD. R. (2002). Model Selection and Multimodel Inference : A Practical Information-Theoretic Approach.New York, NY: Springer.
12
CeaseA. J.FayM.ElserJ. J.HarrisonJ. F. (2016). Dietary phosphate affects food selection, post-ingestive phosphorus fate, and performance of a polyphagous herbivore.J. Exp. Biol.21964–72. 10.1242/jeb.126847
13
De GraafM. C. C.BobbinkR.SmitsN. A. C.Van DiggelenR.RoelofsJ. G. M. (2009). Biodiversity, vegetation gradients and key biogeochemical processes in the heathland landscape.Biol. Conserv.1422191–2201.
14
De GraafM. C. C.VerbeekP. J. M.BobbinkR.RoelofsJ. G. M. (1998). Restoration of species-rich dry heaths: the importance of appropriate soil conditions.Acta Bot. Neerl.4789–111.
15
DiemontW. H. (1996). Survival of Dutch Heathlands. “IBN Scientific Contributions 1.”.Wageningen: DLO Institute for Forestry and Nature Research (IBN-DLO).
16
DorlandE.BobbinkR.MesselinkJ. H.VerhoevenJ. T. A. (2003). Soil ammonium accumulation after sod cutting hampers the restoration of degraded wet heathlands.J. Appl. Ecol.40804–814. 10.1046/j.1365-2664.2003.00845.x
17
DorlandE.van den BergL. J. L.van de BergA. J.VermeerM. L.RoelofsJ. G. M.BobbinkR. (2004). The effects of sod cutting and additional liming on potential net nitrification in heathland soils.Plant Soil265267–277.
18
ElserJ. J.DobberfuhlD. R.MacKayN. A.SchampelJ. H. (1996). Organism size, life history, and N:P Stoichiometry: toward a unified view of cellular and ecosystem processes.BioScience46674–684. 10.2307/1312897
19
ElserJ. J.FaganW. F.DennoR. F.DobberfuhlD. R.FolarinA.HubertyA.et al (2000a). Nutritional constraints in terrestrial and freshwater food webs.Nature408578–580.
20
ElserJ. J.SternerR. W.GorokhovaE.FaganW. F.MarkowT. A.CotnerJ. B.et al (2000b). Biological stoichiometry from genes to ecosystems.Ecol. Lett.3540–550. 10.1046/j.1461-0248.2000.00185.x
21
FilipiakM.KuszewskaK.AsselmanM.DenisowB.StawiarzE.WoyciechowskiM.et al (2017). Ecological stoichiometry of the honeybee: pollen diversity and adequate species composition are needed to mitigate limitations imposed on the growth and development of bees by pollen quality.PLoS One12:e0183236. 10.1371/journal.pone.0183236
22
FottnerS.HardtleW.NiemeyerM.NiemeyerT.von OheimbG.MeyerH.et al (2007). Impact of sheep grazing on nutrient budgets of dry heathlands.Appl. Veg. Sci.10391–398. 10.1111/j.1654-109X.2007.tb00438.x
23
GeiserD. L.WinzerlingJ. J. (2012). Insect transferrins: multifunctional proteins.Biochim. Biophys. Acta1820437–451. 10.1016/j.bbagen.2011.07.011
24
GjengedalE.MartinsenT.SteinnesE. (2015). Background levels of some major, trace, and rare earth elements in indigenous plant species growing in Norway and the influence of soil acidification, soil parent material, and seasonal variation on these levels.Environ. Monit. Assess.187:386. 10.1007/s10661-015-4479-5
25
GravelandJ. (1998). Effects of acid rain on bird populations.Environ. Rev.641–54. 10.1139/a98-003
26
GravelandJ.vanderWalR. (1996). Decline in snail abundance due to soil acidification causes eggshell defects in forest passerines.Oecologia105351–360. 10.1007/bf00328738
27
GravelandJ.VanderwalR.VanbalenJ. H.VannoordwijkA. J. (1994). Poor reproduction in forest passerines from decline of snail abundance on acidified soils.Nature368446–448. 10.1038/368446a0
28
HaddadN. M.HaarstadJ.TilmanD. (2000). The effects of long-term nitrogen loading on grassland insect communities.Oecologia12473–84. 10.1007/s004420050026
29
HärdtleW.NiemeyerM.NiemeyerT.AssmannT.FottnerS. (2006). Can management compensate for atmospheric nutrient deposition in heathland ecosystems?J. Appl. Ecol.43759–769. 10.1111/j.1365-2664.2006.01195.x
30
HärdtleW.Von OheimbG.GerkeA. K.NiemeyerM.NiemeyerT.AssmannT.et al (2009). Shifts in N and P budgets of heathland ecosystems: effects of management and atmospheric inputs.Ecosystems12298–310. 10.1007/s10021-008-9223-3
31
HarrisonS. J.RaubenheimerD.SimpsonS. J.GodinJ. G. J.BertramS. M. (2014). Towards a synthesis of frameworks in nutritional ecology: interacting effects of protein, carbohydrate and phosphorus on field cricket fitness.Proc. R. Soc. B Biol. Sci.281:20140539. 10.1098/rspb.2014.0539
32
HaynesR. (1982). Effects of liming on phosphate availability in acid soils.Plant Soil68289–308. 10.1007/bf02197935
33
IngrischS.KöhlerG. (1998). Die Heuschrecken Mitteleuropas.Magdeburg: Westarp Wissenschaften.
34
JoernA.ProvinT.BehmerS. T. (2011). Not just the usual suspects: insect herbivore populations and communities are associated with multiple plant nutrients.Ecology931002–1015. 10.1890/11-1142.1
35
JonesL.StevensC.RoweE. C.PayneR.CapornS. J. M.EvansC. D.et al (2017). Can on-site management mitigate nitrogen deposition impacts in non-wooded habitats?Biol. Conserv.212464–475. 10.1016/j.biocon.2016.06.012
36
KappesH.ToppW. (2014). Responses of forest snail assemblages to soil acidity buffer system and liming.Ecol. Res.29757–766. 10.1007/s11284-014-1167-9
37
KleijnD.BekkerR. M.BobbinkR.De GraafM. C. C.RoelofsJ. G. M. (2008). In search for key biogeochemical factors affecting plant species persistence in heathland and acidic grasslands: a comparison of common and rare species.J. Appl. Ecol.45680–687. 10.1111/j.1365-2664.2007.01444.x
38
KleukersR. M. J. C.van NieukerkenE. J.OdéB.WillemseL. P. M.van WingerdenW. K. R. E.van NieukerkenE. J.et al (1997). De Sprinkhanen en Krekels van Nederland (Orthoptera).Utrecht: KNNV.
39
KorenkoS.KulaE.HolecM.JarabM.MichalkovaV. (2008). Influence of liming on the epigeic spider (Araneae) community of the Krusne hory Mts (Czech Republic).Eur. J. Soil Biol.44559–566. 10.1016/j.ejsobi.2008.08.001
40
LockeM.NicholH. (1992). Iron economy in insects - transport, metabolism, and storage.Annu. Rev. Entomol.37195–215.
41
MaesD.DecleerK.De KeersmaekerL.Van UytvanckJ.LouetteG. (2017). Intensified habitat management to mitigate negative effects of nitrogen pollution can be detrimental for faunal diversity: a comment on Jones et al. (2017).Biol. Conserv.212493–494. 10.1016/j.biocon.2017.03.001
42
MandilarasK.PathmanathanT.MissirlisF. (2013). Iron absorption in Drosophila melanogaster.Nutrients51622–1647. 10.3390/nu5051622
43
McKieB. G.PetrinZ.MalmqvistB. (2006). Mitigation or disturbance? Effects of liming on macroinvertebrate assemblage structure and leaf-litter decomposition in the humic streams of northern Sweden.J. Appl. Ecol.43780–791. 10.1111/j.1365-2664.2006.01196.x
44
MichalkoR.KulaE.KosulicO. (2018). Liming alters body size distribution in a community of epigeic spiders in birch forest (Betula pendula Roth).Ann. For. Sci.75:92. 10.1007/s13595-018-0769-8
45
MillaleoR.Reyes- DiazM.IvanovA. G.MoraM. L.AlberdiM. (2010). Manganese as essential and toxic element for plants: transport, accumulation and resistance mechanisms.J. Soil Sci. Plant Nutr.10470–481.
46
NicholH.LawJ. H.WinzerlingJ. J. (2002). Iron metabolism in insects.Annu. Rev. Entomol.47535–559. 10.1146/annurev.ento.47.091201.145237
47
NijssenM. E.WallisDeVriesM. F.SiepelH. (2017). Pathways for the effects of increased nitrogen deposition on fauna.Biol. Conserv.212423–431. 10.1016/j.biocon.2017.02.022
48
OrgadS.NelsonH.SegalD.NelsonN. (1998). Metal ions suppress the abnormal taste behavior of the Drosophila mutant malvolio.J. Exp. Biol.201(Pt. 1), 115–120.
49
PabianS. E.BrittinghamM. C. (2007). Terrestrial liming benefits birds in an acidified forest in the northeast.Ecol. Appl.172184–2194. 10.1890/07-0394.1
50
PabianS. E.ErmerN. M.TzilkowskiW. M.BrittinghamM. C. (2012). Effects of liming on forage availability and nutrient content in a forest impacted by acid rain.PLoS One7:e39755. 10.1371/journal.pone.0039755
51
PenuelasJ.JanssensI. A.CiaisP.ObersteinerM.SardansJ. (2020). Anthropogenic global shifts in biospheric N and P concentrations and ratios and their impacts on biodiversity, ecosystem productivity, food security, and human health.Glob. Change Biol.261962–1985. 10.1111/gcb.14981
52
PoteatM. D.Díaz-JaramilloM.BuchwalterD. B. (2012). Divalent metal (Ca, Cd, Mn, Zn) uptake and interactions in the aquatic insect Hydropsyche sparna.J. Exp. Biol.2151575–1583. 10.1242/jeb.063412
53
RaubenheimerD.SimpsonS. J. (1993). The geometry of compensatory feeding in the locust.. Behav.45953–964. 10.1006/anbe.1993.1114
54
ReddyS. R. R.CampbellJ. W. (1969). Arginine metabolism in insects: properties of insect fat body arginase.Comp. Biochem. Physiol.28515–534. 10.1016/0010-406X(69)92087-8
55
RoemW. J.BerendseF. (2000). Soil acidity and nutrient supply ratio as possible factors determining changes in plant species diversity in grassland and heathland communities.Biol. Conserv.92151–161. 10.1016/s0006-3207(99)00049-x
56
RoemW. J.KleesH.BerendseF. (2002). Effects of nutrient addition and acidification on plant species diversity and seed germination in heathland.J. Appl. Ecol.39937–948.
57
ShawkiA.MackenzieB. (2010). Interaction of calcium with the human divalent metal-ion transporter-1.Biochem. Biophys. Res. Commun.393471–475. 10.1016/j.bbrc.2010.02.025
58
SiepelH.BobbinkR.van de RietB. P.van den BurgA. B.JongejansE. (2019). Long-term effects of liming on soil physico-chemical properties and micro-arthropod communities in Scotch pine forest.Biol. Fertil. Soils55675–683. 10.1007/s00374-019-01378-3
59
SiepelH.VogelsJ.BobbinkR.BijlsmaR.-J.JongejansE.de WaalR.et al (2018). Continuous and cumulative acidification and N deposition induce P limitation of the micro-arthropod soil fauna of mineral-poor dry heathlands.Soil Biol. Biochem.119128–134. 10.1016/j.soilbio.2018.01.025
60
SternerR. W.ElserJ. J. (2002). Ecological Stoichiometry: The Biology of Elements from Molecules to the Biosphere.New Jersey, NJ: Princeton University Press.
61
StevensC. J.DiseN. B.GowingD. J. (2009). Regional trends in soil acidification and exchangeable metal concentrations in relation to acid deposition rates.Environ. Pollut.157313–319. 10.1016/j.envpol.2008.06.033
62
StevensC. J.DuprèC.DorlandE.GaudnikC.GowingD. J. G.BleekerA.et al (2010). Nitrogen deposition threatens species richness of grasslands across Europe.Environ. Pollut.1582940–2945. 10.1016/j.envpol.2010.06.006
63
TianD.NiuS. (2015). A global analysis of soil acidification caused by nitrogen addition.Environ. Res. Lett.10:024019. 10.1088/1748-9326/10/2/024019
64
TylerG.OlssonT. (2001). Plant uptake of major and minor mineral elements as influenced by soil acidity and liming.Plant Soil230307–321. 10.1023/A:1010314400976
65
VisanuvimolL.BertramS. M. (2011). How dietary phosphorus availability during development influences condition and life history traits of the cricket, Acheta domesticus.J. Insect Sci.11:63.
66
VogelsJ. J.VerberkW. C. E. P.LamersL. P. M.SiepelH. (2017). Can changes in soil biochemistry and plant stoichiometry explain loss of animal diversity of heathlands?Biol. Conserv.212432–447. 10.1016/j.biocon.2016.08.039
67
VogelsJ. J.WebbN. R.SiepelH. H. (2013). “Impact of changed plant stoichiometric quality on heathland fauna composition,” in Economy and Ecology of Heathlands, edsDiemontW. H.HeijmanW. J. M.SiepelH. H.WebbN. R. (Zeist: KNNV publishing).
68
VogelsJ. J.WeijtersM. J.BobbinkR.BijlsmaR.-J.LamersL. P. M.VerberkW. C. E. P.et al (2020). Barriers to restoration: soil acidity and phosphorus limitation constrain recovery of heathland plant communities after sod cutting.Appl. Veg. Sci.2394–106. 10.1111/avsc.12471
69
WaldbauerG. P.FriedmanS. (1991). Self-selection of optimal diets by insects.Annu. Rev. Entomol.3643–63.
70
WebbN. R. (1998). The traditional management of European heathlands.J. Appl. Ecol.35987–990.
Summary
Keywords
acidification, restoration management, nitrogen deposition, elemental ecology, invertebrate diversity, insect decline, sod-cutting, ecological stoichiometry
Citation
Vogels JJ, Verberk WCEP, Kuper JT, Weijters MJ, Bobbink R and Siepel H (2021) How to Restore Invertebrate Diversity of Degraded Heathlands? A Case Study on the Reproductive Performance of the Field Cricket Gryllus campestris (L.). Front. Ecol. Evol. 9:659363. doi: 10.3389/fevo.2021.659363
Received
27 January 2021
Accepted
22 April 2021
Published
19 May 2021
Volume
9 - 2021
Edited by
Anita Diaz, Bournemouth University, United Kingdom
Reviewed by
Tom Tregenza, University of Exeter, United Kingdom; Piotr Sewerniak, Nicolaus Copernicus University in Toruń, Poland
Updates

Check for updates
Copyright
© 2021 Vogels, Verberk, Kuper, Weijters, Bobbink and Siepel.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joost J. Vogels, j.vogels@science.ru.nl
This article was submitted to Conservation and Restoration Ecology, a section of the journal Frontiers in Ecology and Evolution
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.