Introduction {#s1}
============

The ecological success and life history of organisms is constrained by a variety of intrinsic and extrinsic factors [@pone.0103697-Stearns1]. For herbivores, diet quality is believed to be one of the most influential extrinsic conditions [@pone.0103697-Chapman1]. Increasingly, nutritional studies have demonstrated that dietary nutrient ratios (*e.g*. N:P) have major effects on growth, reproductive success, and other performance characters in animals [@pone.0103697-Elser1]. These studies indicate that nutrient balance may be equally or perhaps even more important than the absolute amount of any one nutrient [@pone.0103697-Elser1], [@pone.0103697-Simpson1]. These findings are central to the theory of ecological stoichiometry, which focuses on the balance of chemical elements and energy in living systems from molecular to global scales [@pone.0103697-Sterner1]--[@pone.0103697-Anderson1], and the related framework, nutritional geometry [@pone.0103697-Simpson1]. While several recent studies have demonstrated that stoichiometric homeostasis may be important in influencing ecological dynamics in various systems [@pone.0103697-Sterner1], [@pone.0103697-Persson1], [@pone.0103697-Yu1], and researchers have identified some important pre-ingestive and post-ingestive regulatory mechanisms in insects to achieve their nutritional balances [@pone.0103697-Behmer1], information is still not sufficient on how organisms achieve stoichiometric homeostasis in the face of elementally imbalanced foods. This is especially true for terrestrial consumers under realistic field conditions.

Body N:P ratio has been considered as a central parameter linking key aspects of a species, such as its structural investments, growth rate, and biogeography [@pone.0103697-Elser3], [@pone.0103697-Zhang1]. Indeed, compared to C:N and C:P ratios, N:P ratios of a species usually have a considerably narrower range of variation [@pone.0103697-Elser4]--[@pone.0103697-Zhang2]. Importantly, stoichiometric mismatches between consumers and plants due to higher C:nutrient ratios in plants could result in N or P limitation of animal performance, depending on the animal\'s N:P ratio [@pone.0103697-Joern1], [@pone.0103697-Hessen1]; thus, regulation of N and P becomes a central issue for herbivorous animals experiencing imbalanced chemical environments [@pone.0103697-Yu1], [@pone.0103697-Behmer1]. However, N and P differ in their functions and cycling pathways not only within ecosystems but also within organisms [@pone.0103697-Sterner1]. For example, previous work has demonstrated that excess nitrogen eaten by herbivores may be stored within the organisms or excreted as uric acid or in the volatile form of ammonium, whereas P remains in phosphate form [@pone.0103697-Simpson2]--[@pone.0103697-Han1]. Furthermore, the relative importance of frass excretion of N and P for regulating N:P stoichiometry in insects remains unclear. We also do not know how the intensity and direction of ecosystem nutrient regulation are affected by grasshopper density, an insight that may be important in understanding when and how insect outbreaks occur.

Grasshoppers are a major component in grassland food webs and are important for controlling plant biodiversity and productivity in grassland ecosystems [@pone.0103697-Chapman1], [@pone.0103697-Ritchie1], [@pone.0103697-Belovsky1]. Similar to most animals, grasshoppers exhibit greater stoichiometric homeostasis compared to their food plants [@pone.0103697-Zhang3]. Herbivores and their food plants can have interactive effects on each other\'s C:N:P ratios through various direct and indirect mechanisms [@pone.0103697-Elser6]. Indeed, our recent work has shown that increasing density in grasshoppers reduces plant standing biomass and dramatically changes the C:N:P stoichiometry of both plants and litterfall [@pone.0103697-Zhang2], possibly due to effects of food plant preference and availability [@pone.0103697-Elser6]. However, we do not yet know if such density-dependent effects of grasshopper herbivory on plant chemical stoichiometry feedback on grasshopper stoichiometry.

The objective of this study was to answer the following questions: 1) Does N:P stoichiometry in grasshoppers change in response to their own density and plant N:P ratios? If not, then 2) How do grasshoppers achieve this N:P stoichiometric homeostasis? Finally, we assessed the implications of these findings in the context of the overall dynamics of grasshopper and grassland plant interactions.

Materials and Methods {#s2}
=====================

Ethics Statement {#s2a}
----------------

No specific permissions were required for all these locations, for these lands were not private and the study did not involve endangered or protected species.

Study site {#s2b}
----------

The field manipulation experiment was carried out in a temperate steppe in the Xilin River Basin in northern China (43°33\' N,116°40\' E, 1,250 m a.s.l.), which is located at the southeastern part of the Eurasian Steppe. The mean annual temperature in this area is 0.3 °C with mean monthly temperatures ranging from −21.6°C in January to 19.0 °C in July. The mean annual precipitation is 346 mm, more than 80% of which occurs during the growing season from May to September. The soil type is dark chestnut (namely Calcic Luvisols in the FAO classification system) with an average depth of ∼100 cm [@pone.0103697-Zhang2]. Vegetation at the study site is dominated by *Leymus chinensis* and *Stipa grandis* and accompanied by other common species such as *Agropyron cristatum*, *Achnatherum sibiricum*, *Cleistogenes squarrosa*, *Carex korshinskyi*, *Koeleria cristata* and *Allium* sp. Aboveground biomass of these seven species accounts for more than 80% of the total community primary production.

Study organism {#s2c}
--------------

The grasshopper *Oedaleus asiaticus* (formerly *O. decorus asiaticus*), a member of the subfamily Oedipodinae (Orthoptera: Acrididae: Oedipodinae), is a dominant herbivorous insect in the study area, with a population density ranging from \<10 individuals m^−2^ before 1990 [@pone.0103697-Kang1] to \>50 individuals m^−2^ after 1990 [@pone.0103697-Xu1]. It is a mid-sized Acridid species (2.64±0.10 cm in length for males, 3.56±0.24 cm for females). *O. asiaticus* has considerable migratory capacity in its brown-colored phase, making it the dominant locust of north Asian grasslands [@pone.0103697-Cease1]. Male and female adults of *O. asiaticus* can be visually identified [@pone.0103697-Zhang2].

Field experiment and sample collecting {#s2d}
--------------------------------------

During July and August 2007 we established 72 cages within a long-term experimental plot (fenced in 2005 in order to exclude sheep grazing) of the Inner Mongolia Grassland Ecosystem Research Station, Chinese Academy of Sciences. The distance between cages was 1 m. Each cage consisted of a steel-framed and nylon-netted cube (1×1×1 m) with the lower edge sealed against the soil to prevent escape of experimental animals and intrusion of other species. A 0.5×0.5 m zippered window on the upper side of each cage allowed experimental operations such as introduction of grasshoppers, checking of grasshopper densities, and removal of predators (mostly spiders).

We designed a density gradient experiment to simulate grasshopper outbreaks with different intensities. On 20 July 2007, adult *O. asiaticus* were introduced to cages at 6 treatment levels: 0 (control), 2, 4, 8, 16, 32 grasshoppers m^−2^. The sex ratio of each treatment was 1∶1 (male to female, consistent with its natural sex ratio). There were 12 replicates for each treatment, the experiment was based on orthogonal design and had two set of the same number (6×6) of cages. After the experimental grasshoppers were introduced, we checked each cage every three days to confirm the number of living grasshoppers (dead insects were removed and replaced with a similarly-sized animal of the same sex), repair damaged cages, and eliminate spiders and other species of grasshoppers that occurred naturally within the cages. Grasshoppers introduced at the very beginning of grazing were labeled in the pronotum with red paint in order to ensure the sampled insects were those that lived in the cages at a given density condition for 30 days.

On 20 August 2007 (one month after the start), cages were removed. The aboveground standing biomass was harvested by clipping off all plants at the soil surface. According to the methods reliable in our previous work, plants were sorted by species, grasshoppers with red paint markers were carefully removed according to sex and then kept in small cages in lab without feeding for 24 h in order to collect their frass, which was then immediately dried. After that, grasshoppers of both sexes were dried and weighed individually. Each sample of grasshopper used for further chemical analysis was composed of male and female insects with the sex ratio of 1∶1 (male to female). All samples (grasshoppers and plants) were oven dried at 65°C for 48 h and then milled to a fine powder and screened with 0.1-mm mesh before analyzing their N and P contents.

Chemical Analyses {#s2e}
-----------------

Total N of plant and insect samples was determined by using a micro-Kjeldahl procedure with sulfuric acid and digestion catalyst. Total phosphorus was analyzed colorimetrically by applying the ammonium molybdate method [@pone.0103697-Bao1]. All element concentrations in this work were expressed as g kg^−1^ (dry material), and elemental stoichiometry were molar ratios. In the mean time, all analyses are for samples pooled at the plot level.

Statistical analyses {#s2f}
--------------------

SPSS version 11.5 (SPSS) was used to analyze our experimental data. T-tests and one-way ANOVAs were applied to compare the N:P stoichiometric difference between different treatments. The relationships between each couple of variables were examined by calculating the Pearson Coefficient, and two-tailed tests of significances were used. Regression analysis was employed to estimate the curves fitted to each couple of experimental variables.

Results {#s3}
=======

Relationships between N and P concentrations in grasshopper body, frass, and food plants {#s3a}
----------------------------------------------------------------------------------------

N and P concentrations in grasshopper bodies were not correlated. In contrast, significant positive linear relationships were found between N and P concentrations in grasshopper frass (*n* = 44, *r* = 0.569, *P*\<0.001) and in food plants (*n* = 44, *r* = 0.41, *P* = 0.002) ([Fig. 1](#pone-0103697-g001){ref-type="fig"}).

![Correlations between N and P concentrations in grasshopper body (a), frass (b), and food plants (c) of grasshopper *Oedaleus asiaticus.*](pone.0103697.g001){#pone-0103697-g001}

Contrasting N:P stoichiometric responses in grasshopper body and frass to grasshopper density {#s3b}
---------------------------------------------------------------------------------------------

Over the density treatments, *O. asiaticus* adults maintained consistent N and P concentrations and N:P ratios, showing strong N:P homeostasis ([Fig. 2d--f](#pone-0103697-g002){ref-type="fig"} and [Table 1](#pone-0103697-t001){ref-type="table"}). Along the experimental gradient of grasshopper density, adult *O. asiaticus* N (*n* = 5, *r* = 0.421, *P* = 0.480) and P (*n* = 5, *r* = 0.685, *P* = 0.202) concentrations as well as N:P ratios (*n* = 5, *r* = −0.545, *P* = 0.342) remained similar ([Fig. 2d--f](#pone-0103697-g002){ref-type="fig"}). However, frass P concentrations declined as grasshopper density increased, and the correlation between frass P and grasshopper density was marginally significant (*n* = 5, *r* = −0.809, *P* = 0.097) at treatment-level ([Fig. 2h](#pone-0103697-g002){ref-type="fig"}) and highly significant (*n* = 45, *r* = −0.467, *P* = 0.001) at quadrat (cage) level ([Fig. S1b](#pone.0103697.s001){ref-type="supplementary-material"}), and further analysis show that phosphorus decline in frass was occurred when grasshopper density decreased from 2 to 8 individuals per square meter ([Fig. S2d--f](#pone.0103697.s002){ref-type="supplementary-material"}). In contrast to changes in P, N concentrations of frass did not change appreciably ([Fig. 2g](#pone-0103697-g002){ref-type="fig"} and [Fig. S1a](#pone.0103697.s001){ref-type="supplementary-material"}). Reflecting the changes in N and P contents, frass N:P ratios (*n* = 45, *r* = 0.474, *P* = 0.001) increased with grasshopper densities ([Fig. S1c](#pone.0103697.s001){ref-type="supplementary-material"}). For food plants, N (*n* = 6, *r* = −0.963, *P* = 0.002) and P (*n* = 6, *r* = −0.925, *P* = 0.008) concentrations both decreased significantly and to similar extents as grasshopper density increased, thus, N:P ratios of food plants (*n* = 6, *r* = −0.358, *P* = 0.485) did not exhibit a strong trend with grasshopper density ([Fig. 2a--c](#pone-0103697-g002){ref-type="fig"}).

![Responses of N and P contents and N:P stoichiometry, of food plants (a--c), grasshopper body (d--f) and frass (g--i), to increasing grasshopper density.\
Error bars indicate ±1 SE. In (a), a strong negative relationship was found between host plant N concentrations and grasshopper density (ANOVA: *r^2^* = 0.93, *F* = 51.16, *P* = 0.002) and yielded the following equation: y = −0.13 x+14.52. In (b), a strong negative relationship was found between host plant N concentrations and grasshopper density (ANOVA: *r^2^* = 0.86, *F* = 23.59, *P* = 0.008) and yielded the following equation: y = −0.01 x+1.17. In (h), a marginally negative relationship was found between frass P concentrations and grasshopper density (ANOVA: *r^2^* = 0.65, *F* = 5.69, *P* = 0.097) and yielded the following equation: y = −0.025 x+2.14.](pone.0103697.g002){#pone-0103697-g002}

10.1371/journal.pone.0103697.t001

###### Results (*F* and *P* values) of one-way ANOVAs on the effects of grasshopper density on N:P ratio of food plant, grasshopper body and frass, corresponding to [Figure 2](#pone-0103697-g002){ref-type="fig"}.

![](pone.0103697.t001){#pone-0103697-t001-1}

  Response variable       Driver: Grasshopper density           
  ---------------------- ----------------------------- -------- -------
  Food plant                           N                51.163   0.002
                                       P                23.589   0.008
                                      N:P               0.589    0.486
  Grasshopper body                     N                0.648    0.480
                                       P                2.648    0.202
                                      N:P               1.270    0.342
  Frass of grasshopper                 N                0.033    0.868
                                       P                5.691    0.097
                                      N:P               4.166    0.134

Food plant phosphorus content and N:P stoichiometric variation in grasshopper frass {#s3c}
-----------------------------------------------------------------------------------

After performing the corresponding correlation analyses between the indices of N:P stoichiometry of grasshopper and frass with those of food plants, we found that there were no clear relationships (the correlation was not significant at 0.05 level) for most of these variables with only one exception: frass P concentrations were positively and significantly related to food plant P concentrations (*n* = 5, *r* = 0.944, *P* = 0.016) ([Fig. 3e](#pone-0103697-g003){ref-type="fig"}).

![The effects of food plant N:P stoichiometry on N concentrations, P concentrations, N:P ratios in grasshopper body and frass.\
Error bars indicate ±1 SE. In (e), a strong positive relationship was found between frass and food plant P concentrations (ANOVA: *r^2^* = 0.85, *F* = 17.65, *P* = 0.025) and yielded the following equation: y = 3.13 x−1.43.](pone.0103697.g003){#pone-0103697-g003}

Body size of the grasshopper *O. asiaticus* declined linearly with increasing population density {#s3d}
------------------------------------------------------------------------------------------------

In this study, grasshopper body size (expressed as body dry weight) declined significantly in response to increased density ([Fig. 4](#pone-0103697-g004){ref-type="fig"}). As grasshopper density increased from 2 to 32 numbers per square meter, dry weight of individual female decreased ∼22% (from 0.271 to 0.221 g) and showed a significant linear relationship (*n* = 5, *r* = 0.896, *P* = 0.039), while dry weight of individual males declined ∼17% (from 0.089 to 0.074 g) and was also negatively related to density (*n* = 5, *r* = 0.883, *P* = 0.047).

![The relationships between female and male body size (expressed by body dry weight) and experimental density in grasshopper *Oedaleus asiaticus*.\
Error bars indicate ±1 SE. In (a), a negative relationship was found between female body size and grasshopper density (ANOVA: *r^2^* = 0.80, *F* = 12.26, *P* = 0.04) and yielded the following equation: y = −0.0016 x+0.274. In (b), a negative relationship was found between male body size and grasshopper density (ANOVA: *r^2^* = 0.78, *F* = 10.65, *P* = 0.047) and yielded the following equation: y = −0.0004 x+0.085.](pone.0103697.g004){#pone-0103697-g004}

Discussion {#s4}
==========

Distinctive responses of consumer versus producer N and P stoichiometry to grasshopper density {#s4a}
----------------------------------------------------------------------------------------------

The grasshopper *O. asiaticus* in the present study exhibited strong N:P stoichiometric homeostasis under changing food conditions, in this case those generated by experimental manipulations of its own density ([Fig. 2d--f](#pone-0103697-g002){ref-type="fig"}). As a grass-feeder, *O. asiaticus* has plentiful food resource in the Inner Mongolian grassland, which is dominated by several species of plants belonging to Family Poaceae [@pone.0103697-Li1]. Due to its ecological and evolutionary adaptation to environmental variation, this grasshopper species may be able to obtain optimal nutrient balance to meet its nutritional needs and maintain its body N:P within a narrow range, as seen in other insects [@pone.0103697-Bernays1], [@pone.0103697-Kay1]. *O. asiaticus* locusts likely adjust their nutrient intake locally by selecting between different plant species and different parts within one individual plant. *O. asiaticus* is a migratory species [@pone.0103697-Cease1], [@pone.0103697-Jiang1], [@pone.0103697-Cease2] and thus also has the capacity to balance food intake regionally by large-scale movement. However, our study insects were confined within 1-m^3^ experimental cages (*see* [Materials and Methods](#s2){ref-type="sec"}) and thus had few opportunities to select food plants. This limitation in food selection may have been a contributing factor to decreased body size at high density ([Fig. 4](#pone-0103697-g004){ref-type="fig"}).

Contrary to the strong N and P homeostasis in the herbivore, food plant N and P declined with increasing grasshopper density ([Fig.2a,b](#pone-0103697-g002){ref-type="fig"}). However, the ratio of N:P in food plants remained unchanged ([Fig.2c](#pone-0103697-g002){ref-type="fig"}), indicating again a tightly coupled relationship between N and P in the primary producers in this system, consistent with our previous report [@pone.0103697-Zhang1].

Consumer P and N regulation {#s4b}
---------------------------

As mentioned above, when grasshopper density increased in the experiment, food plant N and P concentrations both declined substantially ([Fig.2a, b](#pone-0103697-g002){ref-type="fig"}). Accordingly, P concentrations in grasshopper frass decreased with decreasing plant P content (and increasing grasshopper density) ([Fig.2h](#pone-0103697-g002){ref-type="fig"} and [Fig.3e](#pone-0103697-g003){ref-type="fig"}). However, frass N concentrations changed little, even when grasshoppers were eating plants with different N concentrations ([Fig.2g](#pone-0103697-g002){ref-type="fig"} and [Fig.3d](#pone-0103697-g003){ref-type="fig"}), and the N:P stoichiometry in grasshopper frass was not the functional consequences of that in frass ([Fig. S3](#pone.0103697.s003){ref-type="supplementary-material"}). These results suggest that *O. asiaticus* maintains P homeostasis primarily by changing P absorption and excretion rates, but that other mechanisms may be more important for regulating N homeostasis.

When confined to no-choice diets, grasshoppers and other insects can post-ingestively regulate dietary N [@pone.0103697-Kagata1]--[@pone.0103697-Harrison1]. Zanotto et al (1993) showed that, when confined to synthetic diets containing different concentrations and balances of protein and carbohydrate, fifth instar *Locusta migratoria* nymphs retained a greater portion of ingested N when on low N diets [@pone.0103697-Zanotto1] However, when foraging, herbivores can pre-ingestively regulate protein and carbohydrate balance more strongly than other dietary components, e.g., salts [@pone.0103697-Simpson1], [@pone.0103697-Trumper1]. In our study, *O. asiaticus* may have prioritized regulating N balance pre-ingestively by selecting leaves or leaf parts to meet their macronutrient demands (thus minimizing any change in frass N content) and then post-ingestively regulated P balance. This fits with prior studies on this species that demonstrated changing dietary N and protein content has powerful effects on the growth and viability of *O. asiaticus* [@pone.0103697-Cease2], suggesting that this locust may have a low capacity to post-ingestively regulate for N imbalance and should prioritize behaviorally regulating N balance when foraging.

Density-dependent effects on grasshopper body size, and grassland ecosystem stability {#s4c}
-------------------------------------------------------------------------------------

Many insects respond to high population density by decreasing body size; possible mechanisms include increased stress, competition for resources, and/or hypoxia [@pone.0103697-Peters1]. While our study cannot elucidate the mechanism specifically, our data suggest potential for negative feedbacks of high population density to grasshopper performance through decreased dietary P. Our study measured changes in adult body size, so differences based on population density likely reflect differences in structural material and energy stores and/or egg and sperm production. Here in the present work we did not expand our investigation to these physiological variations in grasshoppers, although some ecosystem level implications may be disentangled.

Increasing grasshopper density has been demonstrated to exert considerable pressure on plant growth and ecosystem productivity [@pone.0103697-Zhang2], [@pone.0103697-Lomer1]. Under a situation of grasshopper population outbreaks, such as the high density treatment in this experiment, decreasing P concentrations and thus increasing N:P ratios in frass ([Fig.2i](#pone-0103697-g002){ref-type="fig"}) may lead to more lower availability of soil P relative to N (skewed consumer-driven nutrient recycling). This high N:P frass may stimulate plant growth, which is generally limited by N availability in Inner Mongolia grassland [@pone.0103697-Xia1]. However, it is important to note that our findings are based on a single field season and involved only adult locusts, limiting our ability to generalize to other conditions. Nevertheless, given that grasshopper and locust plagues can occur in even shorter time scales at given locality [@pone.0103697-Uvarov1]--[@pone.0103697-Chen1] and may result from a swarm of adults settling in a new area (as was simulated by our study), our findings provide useful insights into grasshopper-plant-soil interactions. In addition, global disruption of biogeochemical cycles due to various human activities may induce significant alterations in N and P cycling and their coupling relationships in terrestrial ecosystems [@pone.0103697-Suttle1]--[@pone.0103697-Peuelas1], potentially disrupting grassland plant-herbivores relationships. Thus, a more complete understanding of the stoichiometric impacts of grasshopper-plant interactions on grassland biogeochemical cycling is needed and may be informed by our study.

Conclusions and Implications {#s5}
============================

Our experiment sheds some light on how at least one grasshopper species, *O. asiaticus*, maintains stoichiometric homeostasis in the face of changes in food abundance and elemental composition, highlighting changes in rates of absorption and excretion as an important physiological mechanism in maintaining the stability of body P content. While the ecological significance of frass N production has been reported by several authors [@pone.0103697-Suttle1], we still know little about the effects of P release to soil via frass. To our knowledge, this study is the first to investigate the role of excretion in regulating the P homeostasis in an Acridid herbivore under ecologically realistic conditions. The decreasing amount of frass P input to soil nutrient pools that we document may have significant implications for plant-herbivore dynamics and ecosystem functioning.

Supporting Information {#s6}
======================

###### 

Effects of grasshopper density on frass N:P stoichiometry. In (b), a negative relationship was found between frass P concentration and grasshopper density (ANOVA: *r^2^* = 0.218, *F* = 11.98, *P* = 0.001) and yielded the following equation: y = −0.029 x+2.23. In (c), a positive relationship was found between frass N:P ratio and grasshopper density (ANOVA: *r^2^* = 0.225, *F* = 12.48, *P* = 0.001) and yielded the following equation: y = 0.48 x+25.87.

(DOC)

###### 

Click here for additional data file.

###### 

Comparision of the effects of different grasshopper densities on N:P ratio of food plant, grasshopper body and frass. Error bars indicate ±1 SE. Different letters above different bars mean significant difference at 0.05 level, based on one-way ANOVAs.

(DOC)

###### 

Click here for additional data file.

###### 

Relationship of N:P stoichiometry between grasshoppers body and frass. Error bars indicate ±1 SE. In (a) for N analysis, *r* = 0.378, *P* = 0.530; In (b) for P analysis, *r* = −0.258, *P* = 0.675; In (c) for N:P analysis, *r* = −0.279, *P* = 0.650.

(DOC)

###### 

Click here for additional data file.

We thank Weijun Wu, Yanru Liu, Yuanli Shang and Nanjiao Yu for their assistance during field and lab experiments. We appreciate the Inner Mongolia Grassland Ecosystem Research Station (IMGERS) for providing the experimental sites and other supporting.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: GMZ XGH JJE. Performed the experiments: GMZ ZJZ XMZ QY. Analyzed the data: GMZ ZJZ. Wrote the paper: GMZ ZJZ AJC JJE XGH.
