Nutritional quality modulates trait variability

Background Trait based functional and community ecology is en vogue. Most studies, however, ignore phenotypical diversity by characterizing entire species considering only trait means rather than their variability. Phenotypical variability may arise from genotypical differences or from ecological factors (e.g., nutritionally imbalanced diet), and these causes can usually not be separated in natural populations. We used a single genotype from a parthenogenetic model system (the oribatid mite Archegozetes longisetosus Aoki) to exclude genotypical differences. We investigated patterns of dietary (10 different food treatments) induced trait variation by measuring the response of nine different traits (relating to life history, morphology or exocrine gland chemistry). Results Nutritional quality (approximated by carbon-to-nitrogen ratios) influenced all trait means and their variation. Some traits were more prone to variation than others. Furthermore, the “threshold elemental ratio”- rule of element stoichiometry applied to phenotypic trait variation. Imbalanced food (i.e. food not able to fully meet the nutritional demands of an animal) led to lower trait mean values, but also to a higher variation of traits. Conclusion Imbalanced food led not only to lower trait value averages, but also to higher trait variability. There was a negative relationship between both parameters, indicating a direct link of both, average trait levels and trait variation to nutritional quality. Hence, variation of trait means may be a predictor for general food quality, and further indicate trade-offs in specific traits an animal must deal with while feeding on imbalanced diets. Electronic supplementary material The online version of this article (10.1186/s12983-018-0297-2) contains supplementary material, which is available to authorized users.


Background
The theory of phenotypic plasticity predicts that virtually all biological systems inherently include variation of characters at all organization levels [1][2][3]. This variation can be induced by virtually any external factor, leading to a nearly infinite number of possible phenotypic specificities [4,5]. Phenotypic plasticity can therefore be defined as the degree of a single genotype to express variation in multiple traits leading to different phenotypes in response to environmental gradients [3,6]. Trait based approaches are en vogue in all fields of ecology [7,8], and it has been assumed that trait means of a species can be used to sufficiently characterize communities and affiliated functionality [9,10]. Natural populations, however, consists of phenotypically diverse individuals, possessing variable traits [3,11,12]. Furthermore, trait variability of individuals within populations is essential for evolutionary processes via natural selection, as the latter acts on the individuals' phenotypes influenced by a given ecological setting [1,6,11]. That is why understanding causes and mechanism leading to such trait variation of individuals can give unprecedented resolution to explain evolutionary ecology dynamics underpinning community structure, trait assemblages, and stability (e.g., [11,13,14]).
Theory [15], but also some experimental evidence (e.g., [16,17]) suggest that phenotypic plasticity may arises from trade-offs among multiple traits within a species that are ultimately caused by imbalances in elemental stoichiometry or nutritional components of food. For instance, nutrients influence fitness related traits (e.g., [18][19][20][21]), but also body shapes [22], or intraspecific chemical communication (e.g. [23]) of an animal, and may even alter biological rules [24]. In this context, a resource can be considered optimal if i) it enables an animal to respond with the highest possible trait mean (= balanced diet hypothesis; [25]), but also ii) can sustain phenotypes with low trait variability and thus a stable performance (see [18,20,26]). For instance, two recent meta-analyses suggest that mixed-, rather than single-food diets result in significantly higher average fitness [27] and reduce intraspecific variance in fitness [28], because mixed foods are more likely to represent an optimal, balanced diet close to an animals intake target [25]. Bunning et al. [29] as well as Han & Dingemanse [30] tested these hypotheses in the context of reproductive performance and behavioral stability, respectively, and found that optimal dietary choice indeed led to an improved average expression and a reduced variance of their selected traits. Still, reaction norms of only a relatively limited number of traits towards a broad gradient of macronutrients [20] or elements [31] have been tested for a few species and the relationship of means, variation and elemental balance across multiple traits of different types (e.g. adaptive defense, morphological or reproductive traits) has, to our best knowledge, not been investigated in a single species under genotypic control.
Here, we harness the power of a unique model system the all-female parthenogenetic and mono-genotypic oribatid mite Archegozetes longisetosus ran lineage, which reproduces via automixis with an inverted meiosis resulting in clonal offspring [32,33]. While heritable genetic variation of traits cannot clearly be distinguished from environmentally induced trait variation (phenotypic plasticity) in sexual species, clonal systems provide the possibility to quantify reaction norms of traits and their variation (e.g., [34][35][36][37]). We selected a variety of different individual and demographic life-history, morphological and defensive traits to quantify their response to food with different nutritional quality. We fed ten semi-natural resources, differing in C/N ratios, to the mites (all are descendants from a single mother and share an identical genotype), and measured the mean and variation of traits along this food quality gradient. Thus, we exploit our "phenotypically-tractable" model species (i.e. a model system without genotypic, but only phenotypic variation), to decipher patterns and mechanisms of dietary induced trait plasticity. More specifically we ask how nutritional quality (approximated by C/N ratios) influences the expression of multiple traits regarding i) trait means, ii) trait variation, and iii) differences among traits.
Based on recent meta-analyses [27,28] about the effect of diet mixing, we expect similar patterns for the semi-natural single diets; mites which fed on resources with a balanced C/N-ratio (defined as a ratio equal or close to the mites' stoichiometrical demands) should express a higher average and a lower variance across multiple traits. As no intake target (i.e. species specific optimal elemental/nutritional composition of food) is known for any mite species or soil decomposer it is not possible to predict an optimal C/N intake ratio. However, even the effects of single diets on our phenotypically tractable model species should be strong, as the clonal reproduction eliminates genotypic variation. Hence, we also expect to find a direct negative relationship of the expressed means and variances of traits across resources, which would further support the balanced diet hypothesis [20,25].
As traits we selected life-history, morphological and chemical characters (Fig. 1, Table 1). The life-history traits were selected to quantify the reproductive fitness of one female based on its offspring [35,37,38], while morphological traits were used to describe changes in size and body shape [2]. All chemical traits were related to defensive exocrine opisthonotal oil-glands [39,40]. Oil-gland chemistry is assumed to be highly adaptive because it plays an important role in structuring feeding interactions in soil food webs and enables oribatid mites to live in an "enemy-free space" [41,42], i.e. a conceptual way of living that reduces or eliminate a species' vulnerability against predators [43].
At the start of the experiment we selected young (approximately one week after eclosing) adult individuals from their original culture-plates. We directly used 130 specimens per resource (= 1300 in total) for chemical experiments, while 25 specimens per resource (= 250 in total) were individually redistributed into smaller culture boxes (45 × 40 × 35 mm; grounded with the plaster of Paris mixture) for further experiments.

Chemical experiments and analysis
Oil-gland extractions of 30 specimens from each of the resources (n = 10*30) were accomplished by immersing individuals in 50 μl hexane (GC grade, 98% purity purchased from Merck, Darmstadt, Germany) for 3 min (see Fig. 1). The extraction solvent also contained tetradecane (1 ng/μl; ≥99.8%, analytical standard, purchased from Sigma-Aldrich, Munich, Germany) as an internal standard. Crude extracts were used for gas chromatography-mass spectrometry (GC/MS) analyses and mite specimens were stored at − 20°C for further measurements. Samples were analyzed with a QP 2010ultra GC-MS (Shimadzu, Duisburg, Germany) according to a protocol given elsewhere [44]. The gas chromatograph was equipped with a ZB-5MS fused silica capillary column (30 m × 0.25 mm ID, df = 0.25 μm) from Phenomenex (Aschaffenburg, Germany). The extracted mites were dried at 60°C until weight constancy and the dry weight was determined with a microbalance (Mettler Toledo, XS3DU, 0.1 μg readability and 1 μg repeatability) to calculate the secretion amount per animal/dry weight in [ng/μg], see [44]. Furthermore, we used the chemical data to quantify the relative amounts in [%] of all oil gland compounds based on peak areas. All compounds have been identified previously (for Fig. 1 Conceptual scheme depicting the experimental design of this study. Color code corresponds to the ten resources and is used throughout the study. Circles represent sample replicates; sample sizes are also denoted in the figure details see [44]), and were assigned using their retention indices and diagnostic ions.
Additionally, 100 specimens per resource were chemically disarmed using an established hexane -recovery-hexane (HRH) protocol which does not influence the mites' life-history, but results in the complete depletion of the glands ( [38]; see Fig. 1). Mites were redistributed to culture plates (n = 10 plates per resource, with 10 individuals per plate resulting in 100 samples) and fed with the same resources mentioned above. After 28 days, specimens were individually extracted in hexane to measure the proportion of specimens that had regenerated oil-gland secretions (see [45]; Table 1).

Life history experiments
Individual mites (= mothers, n = 10*25) could lay eggs for ten days and the same food and water was provided ad libitum three times a week (see Fig. 1). Every box was checked daily, and we counted the number of eggs, juvenile instars (i.e. larva, protonymph, deutonymph, tritonymph) and adults for a period of up to 12 weeks. We removed freshly hatched adults from the culture plates to ensure no new egg deposition and froze them for further analysis. The total developmental time [days] for each mother's offspring was calculated as weighted arithmetic mean (developmental time = ∑ [d i *p i ]; where d i is the experiment day and p i is the proportion of new adult specimens on d i ). Replicates where all offspring died In reservoir based chemical defense the amount is a primary factor to predict how often an animal is able to defend itself against predators and competitors [66,67] Composition Relative composition (%) of the defensive secretion of one individual The composition of a defensive chemical blend can determine its effectivity against predators, but can also be a consequence of physiological changes/stress of an individual [66,68,69] Regeneration Percentage of individuals (%) per group (n = 10) which regenerated their defensive secretions over time The regeneration of defensive secretion is essential to be defended against predators at all, but also to understand the costs of secretion regeneration [45,70,71] Life-history The biomass output may also describe the fitness, but also quantifies the ability of a mother to translocate biomass from the resource to her offspring [73,74] Morphology Body mass Dry weight of the initially used females (mothers) Body mass is a universal predictor of many ecological processes (e.g., metabolism, abundance, or predation) [75][76][77] Morphometry Eleven morphometric characters of the initially used females (for details see Additional file 2) The shape of an individual plays an important role in basic physiological processes, but also influences predation (by altered handling approaches by predators) and may be used to predict other characters [22,78,79] before reaching adulthood were not used to calculate the developmental time (n = 72). In addition to the counted data (see Table 1) we dried all individuals at 60°C to determine the total biomass output of individual mothers, expressed as dry weight [mg] per female.

Morphological analysis
The removed mothers from the life-history experiments were used for morphological analysis (see Fig. 1, Table 1), except for damaged individuals that were excluded from the measurements. We overall measured eleven continuous variables of individual specimen (five dorsal and six ventral distances) for the morphometric analysis (see Additional file 2 for details) using a VHX-5000 microscope (Keyence Deutschland GmbH, Neu-Isenburg, Germany) equipped with the VH-Z50L lens. To remove the effect of isometric body size scaling on morphometrical measures we standardized all values of one individual by its notogaster length (see Additional file 2). Afterwards, all specimens (mothers) were dried at 60°C and the individual body mass, expressed as dry weight [μg] was determined as described above.

Analyses of nutritional quality
For C/N analyses dried resource powders (5 ± 1 mg) were weighed into tin capsules. Total organic carbon and nitrogen contents were measured by an elemental analyzer (EA 1108 Elemental Analyser, Carlo Erba, Milan, Italy). Acetanilide (Merck, Darmstadt, Germany) was used as standard. Carbon and nitrogen amounts were calculated based on the standard and the initial dry weight and expressed as C/N ratios ( Table 2).

Data analysis
We analyzed the univariate traits (see Table 1) using Kruskal-Wallis tests and Levene tests to access the overall differences and the variance among resources, respectively. Scatter plots showing all data points and individual posthoc comparisons (Dunn's test [46] and false discovery rate [47] to correct for multiple tests as implemented in "PMCMR") for each trait-resource combination can be found in the Additional files 3 and 4. For multivariate traits [=chemical composition (as Bray-Curtis similarities) and morphometric measures (as Euclidean distances)] we also analyzed the differences and variances among resources using PERMANOVA [48] and PERMDISP [48], respectively as implemented in "vegan". Both multivariate traits were ordinated using discriminant analysis of principal components (DAPC; see [49]) using the package "adgenet". DAPC transforms the original data by principal component analysis (PCA) prior to the discriminant analysis. We retained 6 (for chemical composition) and 7 (for morphometry) PC-axes based on their Eigenvalues (> 1) and the explained variance (total cumulative variance > 95%).
To statistically access the optimality response of a certain trait to the resource quality parameter (i.e., C/N-ratio) we fitted a non-linear regression model (t'~C/N + C/N 2 ; where t' is a given trait value), using the lm-function in R Prior to this statistical analysis we had to adjust our data in several ways: first we transformed the trait data (see results section) to ensure the normality of residuals of every regression and to break potential mean-variance relationships; secondly, we rescaled all the transformed trait values between 0 and 1 using the minimum-maximum normalization (eq. 1) where t is any transformed trait value, t min and t max are the lowest or highest expression level of a given trait, respectively and t' is the min-max normalized trait value. The rescaling was necessary to plot all the optima curves simultaneously and elucidate potential trade-offs among the trait-resource combinations.
To quantify the trait (Table 1) variation for all ten resource treatments and every trait we calculated ln CV (coefficient of variation) according to Nakagawa et al. [50] using eq. (2) where x is the sample mean, σ is the standard deviation and n is the sample size, yielding nine different trait ln CVs for each resource. To summarize the multivariate traits (chemical composition and morphometry) we first calculated the ln CVs for each single component and subsequently averaged the single ln CVs to obtain the mean ln CV for both traits. We used a Kruskal-Wallis Table 2 Carbon to nitrogen ratios (mean ± standard deviation; n = 3 replications) of the ten food resources sorted by C/N-ratio, beginning with the lowest (= highest nitrogen content)  test to check whether the nine trait ln CVs across all resources are prone to more/less variation compared to others and if trait types (chemical, life history, morphology) possess different variability using the mean across resource ln CVs. We further analyzed the mean trait variation across traits to test whether some resources produce more variable traits than others and if the variation of trait variation differs across resources using a Kruskal-Wallis and Levene test, respectively. Afterwards, we used these values to calculate the mean ln CV for each resource (see last row of Table 3), to eventually test for optimality of C/N-ratios using a non-linear regression model (using the package "Rcmdr"). Additionally, we used the same analysis to test the influence of the C/ N-ratios on the ln CVs of the nine individual traits. Finally, we accessed the relationship of mean and variance of the different traits across all resources by using the mean normalized trait value t' mean and the ln CV of the univariate traits (because multivariate traits do not have a mean t'). Because the different traits within one resource treatment are not independent from each other we used a linear mixed effect model (LMM) with trait variation as dependent, trait mean as independent and resource treatment as random variable using the R packages "lme4" and "car". The model was checked for normality of residuals and variance homogeneity, and no further transformation was necessary.

Nutritional quality
Some food resources were characterized by an extreme compositione.g. blood meal had a very high C/N ratio, which corresponds to a high nitrogen content (13.3 ± 0.1% N; mean ± SD), while pollen was nitrogen poor, but comparatively carbon rich (3.6 ± 0.1% N; mean ± SD). Yet, most of the remaining resources had a quite equilibrated stoichiometry (Table 2). Thus, the variability of the C/N ratios across the food was high (CV = 47%; 7.5 ± 3.5; mean ± SD), and the spanned from 3.3 (blood meal) to 13.6 (pollen). While we focused our analyses on C/N ratios only, the foods also differed in other nutritional dimensions (i.e. elements and macronutrients; see Additional file 1), but their influence is beyond this study.

Discussion
Dietary induced intraspecific reaction norms and variation are a mostly neglected topic in evolutionary ecology [2,34]. Over 100 years ago, Woltereck [34] started to discuss underlying mechanisms, but also noted that it is hardly possible to quantify the range of a species' variability caused by nutritional quality, because even in pure linages of facultative parthenogenetic species, like Daphnia, a full control of the genotype is not possible. He also suggested that an obligate asexual species would be a perfect model to study nutritional reaction norms of traits and their plasticity. That is why parthenogenetic generalists (like some oribatid mites), which can tolerate a wide range of environmental conditions (i.e. possess a general purpose genotype; [57]), are the models that Woltereck [34] proposed. Correspondingly, we used the single genotype of a "phenotypically-tractable" model lineage A. longisetosus ran to unravel the relationship of food quality and variation of multiple traits (i.e. phenotypic plasticity).

Nutritional effects on traits and variation
Most studies on a single species so far focused on the relationship of nutrients to one or a distinct set of traitsoften times related to fitness (e.g., [18,20,58,59]). We think that our multiple traits approachalso including presumable neutral charactersbetter allows to disentangle patterns and mechanisms of nutrient influence on the mean and variance as well as the inherent plasticity of traits. In our experiment, all traits responded to diet and showed quadratic reaction norms (= optima) of different strength (Fig. 2). Physiological theory formally conceptualized this quadratic (or concave) response of fitness relevant performance traits as "threshold elemental ratio" [for element stoichiometry; 26] or as "Bertrand's rule" for the concentration of essential  [60] and also macronutrients [18]. While the first concept is a stoichiometric approach based on the proportional relationship of elements, the latter uses absolute concentrations or the amount of nutrient intake to explain reaction norms of individuals [18,20,26,31]. Yet, both concepts are unified by a similar prediction: if a food is too unbalanced (containing too low levels of a certain nutrient or element) to reach the intake target, the benefits gained from increasing the amount or ratio of this nutrient or element also increases until they reach an equilibrium. More nutrients/elements beyond this threshold (the intake target) are associated with increasing costs for the regulatory mechanisms resulting in physiological disadvantage higher than the original benefit [20,31]. Our results for C/N show that a simple threshold elemental ratio (presumably close to our model species intake target) applies to multiple traits and, even more important, also to their variability and thus to phenotypic plasticity. This is because the variability of each trait across resources was heteroscedastic, indicating that food quality not only changed the mean (e.g., [61,62]), but also the variation of a trait [25,[28][29][30]. These findings could help to propose a general eco-physiological mechanism causing dietary related intraspecific trait variation derived for our model system: the low performance and high plasticity of animals a b  (Table 1) across all ten resources (a) and the overall trait variation [ln CV] of each resource across all traits (b). For (a), individual traits are grouped into their respective category (chemistry, life-history or morphology). Circles represent means, lines denote the standard error, colors in (b) correspond to Fig. 1 feeding on resources from the "edges" of an elemental or nutritional gradient ingest imbalanced food with stoichiometric shortcomings which causes stress related costs to deal with the deficits as well as surpluses of elements or nutrients [18,26]. While at low N-content (high C/N ratio) trait performance was limited by a short supply of protein, mites shifted to a C-limited trait performance at high level of N (low C/N). Consequently, high variability of all traits (i.e. high phenotypic plasticity) occurred if either protein (N) or energy (C) limited the formation and performance of traits. More general (see Fig. 6): if there is a considerable trade-off between consumed surplus and deficit nutrients in one resource, the variation across multiple traits within a phenotype is high, while the trait mean is low. Reciprocally, trait variation is low, and the trait mean is high at nutritional optima where no essential nutrient is in short supply and a respective genotype can almost bear its full potential. This connection of trait means and variances across resources (expressed as a direct negative relationship, Fig. 6) is further evidence for the wide applicability of the balanced diet hypothesis [25]. The concave responses of all traits and their variability to overall food quality (C/N ratio) in our experiments not only indicated that the threshold-elemental-ratio-rule applied for a wide range of traits and their variability, but further suggest costs [20] to maintain a high mean and low variance of a trait. This means that mites feeding on "edge" resources must deal with high costs during allocating resources compared to mite consuming the "optimum" food. Despite these costs, however, variation at the "edges" may still be beneficial, because it enables at least a small number of individuals to survive unfavored conditions. For instance, a recent synthesis by Forsman and Wennersten [63] found that variation seems to be more important under stressful circumstances when animals are forced to exist under suboptimal conditions -like stoichiometrically imbalanced food not meeting the intake target [26] -and may enable the survival of a population [64].

Inherent variability of traits
Generally, there are no multicellular organisms without a certain plasticity, because intra-individual trade-offs as reactions towards environmental conditions like temperature, salinity or resource availability, but also biotic factors like predation will lead to variation of traits [15,63,65]. Besides selection, it is still poorly understood why some traits are more prone to variation than others, i.e. bear higher plasticity. Our data indicates that some traits tend to be more variable than other, but this was not related to certain "trait types" (in our case "chemistry", "life history" and "morphology"). This different variability potential allows to derive different hypothesis: a certain trait may react with higher variability to an environmental gradient, because it faces more trade-offs along this gradient than other characters; or the formation/development of a trait may be more "complex" and thus demands a stronger segregation of energy leading to more trade-offs. Consequently, lower trait variation may indicate less trade-offs (or selection) along a certain gradient. Also, a lower variability may be a signal for an inherently lower plasticity of a trait because it is less controlled by the phenotypes' response or selection. For instance, in our experiments, the body size of the mites changed considerably along the C/N gradient, yet the overall body shape (morphometric measurements) only showed a weak response and low variability. This may indicate, thatbesides selection or genetical conservationthe overall body size could be influenced by the phenotypic response to altered nutrients. Yet, the proportions of the body shape isometrically scaled with this phenotypic change, leading to low overall variability.

Conclusions
Overall, we have shown that the threshold elemental ratio rule [18,26,60], applies to a wide range of traits and also to dietary caused intraspecific variation of multiple traits (= phenotypic plasticity) in a model system excluding genotypic variation. Thus, there seems not only to be a nutritional optimum for a traits' mean, but also its variation, which further strengthens the balanced diet hypothesis [25]. Highly imbalanced food results in lower trait means and higher trait variability (Figs. 2, 5, 6), potentially caused by more/stronger trade-offs across resources arising from physiological stress. Additionally, some traits seem to have higher inherent variation than others, irrespective of "trait type". Based on our empirical evidence we propose the following hypotheses, leading to further experimental evaluation; i) the mean trait variation calculated based on all trait variation within a resource could be used as an indicator for the food quality and ii) given that there is a direct link of both mean and variation of traits (Fig. 6), the latter might also be indicative for trade-offs an animal faces when feeding on imbalanced diets [20]. Finally, using a single parthenogenetic lineage excludes genotypic plasticity and creates a phenotypically tractable model species. This is a first step to establish an experimental platform to ultimately disentangle how variation (or the capacity to express plasticity) affects the fitness of individuals [65].  Fig. 1