Linking organismal growth, coping styles, stress reactivity, and metabolism via responses against a selective serotonin reuptake inhibitor in an insect

Linking organismal growth, coping styles, stress reactivity, and metabolism via responses against... www.nature.com/scientificreports OPEN Linking organismal growth, coping styles, stress reactivity, and metabolism via responses against Received: 9 October 2017 a selective serotonin reuptake Accepted: 15 May 2018 Published: xx xx xxxx inhibitor in an insect 1,2,3 4,5 5 6 5 Indrikis Krams , Giedrius Trakimas , Sanita Kecko , Didzis Elferts , Ronalds Krams , 7,8 9 10 10 11 Severi Luoto , Markus J. Rantala , Marika Mänd , Aare Kuusik , Jukka Kekäläinen , 12 11 5,10 Priit Jõers , Raine Kortet & Tatjana Krama Evidence suggests that brain serotonin (5-HT) is one of the central mediators of different types of animal personality. We tested this assumption in field crickets Gryllus integer using a selective serotonin reuptake inhibitor (SSRI). Crickets were selected for slow and rapid development and tested for their coping styles under non-stressful conditions (time spent exploring a novel object). Resting metabolic rate, maximum metabolic rate and latency to resume activity were measured under stressful conditions (stress reactivity). Measurements were taken (i) before and (ii) during the SSRI treatment. Before the SSRI treatment, a strong negative correlation was observed between coping style and stress reactivity, which suggests the existence of a behavioral syndrome. After the SSRI treatment, the syndrome was no longer evident. The results of this study show that 5-HT may be involved in regulating behavior not only along a stress reactivity gradient but also along a coping styles axis. The relationship between personality and the strength and direction of 5-HT treatment on observed behaviors indicates trait-like individual differences in 5-HT signaling. Overall, these findings do not support recent ideas arising from the pace-of-life syndrome (POLS) hypothesis, which predict higher exploration and metabolic rates in rapidly developing bold animals. Phenotypic differences in trait expression are a prerequisite for evolution through natural selection. Although most of the traits of an organism vary around a mean value, research on human and animal behavior has revealed between-individual clusters in behavioral responses, oen a ft ttributed to different personality axes that are called personality types, behavioral predispositions, temperaments or coping styles. Some of these behavioral responses are 1,2 3–7 not heritable , while others demonstrate a negligible or a considerable degree of heritability . Animal person- 5,6,8–13 ality is tightly linked to individual life histories, survival and fitness . e Th pace-of-life syndrome (POLS) hypothesis, rooted in the classic concept of r- and K-selection, is particu- 14–16 larly relevant in this context , as it suggests that life history traits such as age at maturity and growth rate are likely coupled with individually consistent behaviors and immune responses. For example, shy individuals oen ft 1 2 Institute of Ecology and Earth Sciences, University of Tartu, Tartu, Estonia. Department of Zoology and Animal Ecology, Faculty of Biology, University of Latvia, Rīga, Latvia. University of Tennessee, Department of Psychology, 4 5 Knoxville, USA. Institute of Biosciences, Vilnius University, Vilnius, Lithuania. Department of Biotechnology, Institute of Life Sciences and Technology, Daugavpils University, Daugavpils, Latvia. Department of Botany and Ecology, Faculty of Biology, University of Latvia, Rīga, Latvia. English, Drama and Writing Studies, University of Auckland, Auckland, New Zealand. School of Psychology, University of Auckland, Auckland, New Zealand. 9 10 Department of Biology & Turku Brain and Mind Centre, University of Turku, Turku, Finland. Department of Plant Protection, Institute of Agricultural and Environmental Sciences, Estonian University of Life Science, Tartu, Estonia. 11 12 Department of Environmental and Biological Sciences, University of Eastern Finland, Joensuu, Finland. Insttute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia. Correspondence and requests for materials should be addressed to I.K. (email: indrikis.krams@ut.ee) SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 1 www.nature.com/scientificreports/ grow slower, reach larger body size, mature later, invest more in immune system and have longer hiding time in an anti-predator context under familiar conditions than bold individuals . Aggressiveness of bold individuals hypothetically facilitates acquiring and monopolizing resources and may be linked to a fast growth rate and more 5,17 intense reproductive effort early in life . Personality research is a complex field that adjoins behavioral ecology with psychology and molecular biology. It is therefore important to define key terms before testable hypotheses are proposed. The term behavioral syn - 18,19 drome indicates that behavioral trait characteristics involve suites of correlated behaviors that are consistent over time and across contexts. The term coping style is more commonly used to characterize individual trait dif- ferences in biomedicine, and this concept has been useful in explaining behavioral responses and their underlying physiological mechanisms. The coping style concept implies that animals may react to a stimulus with alternative response patterns. According to the two-tier model of animal personality proposed by Steimer et al. , coping styles reflect the quality of the response to a stressor, while stress reactivity reflects the quantity of the response expressed. In this model, bold individuals are those with a strong tendency to act (proactivity) combined with low stress reactivity. Docile individuals can be characterized as a combination of reactive coping associated with low stress reactivity. The combination between high stress reactivity and a high tendency for proactive coping is labeled as panicky. Shy individuals exhibit reactive coping style and high stress reactivity. Although those personality characteristics are common, researchers rarely use all of the above terms because that would require the use of several tests for each individual being studied. For example, stress reactivity needs to be measured by obtaining levels of stress hormones/neuromediator concentrations, or by observing the organ- isms under stressful conditions. In contrast, reactive and proactive coping styles differ in behavioral flexibility, where the proactive individuals act on the basis of previous experiences and develop routines in stable environ- 21,22 ments . The reactively coping animal tends to rely more on a detailed appraisal of their current environment because they react to immediate environmental stimuli and tend to explore any changes in their environment. 23,24 The metabolic rate of an organism varies with its body size and body temperature . The POLS hypothesis proposes that rapidly developing, oen b ft old individuals have life histories involving high daily energy expendi- tures and higher basal or resting metabolic rates (RMR) because they are forced to process more food and inevi- tably excrete more waste products at a faster rate compared to slowly developing individuals of same size. Rapidly developing individuals may, however, remain smaller than slow developers . Thus, personality also has the potential to ae ff ct the metabolic rate of an organism. Another important but oen n ft eglected source of variability in RMR arises because individuals consistently differ in their stress responses to conditions in which metabolism 25,26 measurements are performed . For example, because reactive individuals (shy, docile) oen b ft ecome immobile 12,13 27 under stressful conditions , their metabolic rates may be mistakenly classified as RMR . Recent research shows that while the latency to resume body movements aer h ft andling in a familiar environment is longer for slowly developing crickets, they resume body movements sooner than rapidly developing crickets when in a more stress- ful, unfamiliar environment . Slowly developing crickets were found to have significantly higher RMR com- pared to rapidly developing crickets. It is important to note that according to established experimental procedure, measurements of respiration always require placing an organism in a novel environment . This has the potential to induce significant stress, which may force less stress-resistant animals to stay active during respirometry to explore or escape their respirometry chambers. This demonstrates the importance of n fi ding the proper approach to measurements of stress reactivity and metabolism, as well as making a correct interpretation of the results obtained to explain the relationships between animal personality and underlying physiological mechanisms . Behavioral responses of shy, oen a ft nxious individuals are associated with high concentrations of stress hor - mones , implying that their slowness and caution are caused by physiological stress. In psychology and physiol- ogy, shyness is usually defined as the feeling of anxiety and awkwardness when an individual is under conditions of stress. Stronger forms of shyness are usually referred to as anxiety and elevated neuroticism. In extreme cases, it might be related to a depression-like psychological condition, which includes the experience of fear to the extent 31 20,29 of inducing panic . This calls for reliable tests to discriminate between shy and panicky personality types . In the central nervous system, serotonin (5-hydroxytryptamine: 5-HT), a monoamine neurotransmitter , has 33–36 several important functions, including the regulation of aggression, mood, sleep, dispersal and metabolic rate . Modulation of serotonin at synapses of the brain is thought to be the key mechanism of several classes of pharma- cological antidepressants. Experimental evidence has revealed links between low levels of 5-HT and the develop- 37,38 ment of anxiety and stress-prone behavioral phenotypes . However, an association between low concentrations of 5-HT and proactive coping style suggests that 5-HT may be involved in the regulation of behavior not only 20,29,40 34 along the stress reactivity gradient but also along the coping styles axis (but see ). These findings indicate that 5-HT potentially ae ff cts personality types in a more complex way, which should be taken into account when testing the influence of pharmacological compounds on behavioral responses under stress. 41–43 Field crickets are a suitable and often-used model group for animal personality research . In this study, we compared coping style differences between slowly and rapidly developing western stutter-trilling crickets (Gryllus integer) under non-stressful and stressful conditions. Prior to the study, the crickets had been selected for slow and rapid developmental time over five generations. We also measured their RMR, maximum metabolic rate (MMR) and the latency to resume activity under stressful conditions in the insect chamber of respirometer. This was done twice: before a selective serotonin reuptake inhibitor (SSRI) was added to the food and once again during the anti-depressant treatment. We expected greater effects of anti-depressant treatment on behavioral responses of proactive individuals (bold, panicky) − those paying little attention to minor changes in their environment. Although some earlier evidence suggested that central injection of 5-HT or its precursor 5-hydroxytryptophan (5-HTP) may cause 44,45 46 increased metabolic rates (but see ), other research has shown that serotonin deficiency significantly increases 47–49 amino acid, energy, purine, lipid and gut microflora metabolisms and rates of oxidative stress . Therefore, we predicted an SSRI-induced decrease in the metabolic rate of crickets with high stress reactivity (shy, panicky), SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 2 www.nature.com/scientificreports/ Figure 1. Mean (±95%CI) body mass (a), coping styles (b), time to resume activity (c), RMR (d) and MMR (e) before SSRI (▫, ◦) and aer SS ft RI treatment (▪, ⦁) in crickets of rapid (n = 41), slow (n = 37) and control (n = 28) developmental lines. ◦, ⦁ indicate females, ▫, ▪ indicate males. 50,51 especially in panicky individuals , because 5-HT is likely to influence proactive individuals (panicky, bold) 36,52 more than reactive ones (shy, docile) . Results Body mass. Our results demonstrate that SSRI can have relatively small but significant effects on body mass of crickets (rmANOVA, F = 22.27, P < 0.0001, partial η = 0.18, Figs 1a and 2a). The average body mass (1, 100) of crickets decreased from 0.548 ± 0.004 g to 0.539 ± 0.004 g (unweighted means ± SE) during SSRI treatment. e Th body mass of crickets largely depended on the developmental line (rmANOVA, F = 356.9, P < 0.0001, (2, 100) partial η = 0.85). The effect of sex, although much smaller than that of developmental line, was significant (rmANOVA, F = 5.41, P = 0.02, partial η = 0.05). Before the SSRI treatment, slowly developing crickets (1, 100) were heavier (0.701 ± 0.006 g (weighted means ± SE) than rapidly developing (0.466 ± 0.007 g) and control crick- ets (0.466 ± 0.011 g), while females (0.558 ± 0.016 g) were heavier than males (0.537 ± 0.017 g). Although body mass decreased in males and females of all developmental lines during the SSRI treatment, slowly developing crickets remained the heaviest among the lines, while females were heavier than males. All possible interac- tions between factors were non-significant. No aging and testing apparatus effects were evident (rmANOVA, F (1, = 1.0, P = 0.323, partial η = 0.02). Body mass was repeatable in control (R = 0.69) and highly repeatable in 54) slowly (R = 0.93) and rapidly developing crickets (R = 0.94; Table 1). Coping style. We found significant main effects of developmental lines (rmANOVA, F = 70.45, (2, 100) 2 2 P < 0.0001, partial η = 0.56) and SSRI treatment (rmANOVA, F = 10.19, P = 0.002, partial η = 0.09) to (1, 100) coping styles of the crickets, while sex apparently had no effect on it (rmANOVA, F = 0.1, P = 0.76, partial (1, 100) η = 0.001) (Figs 1b and 2b). However, there was significant interaction between developmental line and antide- pressant treatment (F = 42.03, P < 0.0001, partial η = 0.46), suggesting that simultaneous ee ff cts of those (2, 100) factors on the coping style were not additive. Before the SSRI treatment, slowly developing crickets explored the novel object longer than control (Bonferroni test: P = 0.0002) and rapidly developing crickets (Bonferroni test: P < 0.0001); control crickets explored the object longer than rapidly developing individuals (Bonferroni test: P < 0.0001). During the SSRI treatment, rapidly developing crickets decreased their proactivity by increas- ing exploration time of the novel object (Bonferroni test: P = 0.0008). Slowly developing crickets decreased the proportion of time devoted to exploration (Bonferroni test: P < 0.0001) (Fig. 3), while the control individuals explored the novel object equally long (Bonferroni test: P > 0.05). Aer t ft he SSRI treatment, coping styles of all developmental lines became statistically similar (Bonferroni test: P > 0.05). No aging nor testing apparatus effects were evident (rmANOVA, F = 0.93, P = 0.339, partial η = 0.02). Coping styles were repeatable in the control (1, 54) (R = 0.91) and slowly developing crickets (R = 0.67) but not in the rapidly developing crickets (R = 0.16; Table 1). However, combined repeatability was high (R = 0.94) when calculated for all selection lines together. SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 3 www.nature.com/scientificreports/ Figure 2. Relative change in mean body mass (a), coping styles (b), time to resume activity (c), RMR (d) and MMR (e) before SSRI and aer SS ft RI treatment in crickets of rapid, slow and control developmental lines. Significance based on Tukey’s test: n.s. P > 0.05, *P < 0.05, **P < 0.01, ***P < 0.001. Control Slow Rapid Variables R 95%CI limits F R 95%CI limits F R 95%CI limits F 19,19 19,19 19,19 *** **** **** Body mass 0.69 0.38; 0.86 5.75 0.93 0.82; 0.97 25.1 0.94 0.85; 0.98 29.8 a **** *** ns Coping styles 0.91 0.80; 0.97 21.4 0.67 0.33; 0.85 4.95 0.16 −0.15; 0.50 1.58 **** **** **** Stress reactivity 0.99 0.98; 1.00 226.3 0.92 0.82; 0.97 25.5 0.97 0.93; 0.99 66.4 **** **** **** RMR 0.96 0.90; 0.98 51.1 0.91 0.80; 0.97 22.9 0.79 0.51; 0.92 10.2 **** **** **** MMR 0.94 0.87; 0.98 33.7 0.94 0.86; 0.98 33.5 0.94 0.86; 0.98 32.4 Table 1. Repeatabilities of traits in crickets of rapid (n = 41), slow (n = 37) and control (n = 28) developmental *** **** ns lines estimated using two-way random effects model for single measures. P < 0.001; P < 0.0001; P > 0.05; combined all selection lines: R = 0.94 (0.91; 0.97), F = 33.8, P < 0.001. 59,59 Stress reactivity. We found significant main effects of antidepressant treatment on time to resume activ- ity in the insect chamber (rmANOVA, F = 33.83, P < 0.0001, partial η = 0.25) (Figs 1c and 2c). Time to (1, 100) resume activity increased from 221 ± 12 s (unweighted means ± SE) before the SSRI treatment to 284 ± 11 s aer ft the treatment. Stress reactivity of the crickets depended on developmental line (rmANOVA, F = 61.73, (2, 100) P < 0.0001, partial η = 0.55), while sex had no effect on stress reactivity (rmANOVA, F = 0.12, P = 0.73, (1, 100) partial η = 0.001) either before or after the SSRI treatment. After adding antidepressants to the food, all selec- tion lines increased their time to resume active movements in the insect chamber. No significant interactions were observed. No aging and testing apparatus effects were found (rmANOVA, F = 0.41, P = 0.523, partial (1, 54) η = 0.01). Stress reactivity was highly repeatable in all selection lines (range R = 0.92–0.99; Table 1). Resting metabolic rate. A significant main effect of developmental line on RMR was observed (rmANOVA, F = 34.63, P < 0.0001, partial η = 0.41) (Figs 1d and 2d). However, there was a significant (2, 100) interaction between developmental line and antidepressant treatment (rmANOVA, F = 25.52, P < 0.0001, (2, 100) partial η = 0.38). Before the SSRI treatment, the RMR of slowly developing crickets was higher than that of control (Bonferroni test: P < 0.0001) and rapidly developing crickets (Bonferroni test: P < 0.0001). The RMR of rapidly developing crickets was lower than that of control crickets (Bonferroni test: P = 0.008). After the treat- ment, RMR significantly increased in rapidly developing crickets (P = 0.001), significantly decreased in slowly developing crickets (Bonferroni test: P < 0.0001) (Fig. 4), but remained the same in controls (Bonferroni test: P = 1). Aer t ft he antidepressant treatment, the RMR of slowly developing crickets was statistically similar to the RMR in the control line before the SSRI treatment (Bonferroni test: P = 0.29). As the mean RMR data show, all the SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 4 www.nature.com/scientificreports/ Figure 3. Change in coping styles aer SS ft RI treatment in individual crickets of rapid (──, n = 41) and slow (……, n = 37) developmental lines (a) and control (n = 28) line (b). Figure 4. Change in RMR aer SS ft RI treatment in individual crickets of rapid (──, n = 41) and slow (……, n = 37) developmental lines (a) and control (n = 28) line (b). lines became more similar aer t ft he SSRI treatment. However, some of the developmental lines remained signifi- cantly different from each other (i.e., rapidly developing crickets vs slowly developing crickets: P = 0.004; slowly developing crickets vs control group: P = 0.02). RMR was repeatable in rapidly developing (R = 0.79), and highly repeatable in slowly developing (R = 0.91) and control crickets (R = 0.96; Table 1). Maximum metabolic rate. A small significant main effect of developmental line on MMR (rmANOVA, F = 4.03, P = 0.02, par tial η = 0.08) was observed, while sex had no effect on MMR (rmANOVA, (2, 100) F = 2.26, P = 0.07, partial η = 0.03) (Figs 1e,2e). However, there was a significant interaction between (1, 100) developmental line and SSRI treatment (rmANOVA, F = 64.38, P < 0.0001, partial η = 0.56). Before the (2, 100) SSRI treatment, slowly developing crickets had the highest rates of maximum metabolism, followed by control and rapidly developing crickets. Significant differences were found between rapidly and slowly developing lines (Bonferroni test: P < 0.0001) and between rapidly developing and control crickets (Bonferroni test: P = 0.006), while slowly developing crickets and the control group did not differ (Bonferroni test: P = 1). During the SSRI treatment, MMR significantly increased in rapidly developing crickets (Bonferroni test: P < 0.0001) and decreased in slowly developing individuals (Bonferroni test: P < 0.0001), while it did not change significantly in the control group (Bonferroni test: P = 1). Aer t ft he SSRI treatment, the MMR of all developmental lines was sta- tistically similar (Bonferroni test: P = 1). MMR was highly repeatable in all selection lines (all R = 0.94; Table 1). Associations among traits. Before the SSRI treatment, the crickets’ body mass was positively correlated with RMR (r = 0.63, P < 0.0001), MMR (r = 0.34, P < 0.01) and coping styles (r = 0.69, P < 0.0001), while being SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 5 www.nature.com/scientificreports/ Traits Body mass Coping styles Stress reactivity RMR MMR **** *** −0.14 −0.63 0.40 −0.08 Body mass (−0.34; 0.06) (−0.73; −0.53) (0.22; 0.55) (−0.29; 0.13) **** ** 0.69 0.00 0.33 0.09 Coping styles (0.59; 0.78) (−0.17; 0.17) (0.8; 0.53) (−0.09; 0.26) −0.32 **** **** −0.66 −0.60 −0.03 Stress reactivity (−0.48; (−0.74; −0.57) (−0.73; −0.47) (−0.21; 0.16) −0.15) **** **** **** 0.63 0.69 −0.74 0.09 RMR (0.50; 0.74) (0.56; 0.79) (−0.80; −0.67) (−0.10; 0.27) ** ** *** *** 0.34 0.37 −0.40 0.40 MMR (0.13; 0.52) (0.20; 0.53) (−0.55; −0.24) (0.21; 0.56) Table 2. Pearson correlation coefficients between traits in crickets before SSRI treatment (below diagonal) and aer SS ft RI treatment (above diagonal); lower and upper 95% CI values are presented in parentheses; all P-values * ** **** are Bonferroni corrected. P < 0.05; P < 0.01; ***P < 0.001; P < 0.0001; (all Bonferroni corrected). Figure 5. e t Th wo-tier model of cricket personality before the SSRI treatment (a), and aer t ft he SSRI treatment (b). ○ = rapid developmental line (n = 41), ◊ = slow developmental line (n = 37) and □ = control developmental line (n = 28). negatively associated with stress reactivity (i.e. time to resume activity) (r = −0.66, P < 0.0001, Table 2). Aer t ft he SSRI treatment, the relationships between the crickets’ body mass, coping styles and MMR disappeared, while the associations were still apparent between body mass and stress reactivity (r = −0.63, P < 0.0001) and body mass and RMR (r = 0.40, P < 0.001). RMR before the SSRI treatment was positively related with MMR (r = 0.40, P < 0.001), coping styles (r = 0.69, P < 0.0001) and negatively with stress reactivity (r = −0.74, P < 0.0001). Aer ft the SSRI treatment, the relationships between RMR, coping style and stress reactivity were weaker (r = 0.33, P < 0.01; r = −0.32, P < 0.05, respectively, Table 2). Before the SSRI treatment, MMR was also associated with coping styles (r = 0.37, P < 0.01) and stress reactivity (r = −0.40, P < 0.001), while aer t ft he treatment, MMR was not correlated with behaviors or RMR (Table 2). Prior to the SSRI treatment, a strong negative correlation was observed between the coping style and stress reactivity (r = −0.60, P < 0.0001), suggesting a behavioral syndrome (Fig. 5a). Aer t ft he treatment with antidepressants, however, the syndrome disappeared (Table  2, Fig. 5b). Discussion Evidence shows that slow development is associated with upregulation of stress-related genes . In humans, a predisposition toward experiencing negative emotions and low inhibitory control are linked to sickness and many 54–56 psychiatric conditions . It is likely that slower development in crickets confers higher levels of anxiety/neu- roticism along the stress reactivity axis since the SSRI treatment increased the time to resume movements (i.e. decreased anxiety) under the stressful conditions of the insect chamber. The relaxed stress reactivity under the SSRI treatment was associated with significant decreases in the RMR and MMR of slowly developing crickets (Fig. 6). The SSRI treatment presumably increased the concentration of 5-HT in the brain, while 5-HT could have lowered body temperature by decreasing the metabolic rate (Fig. 6). 5-HT activates glycolysis, forming lactic acid. 57,58 Excess lactic acid tends to decrease efficient energy production by interfering with mitochondrial respiration . On the other hand, the SSRI treatment did not ae ff ct stress reactivity of rapidly developing crickets. Moreover, the treatment that brings higher concentration of brain 5-HT increased RMR, MMR and exploration time of rapidly developing, oen b ft old crickets, while it decreased exploration time of slowly developing, oen s ft hy crickets. High concentration of brain 5-HT may cause a serotonin syndrome that oen in ft cludes high body temperature, indicat- ing a higher metabolic rate . Richmond et al. found that aquatic invertebrates are highly sensitive to pharmaceutical pollution in their environment: dipteran midges emerged significantly earlier and the total biomass of emerged adult insects was SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 6 www.nature.com/scientificreports/ Figure 6. A flow chart showing trait associations and the effects of SSRI treatment in slowly and rapidly developing crickets. Slow development is associated with higher (resting and maximum) metabolic rates (MR), stress reactivity (SR), exploration (Expl) and shy personality. Rapid development is associated with lower (resting and maximum) metabolic rates, stress reactivity, exploration and bold personality. By presumably raising 5-HT levels in the brain, SSRI decreases key trait values (MR, SR, Expl) in slowly developing crickets and raises these values in rapidly developing crickets. This results in the “normal” (intermediate) personality phenotype in slowly and rapidly developing crickets that is typical of the control line. greater in fluoxetine- and citalopram-exposed streams when compared with controls. This suggests that SSRI presence in the environment induces selection for faster development and more intense metabolism. This could mean that (i) rapidly and slowly developing crickets have different 5-HT tolerance, 5-HT sensitivity in brain receptors or different 5-HT concentration in the brain; (ii) 5-HT ae ff cts coping styles rather than responses along the stress reactivity axis and/or (iii) serotonergic metabolism is linked to growth. Under natural conditions, pro- active (active) coping may have higher survival value especially under high density of conspecifics/stronger com- petition than a reactive (passive) coping style . On an evolutionary scale, however, survival values for different behavioral types generally depend on environmental variability and predictability . Prior research has revealed differences in sensitivity to 5-HT levels in D. melanogaster . Ns3 mutant flies (which reach less than 60% of normal size and have fewer and smaller cells, but exhibit normal body proportions) and wild-type flies were treated with 5-HTP (5-hydroxy-L-tryptamine), the immediate precursor of 5-HT. This treatment resulted in a developmental delay in wild-type animals, reminiscent of the ns3 mutant phenotype, while ns3 mutants were more sensitive to 5-HTP than their wild-type counterparts. This suggests that ns3 mutants have already increased serotonergic signaling, as is likely in rapidly developing crickets. Using selective 5-HT1A and 5-HT1B agonists, de Boer and Koolhaas found a manyfold die ff rence between high- and low-aggression male rats. The high-aggression rats had a far more sensitive 5-HT1A auto-receptor and 5-HT1B receptor-mediated inhibition of the serotonergic neurons than low-aggression rats. However, the mechanisms that underlie these effects remain elusive. One reason is that the spatial distribution of serotonergic receptors remains poorly located and characterized in both the vertebrate and insect brain. As pointed out by de Boer and Koolhaas , it is crucial to know whether the receptors are located postsynaptically or presynaptically. At the presynaptic level, 5-HT agonists reduce 5-HT signaling, whereas at the level of the postsynaptic receptors, they mimic the effects of enhanced 5-HT signaling. In future research on 5-HT1A and 5-HT1B, agonists need to be tested in slowly and rapidly developing crickets to reveal any differences in the sensitivity of their 5-HT receptors and its relation to the development of personality along the stress reactivity and coping styles axes. Interestingly, the SSRI treatment decreased stress reactivity (i.e. increased the time to resume activity) only in slowly developing crickets while it ae ff cted exploration time and coping styles of both slowly and rapidly devel - oping individuals, but in the opposite way. Koolhaas et al. suggested a mechanism to explain the opposing effect of 5-HT on coping styles of individuals having to different personalities. Evidence shows that activation of the 5-HT1A auto-receptor leads (i) to a reduction of the preferred behavioral response in a forced swim test in rats and (ii) an increase in an alternative behavioral response rather than to a general reduction in anxiety . It is known that 5-HT1A receptor agonists in the reactively coping males show a reduction of immobility and an increase in escape attempts, whereas the proactively coping animal shows a decrease in escape and an increase in 64 29 immobility behavior in the test . Koolhaas et al . concluded that the effect of 5-HT may be related to the coping styles axis rather than to the emotionality axis. This idea is supported by experimental evidence obtained in the current study, where we show that 5-HT may be involved in regulation of behavior not only along a stress reactiv- ity gradient but also along a coping styles axis. These findings indicate that personality types mediate the strength and direction of 5-HT treatment, which should be taken into account in research on behavior and the underlying metabolism, the influence of pharmacological compounds on behavioral responses and developing new types of antidepressants. The proximate task of future research would be to reveal the types of 5-HT receptors that are involved in cricket behavioral responses and whether different receptor subtypes are engaged in stress-induced redistribution of immune reactions during behavioral responses. It is important to explain the link between developmental speed and body size found in this study, since it implies greater complexity in the structure and development of personality than expected before . We found that slowly developing crickets were significantly larger than rapidly developing crickets, which is in line with earlier SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 7 www.nature.com/scientificreports/ data on this species . Although faster developmental time in rapidly developing individuals may explain their smaller body size, rapidly developing, oen b ft old crickets are less anxious and can allocate more time to feeding, thus compensating for shorter time frames left for growth. The fact that the SSRI treatment ae ff cted slowly and rapidly developing crickets in the opposite way may indicate that rapidly developing crickets have higher baseline 5-HT in their brain than slowly developing crickets. In D. melanogaster, neuromodulatory actions of 5-HT were shown to depress feeding, while decreased neu- ronal 5-HT levels increased appetite . In another dipteran species, the flesh fly Neobellieria bullata , 5-HT injec- tion in the hemolymph decreased feeding . 5-HT and its effects on feeding are necessary to induce behavioral gregarization in migratory desert locusts . The decrease in body size is among the most obvious effects of phase 67,68 transition in morphological appearance from solitary to gregarious forms of locusts . Consistent with these findings, reduction in the brain 5-HT signaling by injection of a serotonergic-specific neurotoxin induces hyper - phagia and increased growth in rats . This may explain the larger body size of the slowly developing, oen s ft hy crickets. The observation that D. melanogaster ns3 mutants consume the same amount of food as control flies during early larval stages also suggests that serotonergic signaling appears to coordinate feeding behavior with 62 62 growth at the organismal level . Importantly, Kaplan et al. found close apposition of the axonal projections of serotonergic neurons to the insulin-producing cells (IPCs) in the brain of D. melanogaster, suggesting that these cells communicate to coordinate growth with 5-HT-modulated behaviors and environment. They showed that in D. melanogaster, insulin-like peptides both regulate blood sugar and act as growth factors, where serotonergic 62,70 neurons control the adult body size by ae ff cting insulin-like peptide secretion . This calls for studies on possi- ble similar synergies in cricket brains and bodies to define the functional patterns of the serotonergic synapse and subsequent effects on the phenotype. This may help to understand how and why 5-HT regulates somatic growth and the development of personality. For example, behaviors and sensitivity of 5-HT receptors might depend not only on the concentration of 5-HT but also on the proximity of the IPCs and serotonergic neurons in the brain. In conclusion, we show that personality-indicating consistent behavioral responses in crickets are tightly linked to developmental speed and body size. Antidepressants ae ff ct exploration and stress reactivity in different ways: we observed a higher SSRI impact on the coping styles axis than on the stress reactivity axis. Moreover, we found the opposite type of reactions to SSRI treatment in slowly and rapidly developing crickets along the coping styles axis. On the other hand, there were no simple one-way decreases or increases in exploration time at the two ends of the bimodal distribution of developmental speed / behavioral responses. This study shows that the SSRI treatment significantly decreased the variability of behaviors along the coping styles and stress reactivity axes, thus suggesting that antidepressants reduce variation in personality-related behaviors. It is also important to explain the link between developmental speed and body size found in this study, since it suggests greater com- plexity in the structure and development of personality than expected before . Overall, our results may have lim- ited explanatory power with regard to the magnitude of the SSRI effects on behaviors. This is because we have not studied SSRI activity in a dose-dependent manner, nor do we have data on possible side effects of SSRI treatment on cricket behavior(s). However, the high repeatability values that we recorded on the observed behaviors along the coping styles and the stress reactivity axes have the potential to show the diversity of directions that SSRI effects can take, and how these vary based on developmental speed and personality of the animals. Taken together, the present findings call for more research in brain neurophysiology and receptor distribution, a better definition of synaptic function, involvement of more personality axes and personality types (not only shy and bold ones), and an improved understanding of the brain concentrations of 5-HT and dose-dependent 5-HT effects in animals. Finally, this study shows that slowly and rapidly developing crickets markedly differ in behav- ioral traits and metabolism. Slowly developing individuals were larger, more explorative, less bold and they had higher RMR and lower MMR compared to rapidly developing crickets before the SSRI treatment. Overall, these findings do not support recent ideas concerning the pace-of-life syndrome (POLS) hypothesis , which predict higher exploration and metabolic rates in rapidly developing bold animals. Methods Insects and selection lines. We studied crickets from a laboratory stock originating from a wild popula- th tion (Davis, California, USA). The animals were selected for developmental speed; crickets of the 5 generation were investigated in this study. To start the selection lines, a total of 903 juvenile wild crickets were used. Until the crickets reached maturity, they were kept individually in plastic containers (280 mm × 115 mm × 105 mm) with a ventilation hole (60 mm dia.) covered with a plastic net. They were kept in a 12:12 h light–dark cycle and 26 ± 1 °C, with food and water ad libitum before and during experimental trials. The animals were selected for rapid and slow development, and we also maintained control groups. For rapid and slow development lines, mated males and females were selected according to their maturation age, with maximum 16.67% of the most rapidly or slowly maturing individuals being used, while in the control line the matings were randomized over the whole maturation age range . The first three generations were reared at the University of Eastern Finland, while the rest of the rearing of the lines occurred in Estonia. After five generations of selection for developmental speed, the development time (the average maturation time ± SD) for rapidly developing individuals was 98.03 ± 11.72 days (n = 238 crickets), 120.28 ± 25.12 days (n = 284) for the control individuals and 138.41 ± 23.05 days (n = 269) for slowly developing crickets. The mean longevity after eclosion for all three groups was 88.50 ± 23.54 days (one-way ANOVA, F = 25.79, (2, 788) P < 0.001). Longevity of males and females of all lines did not die ff r signic fi antly (one-way ANOVA, F = 1.0, (1, 789) P = 0.318). Lifespan aer ft eclosion of rapidly developing crickets (83.95 ± 23.52 days) and the control individuals (84.93 ± 22.06 days) did not differ significantly (Tukey’s test, P > 0.05), while lifespan after eclosion of slowly developing crickets (96.61 ± 25.05 days) significantly differed from the mean longevities of rapidly developing and control crickets (Tukey’s tests, both P < 0.001). SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 8 www.nature.com/scientificreports/ Trials, measurements and antidepressants. We tested whether escitalopram, a selective serotonin reuptake inhibitor (SSRI) (Brand name: CIPRALEX, CFL Pharmaceuticals Ltd) has any effect on cricket behavior and metabolism. Escitalopram is the (S)-stereoisomer (Left-enantiomer) of the citalopram prescribed frequently 60,71–74 to humans, low concentrations of which have been detected in natural habitats . Before the SSRI was added to the food of crickets, coping style (exploration) tests (day 1), and stress reactivity tests, measurements of RMR and MMR (day 2) were carried out for each of the F5 generation crickets (slow, n = 37; rapid, n = 41; control, n = 28) 15 days aer m ft aturation at the imago stage. To assess repeatability, we carried out coping style tests, stress reactivity tests and measurements of RMR and MMR one more time on days 3 and 4. As soon as these trials were done, on day 5 we started adding the SSRI to the food of the crickets (0.5 mg/kg per day added to Tetra Fish Food Flakes: Tetra, Blacksburg, VA) for twelve days (days 5–16). On days 13 and 14, the behavioral tests (coping styles and stress reactivity) and metabolism measurements were done again using the same animals as used during the previous trials. To assess repeatability of the behavioral variables, we repeated the tests on days 15 and 16. e Th insects were weighed before each measurement of metabolic rate using a Kern analytical balance (Kern & Sohn GmbH, Balingen, Germany). RMR and MMR values were corrected for body mass . Behavioral trials were conducted under constant temperature (24 ± 1.5 °C) and soundproof conditions. We used dim red light (25 W red incandescent bulb) since Gryllus spp. cannot see long (red) wavelengths properly , which allowed us to observe the nocturnal insects without disturbing them. The crickets were provided with drinking water before the onset of the trials, while food was removed 3 h prior the beginning of experimental trials. All behavioral trials were performed between 17:00 and 24:00 h. e Th observers of insects’ behavior and metabolic responses were blind to all predictions and treatment . Tests under familiar non-stressful environment: coping styles (exploration). The shelter part and the food area were separated by a corridor within individual containers. e Th crickets were trained for a week to run through a corridor (3 cm wide and 10 cm long) to reach the food aer a 3 ft h fast. During the fast, water was provided ad libitum. We used a one-cent Euro coin as a novel object before the crickets received the SSRI, while a one-dime USD coin was used during the SSRI treatment. The novel objects were placed on the floor of the cor - ridor 5 cm away from the entrance. We randomly chose the side of the corridor in which to place a novel object, and we always placed it on the opposite side during the second trial. As soon as the door between the shelter and the corridor was opened and the animal started to move through the corridor, we recorded the time it took to reach the food area. Since the crickets did not make any stops in the corridor when a novel object was not present, the time spent in the corridor reflects exploration activity / coping style of each individual – the extent of paying attention to changes in the environment. Tests under unfamiliar environments: stress reactivity. Next day after the coping style test, we handled each cricket for 1 min and then gave it an opportunity to escape into a conical plastic Eppendorf test tube (volume 5 ml) kept in dark conditions, which was used as an insect chamber during the measurement of RMR. This insect chamber represented a novel and stressful environment to the crickets. Resuming activity thus reflected each individual’s level of anxiety in this stressful environment. We assumed that the higher the anxiety, the sooner the individual will resume active body movements to avoid the stressful conditions of the windy and 26,27 noisy insect chamber . The latency to resume activity indicates the duration of freezing or immobile state, a 77,78 widespread anti-predator response used by many taxa . The insect chamber was connected to the respirometer by means of rubber tubing . Upon escape to the insect chamber, each of the crickets immediately became com- pletely immobile, as if hiding in a burrow or a crevice. We waited until all crickets resumed active movements, recorded the time between entering the tube and resuming activity, and used it as a measure of stress reactivity. Measuring RMR and MMR. The insects remained in their metabolic chambers until their CO emissions (and therefore their RMRs) reached the lowest stable rates. We also measured MMR at times when the crickets’ movements reached the highest rates while they performed struggling movements. As soon as the measure- ments were over, we returned the crickets back to their plastic housing-boxes. The LI-7000 die ff rential CO /H O analyzer (LiCor, Lincoln, NE, USA) was calibrated by means of calibra- 2 2 12,13,26 tion gases with gas injection . While measuring CO emissions, the insect chamber was perfused with dry −1 CO -free air at a flow rate of 60 ml min . The respirometric device was combined with an infrared optical system 12,13 using IR-emitting and IR-sensor diodes . IR-diodes made it possible to record CO production and to follow the movement intensity of each cricket. Statistical analyses. We used repeated measures ANOVAs (rmANOVAs) to assess changes in body mass, coping style, stress reactivity, resting metabolic rate and maximum metabolic rate aer a ft ntidepressant treatment (SSRI, within-subjects factor) in the same subjects with developmental line and sex used as between-subjects factors. We reported only the main effects when no significant interactions between between-subjects factors and/or interactions between within-subjects factor (antidepressant treatment) and between-subjects factors were found; otherwise the simple main effects using post hoc Bonferroni test were also reported. In order to account for age and testing apparatus effects, we performed rmANOVA with untreated (no antidepressant was given) groups of crickets (slow, n = 20; rapid, n = 20; control, n = 20). Analyses were performed in Statistica 8.0 for Windows (StatSoft Inc., Tulsa, OK, USA). To assess repeatability, we calculated the intraclass correlation coeffi- cient using a two-way random effects model for single measures . Pearson correlations were used to assess asso- ciations between different traits. Lower and upper 95% CI limits of correlation coefficient were calculated using non-stratified bootstrapping sampling with 2000 iterations. SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 9 www.nature.com/scientificreports/ References 1. Kain, J. S., Stokes, C. & de Bivort, B. L. Phototactic personality in fruit flies and its suppression by serotonin and white. Proc. Natl. Acad. Sci. USA 109, 19834–19839 (2012). 2. Buchanan, S. M., Kain, J. S. & de Bivort, B. L. Neuronal control of locomotor handedness in Drosophila. Proc. Natl. Acad. Sci. USA 112, 6700–6705 (2015). 3. Verbeek, M. E. M., Drent, P. J. & Wiepkema, P. R. Consistent individual differences in early exploratory behavior of male great tits. Anim. Behav. 48, 1113–1121 (1994). 4. Drent, P. J., van Oers, K. & van Noordwijk, A. J. Realized heritability of personalities in the great tit (Parus major). Proc. Biol. Sci. 270, 45–51 (2003). 5. Wolf, M., van Doorn, G. S., Leimar, O. & Weissing, F. J. Life-history trade-offs favour the evolution of animal personalities. Nature 447, 581–584 (2007). 6. Wolf, M., van Doorn, G. S. & Weissing, F. J. Evolutionary emergence of responsive and unresponsive personalities. Proc. Natl. Acad. Sci. USA 105, 15825–15830 (2008). 7. Kortet, R., Vainikka, A., Janhunen, M., Piironen, J. & Hyvärinen, P. Behavioral variation shows heritability in juvenile brown trout Salmo trutta. Behav. Ecol. Sociobiol. 6, 927–934 (2014). 8. Wolf, M. & Weissing, F. J. Animal personalities: consequences for ecology and evolution. Trends Ecol. Evol. 27, 452–461 (2012). 9. Smith, B. R. & Blumstein, D. T. Behavioral types as predictors of survival in Trinidadian guppies (Poecilia reticulate). Behav. Ecol. 21, 919–926 (2010). 10. Kortet, R., Hedrick, A. V. & Vainikka, A. Parasitism, predation and the evolution of animal personalities. Ecol Lett. 13, 1449–1458 (2010). 11. Dall, S. R. X., Bell, A. M., Bolnick, D. I. & Ratnieks, F. L. W. An evolution of ecology of individual differences. Ecol. Lett. 15, 1189–1198 (2012). 12. Krams, I. et al. Predation selects for low resting metabolic rate and consistent individual differences in anti-predator behavior in a beetle. Acta Ethol. 16, 163–172 (2013). 13. Krams, I. et al. Predation promotes survival of beetles with lower resting metabolic rates. Entomol Exp. Appl. 148, 94–103 (2013). 14. MacArthur, R. H. & Wilson, E. O. e Th Th eory of Island Biogeography (Princeton University Press, 1967). 15. Pianka, E. R. On r- and K-selection. Am. Nat. 104, 592–597 (1970). 16. Réale, D. et al. Personality and the emergence of the pace-of-life syndrome concept at the population level. Philos. Trans. R. Soc. London B Biol. Sci. 365, 4051–4063 (2010). 17. Biro, P. A. & Stamps, J. A. Are animal personality traits linked to life-history productivity? Trends Ecol. Evol. 23, 361–368 (2008). 18. Sih, A., Bell, A. M., Johnson, J. C. & Ziemba, R. E. Behavioral syndromes: an integrative overview. Q. Rev. Biol. 79, 241–277 (2004). 19. Réale, D., Reader, S. M., Sol, D., McDougall, P. T. & Dingemanse, N. J. Integrating animal temperament within ecology and evolution. Biol. Rev. Camb. Philos. Soc. 82, 291–318 (2007). 20. Steimer, T., La Fleur, S. & Schulz, P. E. Neuroendocrine correlates of emotional reactivity and coping in male rats from the Roman high (RHA/Verh)- and low (RLA/Verh)-avoidance lines. Behav. Genet. 27, 503–512 (1997). 21. Benus, R. F., Koolhaas, J. M. & van Oortmerssen, G. A. Aggression and adaptation to the light–dark cycle: role of intrinsic and extrinsic control. Physiol. Behav. 43, 131–137 (1988). 22. Bolhuis, J. E., Schouten, W. G., de Leeuw, J. A., Schrama, J. W. & Wiegant, V. M. Individual coping characteristics, rearing conditions and behavioural flexibility in pigs. Behav. Brain. Res. 152, 351–360 (2004). 23. West, G. B., Brown, J. H. & Enquist, B. J. A general model for ontogenetic growth. Nature 413, 628–631 (2001). 24. Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Metabolic theory of ecology. Ecology 85, 1771–1789 (2004). 25. Speakman, J. R. Doubly Labelled Water: Theory and Practice (Chapman & Hall, 1997). 26. Lighton, J. R. B. Measuring Metabolic Rate: a Manual for Scientists (Oxford University Press, 2008), 27. Careau, V., Thomas, D., Humphries, M. M. & Réale, D. Energy metabolism and animal personality. Oikos 117, 641–653 (2008). 28. Krams, I. A. et al. Metabolic rate associates with, but does not generate covariation between, behaviours in western-trilling crickets, Gryllus integer. Proc. Biol. Sci. 284, e20162481 (2017). 29. Koolhaas, J. M., de Boer, S. F., Coppens, C. M. & Buwalda, B. Neuroendocrinology of coping styles: Towards understanding the biology of individual variation. Front. Neuroendocrinol. 31, 307–321 (2010). 30. Seltmann, M. W. et al. Stress responsiveness, age and body condition interactively ae ff ct flight initiation distance in breeding female eiders. Anim. Behav. 84, 889–896 (2012). 31. Koolhaas., J. M. et al. Coping styles in animals: current status in behavior and stress-physiology. Neurosci. Biobehav. Rev. 23, 925–935 (1999). 32. Gelder, M., Mayou, R. & Geddes, J. Psychiatry. (Oxford University Press, 2005). 33. Kravitz, E. A. & Huber, R. Aggression in invertebrates. Curr. Opin. Neurobiol. 13, 736–743 (2003). 34. Anstey, M. L., Rogers, S. M., Ott, S. R., Burrows, M. & Simpson, S. J. Serotonin mediates behavioral gregarization underlying swarm formation in desert locusts. Science 323, 627–630 (2009). 35. Donovan, M. H. & Tecott, L. H. Serotonin and the regulation of mammalian energy balance. Front. Neurosci. 7, e36 (2013). 36. De Boer, S. F. & Newman-Tancredi, A. Anti-aggressive effects of the selective high-efficacy ‘biased’ 5-HT1A receptor agonists F15599 and F13714 in male WTG rats. Psychopharmacology 233, 937–947 (2016). 37. Dolan, R. J., Morris, J. S. & de Gelder, B. Crossmodal binding of fear in voice and face. Proc. Natl. Acad. Sci. USA 98, 10006–10010 (2001). 38. Murphy, D. L. & Lesch, K. P. Targeting the murine serotonin transporter: insights into human neurobiology. Nat. Rev. Neurosci. 9, 85–96 (2008). 39. Caramaschi, C., De Boer, S. F. & Koolhaas, J. M. Differential role of the 5-HT1A receptor in aggressive and non-aggressive mice: an across-strain comparison. Physiol. Behav. 90, 590–601 (2007). 40. Steimer, T. & Driscoll, P. Divergent stress responses and coping styles in psychogenetically selected Roman high-(RHA) and low- (RLA) avoidance rats: behavioural, neuroendocrine and developmental aspects. Stress 6, 87–100 (2003). 41. Niemelä, P. T., Vainikka, A., Lahdenperä, S. & Kortet, R. Nymphal density, behavioural development, and life history in a field cricket. Behav. Ecol. Sociobiol. 66, 645–652 (2012). 42. Niemelä, P. T., Vainikka, A., Hedrick, A. V. & Kortet, R. Integrating behaviour with life history: boldness of the field cricket, Gryllus integer, during ontogeny. Funct. Ecol. 26, 450–456 (2012). 43. Niemelä, P. T., Dingemanse, N., Alioravainen, N., Vainikka, A. & Kortet, R. Personality pace-of-life hypothesis: testing genetic associations among personality and life-history. Behav. Ecol. 24, 935–941 (2013). 44. Rothwell, N. J. & Stock, M. J. Effect of diet and fenfluramine on thermogenesis in the rat: possible involvement of serotonergic mechanisms. Int. J. Obes. 11, 319–324 (1987). 45. Le Feuvre, R. A., Aisenthal, L. & Rothwell, N. J. Involvement of corticotrophin releasing factor (CRF) in the thermogenic and anorexic actions of serotonin (5-HT) and related compounds. Brain Res. 555, 245–250 (1991). 46. Laporta, J. et al. Increasing serotonin concentrations alter calcium and energy metabolism in dairy cows. J. Endocrinol. 226, 43–55 (2015). SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 10 www.nature.com/scientificreports/ 47. Lin, M. T., Pang, I. H., Chern, S. I. & Chia, W. Y. Changes in serotonin contents in brain ae ff ct metabolic heat production of rabbits in cold. Am. J. Physiol. 235, R41–47 (1978). 48. Nonogaki, K., Abdallah, L., Goulding, E. H., Bonasera, S. J. & Tecott, L. H. Hyperactivity and reduced energy cost of physical activity in serotonin 5-HT2C receptor mutant mice. Diabetes 52, 315–320 (2003). 49. Weng, R. et al. Metabolomics Approach Reveals Integrated Metabolic Network Associated with Serotonin Deficiency. Sci. Rep. 5, e11864 (2015). 50. Klein, D. F. Delineation of two drug-responsive anxiety syndromes. Psychopharmacologia 5, 397–408 (1964). 51. Graeff, F. G. & Del-Ben, C. M. Neurobiology of panic disorder: from animal models to brain neuroimaging. Neurosci. Biobehav. Rev. 32, 1326–1335 (2008). 52. De Boer, S. F. & Koolhaas, J. M. 5-HT1A and 5-HT1B receptor agonists and aggression: a pharmacological challenge of the serotonin deficiency hypothesis. Eur. J. Pharmacol. 526, 125–139 (2005). 53. Gutiérrez-Adán, A. et al. Effect of speed of development on mRNA expression pattern in early bovine embryos cultured in vivo or in vitro. Mol. Reprod. Dev. 68, 441–448 (2004). 54. Chung, J. Y. Y. & Evans, M. A. Shyness and symptoms of illness in young children. Can. J. Behav. Sci. 32, 49–57 (2000). 55. Rettew, D. C. & McKee, L. Temperament and its role in developmental psychopathology. Harv. Rev. Psychiatry 13, 14–27 (2005). 56. Rantala, M. J., Luoto, S., Krams, I. & Karlsson, H. Depression subtyping based on evolutionary psychiatry: Proximate mechanisms and ultimate functions. Brain Behav. Immun. 69, 603–617 (2018). 57. Koren-Schwartzer, N., Chen-Zion, M., Ben-Porat, H. & Beitner, R. Serotonin-induced decrease in brain ATP, stimulation of brain anaerobic glycolysis and elevation of plasma hemoglobin; the protective action of calmodulin antagonists. Gen. Pharmacol. 25, 1257–1262 (1994). 58. Gardner, A. & Boles, R. G. Beyond the serotonin hypothesis: mitochondria, inflammation and neurodegeneration in major depression and ae ff ctive spectrum disorders. Prog. Neuropsychopharmacol. Biol. Psychiatry 35, 730–743 (2011). 59. Volpi-Abadie, J., Kaye, A. M. & Kaye, A. D. Serotonin syndrome. Ochsner J. 13, 533–540 (2013). 60. Richmond, E. K., Rosi-Marshall, E. J., Lee, S. S., Thompson, R. M. & Grace, M. R. Antidepressants in stream ecosystems: influence of selective serotonin reuptake inhibitors (SSRIs) on algal production and insect emergence. Freshw. Sci. 35, 845–855 (2016). 61. Niemelä, P. T., Vainikka, A., Forsman, J. T., Loukola, O. J. & Kortet, R. How does variation in the environment and individual cognition explain the existence of consistent behavioural differences? Ecol. Evol. 3, 457–464 (2013). 62. Kaplan, D. D., Zimmermann, G., Suyama, K., Meyer, T. & Scott, M. P. A nucleostemin family GTPase, NS3, acts in serotonergic neurons to regulate insulin signaling and control body size. Genes Dev. 22, 1877–1893 (2008). 63. Veenema, A. H. et al. Differences in the effects of 5-HT(1A) receptor agonists on forced swimming behavior and brain 5-HT metabolism between low and high aggressive mice. Psychopharmacology 178, 151–160 (2005). 64. Cryan, J. F., Redmond, A. M., Kelly, J. P. & Leonard, B. E. The effects of the 5-HT1A agonist flesinoxan, in three paradigms for assessing antidepressant potential in the rat. Eur. Neuropsychopharmacol. 7, 109–114 (1997). 65. Neckameyer, W. S. A trophic role for serotonin in the development of a simple feeding circuit. Dev. Neurosci. 32, 217–237 (2010). 66. Dacks, A. M., Nickel, T. & Mitchell, B. K. An examination of serotonin and feeding in the flesh fly Neobellieria bullata (Sarcophagidae: Diptera). J. Insect. Behav. 16, 1–21 (2003). 67. Pener, M. & Simpson, S. Locust phase polyphenism: an update. Adv. Insect. Physiol. 36, 1–272 (2009). 68. Uvarov, B. P. Grasshoppers and Locusts (Cambridge University Press, 1966). 69. Saller, C. F. & Stricker, E. M. Hyperphagia and increased growth in rats aer in ft traventricular injection of 5,7-dihydroxytryptamine. Science 192, 385–387 (1976). 70. Ruaud, A. F. & Thummel, C. S. Serotonin and insulin signaling team up to control growth in Drosophila. Genes Dev. 22, 1851–1855 (2008). 71. Schultz, M. M. & Furlong, E. T. Trace analysis of antidepressant pharmaceuticals and their select degradates in aquatic matrixes by LC/ESI/MS/MS. Anal. Chem. 80, 1756–1762 (2008). 72. Metcalfe, C. D. et al. Antidepressants and their metabolites in municipal wastewater, and downstream exposure in an urban watershed. Environ. Toxicol. Chem. 29, 79–89 (2010). 73. Cardoso, O., Porcher, J.-M. & Sanchez, W. Factory-discharged pharmaceuticals could be a relevant source of aquatic environment contamination: review of evidence and need for knowledge. Chemosphere 115, 20–30 (2014). 74. Larsson, D. G. J. Pollution from drug manufacturing: review and perspectives. Philos. Trans. R. Soc. Lond. B Biol. Sci. 369, e20130571 (2014). 75. Briscoe, A. D. & Chittka, L. The evolution of color vision in insects. Ann. Rev. Entomol. 46, 471–510 (2001). 76. Burghardt, G. M. et al. Perspectives – minimizing observer bias in behavioral studies: a review and recommendations. Ethology 118, 511–517 (2012). 77. Chelini, M.-C., Willemart, R. H. & Hebets, E. A. Costs and benefits of freezing behavior in the harvestman Eumesosoma roeweri (Arachnida, Opiliones). Behav. Processes. 82, 153–159 (2009). 78. Krams, I. et al. High repeatability of anti-predator responses and resting metabolic rate in a beetle. J. Insect. Behav. 27, 57–66 (2014). 79. Mattila, A. L. K. & Hanski, I. Heritability of flight and resting metabolic rates in the Glanville fritillary butterfly. J. Evol. Biol. 27, 1733–1743 (2014). 80. McGraw, K. O. & Wong, S. P. Forming inferences about some intraclass correlation coefficients. Psychol. Methods 1, 30–46 (1996). Acknowledgements We thank Ilkka Sirkka, Sanni Hietaketo, Leena Koponen, Leena Pääkkönen, Julia Keronen, Sirpa Kaunisto, Anssi Vainikka, Inese Kivleniece, Jolanta Vrublevska and Ljubova Sivačova for assistance, and Raivo Mänd for comments on the manuscript. Author Contributions I.K., G.T., Ra.K. and T.K. conceived the study. I.K., A.K., R.K., Ra.K., T.K., R.K., M.M., M.J.R., J.K., P.J. and S.K. contributed to the design of the study and carried out the experiments. D.E., G.T. and I.K. analyzed the data. I.K., S.L. and G.T. wrote the manuscript. All authors participated in results interpretation and they approved the final manuscript. Additional Information Competing Interests: The authors declare no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 11 www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2018 SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 12 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

Linking organismal growth, coping styles, stress reactivity, and metabolism via responses against a selective serotonin reuptake inhibitor in an insect

Free
12 pages
Loading next page...
 
/lp/springer_journal/linking-organismal-growth-coping-styles-stress-reactivity-and-9MOf49d4eK
Publisher
Nature Publishing Group UK
Copyright
Copyright © 2018 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
D.O.I.
10.1038/s41598-018-26722-9
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports OPEN Linking organismal growth, coping styles, stress reactivity, and metabolism via responses against Received: 9 October 2017 a selective serotonin reuptake Accepted: 15 May 2018 Published: xx xx xxxx inhibitor in an insect 1,2,3 4,5 5 6 5 Indrikis Krams , Giedrius Trakimas , Sanita Kecko , Didzis Elferts , Ronalds Krams , 7,8 9 10 10 11 Severi Luoto , Markus J. Rantala , Marika Mänd , Aare Kuusik , Jukka Kekäläinen , 12 11 5,10 Priit Jõers , Raine Kortet & Tatjana Krama Evidence suggests that brain serotonin (5-HT) is one of the central mediators of different types of animal personality. We tested this assumption in field crickets Gryllus integer using a selective serotonin reuptake inhibitor (SSRI). Crickets were selected for slow and rapid development and tested for their coping styles under non-stressful conditions (time spent exploring a novel object). Resting metabolic rate, maximum metabolic rate and latency to resume activity were measured under stressful conditions (stress reactivity). Measurements were taken (i) before and (ii) during the SSRI treatment. Before the SSRI treatment, a strong negative correlation was observed between coping style and stress reactivity, which suggests the existence of a behavioral syndrome. After the SSRI treatment, the syndrome was no longer evident. The results of this study show that 5-HT may be involved in regulating behavior not only along a stress reactivity gradient but also along a coping styles axis. The relationship between personality and the strength and direction of 5-HT treatment on observed behaviors indicates trait-like individual differences in 5-HT signaling. Overall, these findings do not support recent ideas arising from the pace-of-life syndrome (POLS) hypothesis, which predict higher exploration and metabolic rates in rapidly developing bold animals. Phenotypic differences in trait expression are a prerequisite for evolution through natural selection. Although most of the traits of an organism vary around a mean value, research on human and animal behavior has revealed between-individual clusters in behavioral responses, oen a ft ttributed to different personality axes that are called personality types, behavioral predispositions, temperaments or coping styles. Some of these behavioral responses are 1,2 3–7 not heritable , while others demonstrate a negligible or a considerable degree of heritability . Animal person- 5,6,8–13 ality is tightly linked to individual life histories, survival and fitness . e Th pace-of-life syndrome (POLS) hypothesis, rooted in the classic concept of r- and K-selection, is particu- 14–16 larly relevant in this context , as it suggests that life history traits such as age at maturity and growth rate are likely coupled with individually consistent behaviors and immune responses. For example, shy individuals oen ft 1 2 Institute of Ecology and Earth Sciences, University of Tartu, Tartu, Estonia. Department of Zoology and Animal Ecology, Faculty of Biology, University of Latvia, Rīga, Latvia. University of Tennessee, Department of Psychology, 4 5 Knoxville, USA. Institute of Biosciences, Vilnius University, Vilnius, Lithuania. Department of Biotechnology, Institute of Life Sciences and Technology, Daugavpils University, Daugavpils, Latvia. Department of Botany and Ecology, Faculty of Biology, University of Latvia, Rīga, Latvia. English, Drama and Writing Studies, University of Auckland, Auckland, New Zealand. School of Psychology, University of Auckland, Auckland, New Zealand. 9 10 Department of Biology & Turku Brain and Mind Centre, University of Turku, Turku, Finland. Department of Plant Protection, Institute of Agricultural and Environmental Sciences, Estonian University of Life Science, Tartu, Estonia. 11 12 Department of Environmental and Biological Sciences, University of Eastern Finland, Joensuu, Finland. Insttute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia. Correspondence and requests for materials should be addressed to I.K. (email: indrikis.krams@ut.ee) SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 1 www.nature.com/scientificreports/ grow slower, reach larger body size, mature later, invest more in immune system and have longer hiding time in an anti-predator context under familiar conditions than bold individuals . Aggressiveness of bold individuals hypothetically facilitates acquiring and monopolizing resources and may be linked to a fast growth rate and more 5,17 intense reproductive effort early in life . Personality research is a complex field that adjoins behavioral ecology with psychology and molecular biology. It is therefore important to define key terms before testable hypotheses are proposed. The term behavioral syn - 18,19 drome indicates that behavioral trait characteristics involve suites of correlated behaviors that are consistent over time and across contexts. The term coping style is more commonly used to characterize individual trait dif- ferences in biomedicine, and this concept has been useful in explaining behavioral responses and their underlying physiological mechanisms. The coping style concept implies that animals may react to a stimulus with alternative response patterns. According to the two-tier model of animal personality proposed by Steimer et al. , coping styles reflect the quality of the response to a stressor, while stress reactivity reflects the quantity of the response expressed. In this model, bold individuals are those with a strong tendency to act (proactivity) combined with low stress reactivity. Docile individuals can be characterized as a combination of reactive coping associated with low stress reactivity. The combination between high stress reactivity and a high tendency for proactive coping is labeled as panicky. Shy individuals exhibit reactive coping style and high stress reactivity. Although those personality characteristics are common, researchers rarely use all of the above terms because that would require the use of several tests for each individual being studied. For example, stress reactivity needs to be measured by obtaining levels of stress hormones/neuromediator concentrations, or by observing the organ- isms under stressful conditions. In contrast, reactive and proactive coping styles differ in behavioral flexibility, where the proactive individuals act on the basis of previous experiences and develop routines in stable environ- 21,22 ments . The reactively coping animal tends to rely more on a detailed appraisal of their current environment because they react to immediate environmental stimuli and tend to explore any changes in their environment. 23,24 The metabolic rate of an organism varies with its body size and body temperature . The POLS hypothesis proposes that rapidly developing, oen b ft old individuals have life histories involving high daily energy expendi- tures and higher basal or resting metabolic rates (RMR) because they are forced to process more food and inevi- tably excrete more waste products at a faster rate compared to slowly developing individuals of same size. Rapidly developing individuals may, however, remain smaller than slow developers . Thus, personality also has the potential to ae ff ct the metabolic rate of an organism. Another important but oen n ft eglected source of variability in RMR arises because individuals consistently differ in their stress responses to conditions in which metabolism 25,26 measurements are performed . For example, because reactive individuals (shy, docile) oen b ft ecome immobile 12,13 27 under stressful conditions , their metabolic rates may be mistakenly classified as RMR . Recent research shows that while the latency to resume body movements aer h ft andling in a familiar environment is longer for slowly developing crickets, they resume body movements sooner than rapidly developing crickets when in a more stress- ful, unfamiliar environment . Slowly developing crickets were found to have significantly higher RMR com- pared to rapidly developing crickets. It is important to note that according to established experimental procedure, measurements of respiration always require placing an organism in a novel environment . This has the potential to induce significant stress, which may force less stress-resistant animals to stay active during respirometry to explore or escape their respirometry chambers. This demonstrates the importance of n fi ding the proper approach to measurements of stress reactivity and metabolism, as well as making a correct interpretation of the results obtained to explain the relationships between animal personality and underlying physiological mechanisms . Behavioral responses of shy, oen a ft nxious individuals are associated with high concentrations of stress hor - mones , implying that their slowness and caution are caused by physiological stress. In psychology and physiol- ogy, shyness is usually defined as the feeling of anxiety and awkwardness when an individual is under conditions of stress. Stronger forms of shyness are usually referred to as anxiety and elevated neuroticism. In extreme cases, it might be related to a depression-like psychological condition, which includes the experience of fear to the extent 31 20,29 of inducing panic . This calls for reliable tests to discriminate between shy and panicky personality types . In the central nervous system, serotonin (5-hydroxytryptamine: 5-HT), a monoamine neurotransmitter , has 33–36 several important functions, including the regulation of aggression, mood, sleep, dispersal and metabolic rate . Modulation of serotonin at synapses of the brain is thought to be the key mechanism of several classes of pharma- cological antidepressants. Experimental evidence has revealed links between low levels of 5-HT and the develop- 37,38 ment of anxiety and stress-prone behavioral phenotypes . However, an association between low concentrations of 5-HT and proactive coping style suggests that 5-HT may be involved in the regulation of behavior not only 20,29,40 34 along the stress reactivity gradient but also along the coping styles axis (but see ). These findings indicate that 5-HT potentially ae ff cts personality types in a more complex way, which should be taken into account when testing the influence of pharmacological compounds on behavioral responses under stress. 41–43 Field crickets are a suitable and often-used model group for animal personality research . In this study, we compared coping style differences between slowly and rapidly developing western stutter-trilling crickets (Gryllus integer) under non-stressful and stressful conditions. Prior to the study, the crickets had been selected for slow and rapid developmental time over five generations. We also measured their RMR, maximum metabolic rate (MMR) and the latency to resume activity under stressful conditions in the insect chamber of respirometer. This was done twice: before a selective serotonin reuptake inhibitor (SSRI) was added to the food and once again during the anti-depressant treatment. We expected greater effects of anti-depressant treatment on behavioral responses of proactive individuals (bold, panicky) − those paying little attention to minor changes in their environment. Although some earlier evidence suggested that central injection of 5-HT or its precursor 5-hydroxytryptophan (5-HTP) may cause 44,45 46 increased metabolic rates (but see ), other research has shown that serotonin deficiency significantly increases 47–49 amino acid, energy, purine, lipid and gut microflora metabolisms and rates of oxidative stress . Therefore, we predicted an SSRI-induced decrease in the metabolic rate of crickets with high stress reactivity (shy, panicky), SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 2 www.nature.com/scientificreports/ Figure 1. Mean (±95%CI) body mass (a), coping styles (b), time to resume activity (c), RMR (d) and MMR (e) before SSRI (▫, ◦) and aer SS ft RI treatment (▪, ⦁) in crickets of rapid (n = 41), slow (n = 37) and control (n = 28) developmental lines. ◦, ⦁ indicate females, ▫, ▪ indicate males. 50,51 especially in panicky individuals , because 5-HT is likely to influence proactive individuals (panicky, bold) 36,52 more than reactive ones (shy, docile) . Results Body mass. Our results demonstrate that SSRI can have relatively small but significant effects on body mass of crickets (rmANOVA, F = 22.27, P < 0.0001, partial η = 0.18, Figs 1a and 2a). The average body mass (1, 100) of crickets decreased from 0.548 ± 0.004 g to 0.539 ± 0.004 g (unweighted means ± SE) during SSRI treatment. e Th body mass of crickets largely depended on the developmental line (rmANOVA, F = 356.9, P < 0.0001, (2, 100) partial η = 0.85). The effect of sex, although much smaller than that of developmental line, was significant (rmANOVA, F = 5.41, P = 0.02, partial η = 0.05). Before the SSRI treatment, slowly developing crickets (1, 100) were heavier (0.701 ± 0.006 g (weighted means ± SE) than rapidly developing (0.466 ± 0.007 g) and control crick- ets (0.466 ± 0.011 g), while females (0.558 ± 0.016 g) were heavier than males (0.537 ± 0.017 g). Although body mass decreased in males and females of all developmental lines during the SSRI treatment, slowly developing crickets remained the heaviest among the lines, while females were heavier than males. All possible interac- tions between factors were non-significant. No aging and testing apparatus effects were evident (rmANOVA, F (1, = 1.0, P = 0.323, partial η = 0.02). Body mass was repeatable in control (R = 0.69) and highly repeatable in 54) slowly (R = 0.93) and rapidly developing crickets (R = 0.94; Table 1). Coping style. We found significant main effects of developmental lines (rmANOVA, F = 70.45, (2, 100) 2 2 P < 0.0001, partial η = 0.56) and SSRI treatment (rmANOVA, F = 10.19, P = 0.002, partial η = 0.09) to (1, 100) coping styles of the crickets, while sex apparently had no effect on it (rmANOVA, F = 0.1, P = 0.76, partial (1, 100) η = 0.001) (Figs 1b and 2b). However, there was significant interaction between developmental line and antide- pressant treatment (F = 42.03, P < 0.0001, partial η = 0.46), suggesting that simultaneous ee ff cts of those (2, 100) factors on the coping style were not additive. Before the SSRI treatment, slowly developing crickets explored the novel object longer than control (Bonferroni test: P = 0.0002) and rapidly developing crickets (Bonferroni test: P < 0.0001); control crickets explored the object longer than rapidly developing individuals (Bonferroni test: P < 0.0001). During the SSRI treatment, rapidly developing crickets decreased their proactivity by increas- ing exploration time of the novel object (Bonferroni test: P = 0.0008). Slowly developing crickets decreased the proportion of time devoted to exploration (Bonferroni test: P < 0.0001) (Fig. 3), while the control individuals explored the novel object equally long (Bonferroni test: P > 0.05). Aer t ft he SSRI treatment, coping styles of all developmental lines became statistically similar (Bonferroni test: P > 0.05). No aging nor testing apparatus effects were evident (rmANOVA, F = 0.93, P = 0.339, partial η = 0.02). Coping styles were repeatable in the control (1, 54) (R = 0.91) and slowly developing crickets (R = 0.67) but not in the rapidly developing crickets (R = 0.16; Table 1). However, combined repeatability was high (R = 0.94) when calculated for all selection lines together. SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 3 www.nature.com/scientificreports/ Figure 2. Relative change in mean body mass (a), coping styles (b), time to resume activity (c), RMR (d) and MMR (e) before SSRI and aer SS ft RI treatment in crickets of rapid, slow and control developmental lines. Significance based on Tukey’s test: n.s. P > 0.05, *P < 0.05, **P < 0.01, ***P < 0.001. Control Slow Rapid Variables R 95%CI limits F R 95%CI limits F R 95%CI limits F 19,19 19,19 19,19 *** **** **** Body mass 0.69 0.38; 0.86 5.75 0.93 0.82; 0.97 25.1 0.94 0.85; 0.98 29.8 a **** *** ns Coping styles 0.91 0.80; 0.97 21.4 0.67 0.33; 0.85 4.95 0.16 −0.15; 0.50 1.58 **** **** **** Stress reactivity 0.99 0.98; 1.00 226.3 0.92 0.82; 0.97 25.5 0.97 0.93; 0.99 66.4 **** **** **** RMR 0.96 0.90; 0.98 51.1 0.91 0.80; 0.97 22.9 0.79 0.51; 0.92 10.2 **** **** **** MMR 0.94 0.87; 0.98 33.7 0.94 0.86; 0.98 33.5 0.94 0.86; 0.98 32.4 Table 1. Repeatabilities of traits in crickets of rapid (n = 41), slow (n = 37) and control (n = 28) developmental *** **** ns lines estimated using two-way random effects model for single measures. P < 0.001; P < 0.0001; P > 0.05; combined all selection lines: R = 0.94 (0.91; 0.97), F = 33.8, P < 0.001. 59,59 Stress reactivity. We found significant main effects of antidepressant treatment on time to resume activ- ity in the insect chamber (rmANOVA, F = 33.83, P < 0.0001, partial η = 0.25) (Figs 1c and 2c). Time to (1, 100) resume activity increased from 221 ± 12 s (unweighted means ± SE) before the SSRI treatment to 284 ± 11 s aer ft the treatment. Stress reactivity of the crickets depended on developmental line (rmANOVA, F = 61.73, (2, 100) P < 0.0001, partial η = 0.55), while sex had no effect on stress reactivity (rmANOVA, F = 0.12, P = 0.73, (1, 100) partial η = 0.001) either before or after the SSRI treatment. After adding antidepressants to the food, all selec- tion lines increased their time to resume active movements in the insect chamber. No significant interactions were observed. No aging and testing apparatus effects were found (rmANOVA, F = 0.41, P = 0.523, partial (1, 54) η = 0.01). Stress reactivity was highly repeatable in all selection lines (range R = 0.92–0.99; Table 1). Resting metabolic rate. A significant main effect of developmental line on RMR was observed (rmANOVA, F = 34.63, P < 0.0001, partial η = 0.41) (Figs 1d and 2d). However, there was a significant (2, 100) interaction between developmental line and antidepressant treatment (rmANOVA, F = 25.52, P < 0.0001, (2, 100) partial η = 0.38). Before the SSRI treatment, the RMR of slowly developing crickets was higher than that of control (Bonferroni test: P < 0.0001) and rapidly developing crickets (Bonferroni test: P < 0.0001). The RMR of rapidly developing crickets was lower than that of control crickets (Bonferroni test: P = 0.008). After the treat- ment, RMR significantly increased in rapidly developing crickets (P = 0.001), significantly decreased in slowly developing crickets (Bonferroni test: P < 0.0001) (Fig. 4), but remained the same in controls (Bonferroni test: P = 1). Aer t ft he antidepressant treatment, the RMR of slowly developing crickets was statistically similar to the RMR in the control line before the SSRI treatment (Bonferroni test: P = 0.29). As the mean RMR data show, all the SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 4 www.nature.com/scientificreports/ Figure 3. Change in coping styles aer SS ft RI treatment in individual crickets of rapid (──, n = 41) and slow (……, n = 37) developmental lines (a) and control (n = 28) line (b). Figure 4. Change in RMR aer SS ft RI treatment in individual crickets of rapid (──, n = 41) and slow (……, n = 37) developmental lines (a) and control (n = 28) line (b). lines became more similar aer t ft he SSRI treatment. However, some of the developmental lines remained signifi- cantly different from each other (i.e., rapidly developing crickets vs slowly developing crickets: P = 0.004; slowly developing crickets vs control group: P = 0.02). RMR was repeatable in rapidly developing (R = 0.79), and highly repeatable in slowly developing (R = 0.91) and control crickets (R = 0.96; Table 1). Maximum metabolic rate. A small significant main effect of developmental line on MMR (rmANOVA, F = 4.03, P = 0.02, par tial η = 0.08) was observed, while sex had no effect on MMR (rmANOVA, (2, 100) F = 2.26, P = 0.07, partial η = 0.03) (Figs 1e,2e). However, there was a significant interaction between (1, 100) developmental line and SSRI treatment (rmANOVA, F = 64.38, P < 0.0001, partial η = 0.56). Before the (2, 100) SSRI treatment, slowly developing crickets had the highest rates of maximum metabolism, followed by control and rapidly developing crickets. Significant differences were found between rapidly and slowly developing lines (Bonferroni test: P < 0.0001) and between rapidly developing and control crickets (Bonferroni test: P = 0.006), while slowly developing crickets and the control group did not differ (Bonferroni test: P = 1). During the SSRI treatment, MMR significantly increased in rapidly developing crickets (Bonferroni test: P < 0.0001) and decreased in slowly developing individuals (Bonferroni test: P < 0.0001), while it did not change significantly in the control group (Bonferroni test: P = 1). Aer t ft he SSRI treatment, the MMR of all developmental lines was sta- tistically similar (Bonferroni test: P = 1). MMR was highly repeatable in all selection lines (all R = 0.94; Table 1). Associations among traits. Before the SSRI treatment, the crickets’ body mass was positively correlated with RMR (r = 0.63, P < 0.0001), MMR (r = 0.34, P < 0.01) and coping styles (r = 0.69, P < 0.0001), while being SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 5 www.nature.com/scientificreports/ Traits Body mass Coping styles Stress reactivity RMR MMR **** *** −0.14 −0.63 0.40 −0.08 Body mass (−0.34; 0.06) (−0.73; −0.53) (0.22; 0.55) (−0.29; 0.13) **** ** 0.69 0.00 0.33 0.09 Coping styles (0.59; 0.78) (−0.17; 0.17) (0.8; 0.53) (−0.09; 0.26) −0.32 **** **** −0.66 −0.60 −0.03 Stress reactivity (−0.48; (−0.74; −0.57) (−0.73; −0.47) (−0.21; 0.16) −0.15) **** **** **** 0.63 0.69 −0.74 0.09 RMR (0.50; 0.74) (0.56; 0.79) (−0.80; −0.67) (−0.10; 0.27) ** ** *** *** 0.34 0.37 −0.40 0.40 MMR (0.13; 0.52) (0.20; 0.53) (−0.55; −0.24) (0.21; 0.56) Table 2. Pearson correlation coefficients between traits in crickets before SSRI treatment (below diagonal) and aer SS ft RI treatment (above diagonal); lower and upper 95% CI values are presented in parentheses; all P-values * ** **** are Bonferroni corrected. P < 0.05; P < 0.01; ***P < 0.001; P < 0.0001; (all Bonferroni corrected). Figure 5. e t Th wo-tier model of cricket personality before the SSRI treatment (a), and aer t ft he SSRI treatment (b). ○ = rapid developmental line (n = 41), ◊ = slow developmental line (n = 37) and □ = control developmental line (n = 28). negatively associated with stress reactivity (i.e. time to resume activity) (r = −0.66, P < 0.0001, Table 2). Aer t ft he SSRI treatment, the relationships between the crickets’ body mass, coping styles and MMR disappeared, while the associations were still apparent between body mass and stress reactivity (r = −0.63, P < 0.0001) and body mass and RMR (r = 0.40, P < 0.001). RMR before the SSRI treatment was positively related with MMR (r = 0.40, P < 0.001), coping styles (r = 0.69, P < 0.0001) and negatively with stress reactivity (r = −0.74, P < 0.0001). Aer ft the SSRI treatment, the relationships between RMR, coping style and stress reactivity were weaker (r = 0.33, P < 0.01; r = −0.32, P < 0.05, respectively, Table 2). Before the SSRI treatment, MMR was also associated with coping styles (r = 0.37, P < 0.01) and stress reactivity (r = −0.40, P < 0.001), while aer t ft he treatment, MMR was not correlated with behaviors or RMR (Table 2). Prior to the SSRI treatment, a strong negative correlation was observed between the coping style and stress reactivity (r = −0.60, P < 0.0001), suggesting a behavioral syndrome (Fig. 5a). Aer t ft he treatment with antidepressants, however, the syndrome disappeared (Table  2, Fig. 5b). Discussion Evidence shows that slow development is associated with upregulation of stress-related genes . In humans, a predisposition toward experiencing negative emotions and low inhibitory control are linked to sickness and many 54–56 psychiatric conditions . It is likely that slower development in crickets confers higher levels of anxiety/neu- roticism along the stress reactivity axis since the SSRI treatment increased the time to resume movements (i.e. decreased anxiety) under the stressful conditions of the insect chamber. The relaxed stress reactivity under the SSRI treatment was associated with significant decreases in the RMR and MMR of slowly developing crickets (Fig. 6). The SSRI treatment presumably increased the concentration of 5-HT in the brain, while 5-HT could have lowered body temperature by decreasing the metabolic rate (Fig. 6). 5-HT activates glycolysis, forming lactic acid. 57,58 Excess lactic acid tends to decrease efficient energy production by interfering with mitochondrial respiration . On the other hand, the SSRI treatment did not ae ff ct stress reactivity of rapidly developing crickets. Moreover, the treatment that brings higher concentration of brain 5-HT increased RMR, MMR and exploration time of rapidly developing, oen b ft old crickets, while it decreased exploration time of slowly developing, oen s ft hy crickets. High concentration of brain 5-HT may cause a serotonin syndrome that oen in ft cludes high body temperature, indicat- ing a higher metabolic rate . Richmond et al. found that aquatic invertebrates are highly sensitive to pharmaceutical pollution in their environment: dipteran midges emerged significantly earlier and the total biomass of emerged adult insects was SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 6 www.nature.com/scientificreports/ Figure 6. A flow chart showing trait associations and the effects of SSRI treatment in slowly and rapidly developing crickets. Slow development is associated with higher (resting and maximum) metabolic rates (MR), stress reactivity (SR), exploration (Expl) and shy personality. Rapid development is associated with lower (resting and maximum) metabolic rates, stress reactivity, exploration and bold personality. By presumably raising 5-HT levels in the brain, SSRI decreases key trait values (MR, SR, Expl) in slowly developing crickets and raises these values in rapidly developing crickets. This results in the “normal” (intermediate) personality phenotype in slowly and rapidly developing crickets that is typical of the control line. greater in fluoxetine- and citalopram-exposed streams when compared with controls. This suggests that SSRI presence in the environment induces selection for faster development and more intense metabolism. This could mean that (i) rapidly and slowly developing crickets have different 5-HT tolerance, 5-HT sensitivity in brain receptors or different 5-HT concentration in the brain; (ii) 5-HT ae ff cts coping styles rather than responses along the stress reactivity axis and/or (iii) serotonergic metabolism is linked to growth. Under natural conditions, pro- active (active) coping may have higher survival value especially under high density of conspecifics/stronger com- petition than a reactive (passive) coping style . On an evolutionary scale, however, survival values for different behavioral types generally depend on environmental variability and predictability . Prior research has revealed differences in sensitivity to 5-HT levels in D. melanogaster . Ns3 mutant flies (which reach less than 60% of normal size and have fewer and smaller cells, but exhibit normal body proportions) and wild-type flies were treated with 5-HTP (5-hydroxy-L-tryptamine), the immediate precursor of 5-HT. This treatment resulted in a developmental delay in wild-type animals, reminiscent of the ns3 mutant phenotype, while ns3 mutants were more sensitive to 5-HTP than their wild-type counterparts. This suggests that ns3 mutants have already increased serotonergic signaling, as is likely in rapidly developing crickets. Using selective 5-HT1A and 5-HT1B agonists, de Boer and Koolhaas found a manyfold die ff rence between high- and low-aggression male rats. The high-aggression rats had a far more sensitive 5-HT1A auto-receptor and 5-HT1B receptor-mediated inhibition of the serotonergic neurons than low-aggression rats. However, the mechanisms that underlie these effects remain elusive. One reason is that the spatial distribution of serotonergic receptors remains poorly located and characterized in both the vertebrate and insect brain. As pointed out by de Boer and Koolhaas , it is crucial to know whether the receptors are located postsynaptically or presynaptically. At the presynaptic level, 5-HT agonists reduce 5-HT signaling, whereas at the level of the postsynaptic receptors, they mimic the effects of enhanced 5-HT signaling. In future research on 5-HT1A and 5-HT1B, agonists need to be tested in slowly and rapidly developing crickets to reveal any differences in the sensitivity of their 5-HT receptors and its relation to the development of personality along the stress reactivity and coping styles axes. Interestingly, the SSRI treatment decreased stress reactivity (i.e. increased the time to resume activity) only in slowly developing crickets while it ae ff cted exploration time and coping styles of both slowly and rapidly devel - oping individuals, but in the opposite way. Koolhaas et al. suggested a mechanism to explain the opposing effect of 5-HT on coping styles of individuals having to different personalities. Evidence shows that activation of the 5-HT1A auto-receptor leads (i) to a reduction of the preferred behavioral response in a forced swim test in rats and (ii) an increase in an alternative behavioral response rather than to a general reduction in anxiety . It is known that 5-HT1A receptor agonists in the reactively coping males show a reduction of immobility and an increase in escape attempts, whereas the proactively coping animal shows a decrease in escape and an increase in 64 29 immobility behavior in the test . Koolhaas et al . concluded that the effect of 5-HT may be related to the coping styles axis rather than to the emotionality axis. This idea is supported by experimental evidence obtained in the current study, where we show that 5-HT may be involved in regulation of behavior not only along a stress reactiv- ity gradient but also along a coping styles axis. These findings indicate that personality types mediate the strength and direction of 5-HT treatment, which should be taken into account in research on behavior and the underlying metabolism, the influence of pharmacological compounds on behavioral responses and developing new types of antidepressants. The proximate task of future research would be to reveal the types of 5-HT receptors that are involved in cricket behavioral responses and whether different receptor subtypes are engaged in stress-induced redistribution of immune reactions during behavioral responses. It is important to explain the link between developmental speed and body size found in this study, since it implies greater complexity in the structure and development of personality than expected before . We found that slowly developing crickets were significantly larger than rapidly developing crickets, which is in line with earlier SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 7 www.nature.com/scientificreports/ data on this species . Although faster developmental time in rapidly developing individuals may explain their smaller body size, rapidly developing, oen b ft old crickets are less anxious and can allocate more time to feeding, thus compensating for shorter time frames left for growth. The fact that the SSRI treatment ae ff cted slowly and rapidly developing crickets in the opposite way may indicate that rapidly developing crickets have higher baseline 5-HT in their brain than slowly developing crickets. In D. melanogaster, neuromodulatory actions of 5-HT were shown to depress feeding, while decreased neu- ronal 5-HT levels increased appetite . In another dipteran species, the flesh fly Neobellieria bullata , 5-HT injec- tion in the hemolymph decreased feeding . 5-HT and its effects on feeding are necessary to induce behavioral gregarization in migratory desert locusts . The decrease in body size is among the most obvious effects of phase 67,68 transition in morphological appearance from solitary to gregarious forms of locusts . Consistent with these findings, reduction in the brain 5-HT signaling by injection of a serotonergic-specific neurotoxin induces hyper - phagia and increased growth in rats . This may explain the larger body size of the slowly developing, oen s ft hy crickets. The observation that D. melanogaster ns3 mutants consume the same amount of food as control flies during early larval stages also suggests that serotonergic signaling appears to coordinate feeding behavior with 62 62 growth at the organismal level . Importantly, Kaplan et al. found close apposition of the axonal projections of serotonergic neurons to the insulin-producing cells (IPCs) in the brain of D. melanogaster, suggesting that these cells communicate to coordinate growth with 5-HT-modulated behaviors and environment. They showed that in D. melanogaster, insulin-like peptides both regulate blood sugar and act as growth factors, where serotonergic 62,70 neurons control the adult body size by ae ff cting insulin-like peptide secretion . This calls for studies on possi- ble similar synergies in cricket brains and bodies to define the functional patterns of the serotonergic synapse and subsequent effects on the phenotype. This may help to understand how and why 5-HT regulates somatic growth and the development of personality. For example, behaviors and sensitivity of 5-HT receptors might depend not only on the concentration of 5-HT but also on the proximity of the IPCs and serotonergic neurons in the brain. In conclusion, we show that personality-indicating consistent behavioral responses in crickets are tightly linked to developmental speed and body size. Antidepressants ae ff ct exploration and stress reactivity in different ways: we observed a higher SSRI impact on the coping styles axis than on the stress reactivity axis. Moreover, we found the opposite type of reactions to SSRI treatment in slowly and rapidly developing crickets along the coping styles axis. On the other hand, there were no simple one-way decreases or increases in exploration time at the two ends of the bimodal distribution of developmental speed / behavioral responses. This study shows that the SSRI treatment significantly decreased the variability of behaviors along the coping styles and stress reactivity axes, thus suggesting that antidepressants reduce variation in personality-related behaviors. It is also important to explain the link between developmental speed and body size found in this study, since it suggests greater com- plexity in the structure and development of personality than expected before . Overall, our results may have lim- ited explanatory power with regard to the magnitude of the SSRI effects on behaviors. This is because we have not studied SSRI activity in a dose-dependent manner, nor do we have data on possible side effects of SSRI treatment on cricket behavior(s). However, the high repeatability values that we recorded on the observed behaviors along the coping styles and the stress reactivity axes have the potential to show the diversity of directions that SSRI effects can take, and how these vary based on developmental speed and personality of the animals. Taken together, the present findings call for more research in brain neurophysiology and receptor distribution, a better definition of synaptic function, involvement of more personality axes and personality types (not only shy and bold ones), and an improved understanding of the brain concentrations of 5-HT and dose-dependent 5-HT effects in animals. Finally, this study shows that slowly and rapidly developing crickets markedly differ in behav- ioral traits and metabolism. Slowly developing individuals were larger, more explorative, less bold and they had higher RMR and lower MMR compared to rapidly developing crickets before the SSRI treatment. Overall, these findings do not support recent ideas concerning the pace-of-life syndrome (POLS) hypothesis , which predict higher exploration and metabolic rates in rapidly developing bold animals. Methods Insects and selection lines. We studied crickets from a laboratory stock originating from a wild popula- th tion (Davis, California, USA). The animals were selected for developmental speed; crickets of the 5 generation were investigated in this study. To start the selection lines, a total of 903 juvenile wild crickets were used. Until the crickets reached maturity, they were kept individually in plastic containers (280 mm × 115 mm × 105 mm) with a ventilation hole (60 mm dia.) covered with a plastic net. They were kept in a 12:12 h light–dark cycle and 26 ± 1 °C, with food and water ad libitum before and during experimental trials. The animals were selected for rapid and slow development, and we also maintained control groups. For rapid and slow development lines, mated males and females were selected according to their maturation age, with maximum 16.67% of the most rapidly or slowly maturing individuals being used, while in the control line the matings were randomized over the whole maturation age range . The first three generations were reared at the University of Eastern Finland, while the rest of the rearing of the lines occurred in Estonia. After five generations of selection for developmental speed, the development time (the average maturation time ± SD) for rapidly developing individuals was 98.03 ± 11.72 days (n = 238 crickets), 120.28 ± 25.12 days (n = 284) for the control individuals and 138.41 ± 23.05 days (n = 269) for slowly developing crickets. The mean longevity after eclosion for all three groups was 88.50 ± 23.54 days (one-way ANOVA, F = 25.79, (2, 788) P < 0.001). Longevity of males and females of all lines did not die ff r signic fi antly (one-way ANOVA, F = 1.0, (1, 789) P = 0.318). Lifespan aer ft eclosion of rapidly developing crickets (83.95 ± 23.52 days) and the control individuals (84.93 ± 22.06 days) did not differ significantly (Tukey’s test, P > 0.05), while lifespan after eclosion of slowly developing crickets (96.61 ± 25.05 days) significantly differed from the mean longevities of rapidly developing and control crickets (Tukey’s tests, both P < 0.001). SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 8 www.nature.com/scientificreports/ Trials, measurements and antidepressants. We tested whether escitalopram, a selective serotonin reuptake inhibitor (SSRI) (Brand name: CIPRALEX, CFL Pharmaceuticals Ltd) has any effect on cricket behavior and metabolism. Escitalopram is the (S)-stereoisomer (Left-enantiomer) of the citalopram prescribed frequently 60,71–74 to humans, low concentrations of which have been detected in natural habitats . Before the SSRI was added to the food of crickets, coping style (exploration) tests (day 1), and stress reactivity tests, measurements of RMR and MMR (day 2) were carried out for each of the F5 generation crickets (slow, n = 37; rapid, n = 41; control, n = 28) 15 days aer m ft aturation at the imago stage. To assess repeatability, we carried out coping style tests, stress reactivity tests and measurements of RMR and MMR one more time on days 3 and 4. As soon as these trials were done, on day 5 we started adding the SSRI to the food of the crickets (0.5 mg/kg per day added to Tetra Fish Food Flakes: Tetra, Blacksburg, VA) for twelve days (days 5–16). On days 13 and 14, the behavioral tests (coping styles and stress reactivity) and metabolism measurements were done again using the same animals as used during the previous trials. To assess repeatability of the behavioral variables, we repeated the tests on days 15 and 16. e Th insects were weighed before each measurement of metabolic rate using a Kern analytical balance (Kern & Sohn GmbH, Balingen, Germany). RMR and MMR values were corrected for body mass . Behavioral trials were conducted under constant temperature (24 ± 1.5 °C) and soundproof conditions. We used dim red light (25 W red incandescent bulb) since Gryllus spp. cannot see long (red) wavelengths properly , which allowed us to observe the nocturnal insects without disturbing them. The crickets were provided with drinking water before the onset of the trials, while food was removed 3 h prior the beginning of experimental trials. All behavioral trials were performed between 17:00 and 24:00 h. e Th observers of insects’ behavior and metabolic responses were blind to all predictions and treatment . Tests under familiar non-stressful environment: coping styles (exploration). The shelter part and the food area were separated by a corridor within individual containers. e Th crickets were trained for a week to run through a corridor (3 cm wide and 10 cm long) to reach the food aer a 3 ft h fast. During the fast, water was provided ad libitum. We used a one-cent Euro coin as a novel object before the crickets received the SSRI, while a one-dime USD coin was used during the SSRI treatment. The novel objects were placed on the floor of the cor - ridor 5 cm away from the entrance. We randomly chose the side of the corridor in which to place a novel object, and we always placed it on the opposite side during the second trial. As soon as the door between the shelter and the corridor was opened and the animal started to move through the corridor, we recorded the time it took to reach the food area. Since the crickets did not make any stops in the corridor when a novel object was not present, the time spent in the corridor reflects exploration activity / coping style of each individual – the extent of paying attention to changes in the environment. Tests under unfamiliar environments: stress reactivity. Next day after the coping style test, we handled each cricket for 1 min and then gave it an opportunity to escape into a conical plastic Eppendorf test tube (volume 5 ml) kept in dark conditions, which was used as an insect chamber during the measurement of RMR. This insect chamber represented a novel and stressful environment to the crickets. Resuming activity thus reflected each individual’s level of anxiety in this stressful environment. We assumed that the higher the anxiety, the sooner the individual will resume active body movements to avoid the stressful conditions of the windy and 26,27 noisy insect chamber . The latency to resume activity indicates the duration of freezing or immobile state, a 77,78 widespread anti-predator response used by many taxa . The insect chamber was connected to the respirometer by means of rubber tubing . Upon escape to the insect chamber, each of the crickets immediately became com- pletely immobile, as if hiding in a burrow or a crevice. We waited until all crickets resumed active movements, recorded the time between entering the tube and resuming activity, and used it as a measure of stress reactivity. Measuring RMR and MMR. The insects remained in their metabolic chambers until their CO emissions (and therefore their RMRs) reached the lowest stable rates. We also measured MMR at times when the crickets’ movements reached the highest rates while they performed struggling movements. As soon as the measure- ments were over, we returned the crickets back to their plastic housing-boxes. The LI-7000 die ff rential CO /H O analyzer (LiCor, Lincoln, NE, USA) was calibrated by means of calibra- 2 2 12,13,26 tion gases with gas injection . While measuring CO emissions, the insect chamber was perfused with dry −1 CO -free air at a flow rate of 60 ml min . The respirometric device was combined with an infrared optical system 12,13 using IR-emitting and IR-sensor diodes . IR-diodes made it possible to record CO production and to follow the movement intensity of each cricket. Statistical analyses. We used repeated measures ANOVAs (rmANOVAs) to assess changes in body mass, coping style, stress reactivity, resting metabolic rate and maximum metabolic rate aer a ft ntidepressant treatment (SSRI, within-subjects factor) in the same subjects with developmental line and sex used as between-subjects factors. We reported only the main effects when no significant interactions between between-subjects factors and/or interactions between within-subjects factor (antidepressant treatment) and between-subjects factors were found; otherwise the simple main effects using post hoc Bonferroni test were also reported. In order to account for age and testing apparatus effects, we performed rmANOVA with untreated (no antidepressant was given) groups of crickets (slow, n = 20; rapid, n = 20; control, n = 20). Analyses were performed in Statistica 8.0 for Windows (StatSoft Inc., Tulsa, OK, USA). To assess repeatability, we calculated the intraclass correlation coeffi- cient using a two-way random effects model for single measures . Pearson correlations were used to assess asso- ciations between different traits. Lower and upper 95% CI limits of correlation coefficient were calculated using non-stratified bootstrapping sampling with 2000 iterations. SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 9 www.nature.com/scientificreports/ References 1. Kain, J. S., Stokes, C. & de Bivort, B. L. Phototactic personality in fruit flies and its suppression by serotonin and white. Proc. Natl. Acad. Sci. USA 109, 19834–19839 (2012). 2. Buchanan, S. M., Kain, J. S. & de Bivort, B. L. Neuronal control of locomotor handedness in Drosophila. Proc. Natl. Acad. Sci. USA 112, 6700–6705 (2015). 3. Verbeek, M. E. M., Drent, P. J. & Wiepkema, P. R. Consistent individual differences in early exploratory behavior of male great tits. Anim. Behav. 48, 1113–1121 (1994). 4. Drent, P. J., van Oers, K. & van Noordwijk, A. J. Realized heritability of personalities in the great tit (Parus major). Proc. Biol. Sci. 270, 45–51 (2003). 5. Wolf, M., van Doorn, G. S., Leimar, O. & Weissing, F. J. Life-history trade-offs favour the evolution of animal personalities. Nature 447, 581–584 (2007). 6. Wolf, M., van Doorn, G. S. & Weissing, F. J. Evolutionary emergence of responsive and unresponsive personalities. Proc. Natl. Acad. Sci. USA 105, 15825–15830 (2008). 7. Kortet, R., Vainikka, A., Janhunen, M., Piironen, J. & Hyvärinen, P. Behavioral variation shows heritability in juvenile brown trout Salmo trutta. Behav. Ecol. Sociobiol. 6, 927–934 (2014). 8. Wolf, M. & Weissing, F. J. Animal personalities: consequences for ecology and evolution. Trends Ecol. Evol. 27, 452–461 (2012). 9. Smith, B. R. & Blumstein, D. T. Behavioral types as predictors of survival in Trinidadian guppies (Poecilia reticulate). Behav. Ecol. 21, 919–926 (2010). 10. Kortet, R., Hedrick, A. V. & Vainikka, A. Parasitism, predation and the evolution of animal personalities. Ecol Lett. 13, 1449–1458 (2010). 11. Dall, S. R. X., Bell, A. M., Bolnick, D. I. & Ratnieks, F. L. W. An evolution of ecology of individual differences. Ecol. Lett. 15, 1189–1198 (2012). 12. Krams, I. et al. Predation selects for low resting metabolic rate and consistent individual differences in anti-predator behavior in a beetle. Acta Ethol. 16, 163–172 (2013). 13. Krams, I. et al. Predation promotes survival of beetles with lower resting metabolic rates. Entomol Exp. Appl. 148, 94–103 (2013). 14. MacArthur, R. H. & Wilson, E. O. e Th Th eory of Island Biogeography (Princeton University Press, 1967). 15. Pianka, E. R. On r- and K-selection. Am. Nat. 104, 592–597 (1970). 16. Réale, D. et al. Personality and the emergence of the pace-of-life syndrome concept at the population level. Philos. Trans. R. Soc. London B Biol. Sci. 365, 4051–4063 (2010). 17. Biro, P. A. & Stamps, J. A. Are animal personality traits linked to life-history productivity? Trends Ecol. Evol. 23, 361–368 (2008). 18. Sih, A., Bell, A. M., Johnson, J. C. & Ziemba, R. E. Behavioral syndromes: an integrative overview. Q. Rev. Biol. 79, 241–277 (2004). 19. Réale, D., Reader, S. M., Sol, D., McDougall, P. T. & Dingemanse, N. J. Integrating animal temperament within ecology and evolution. Biol. Rev. Camb. Philos. Soc. 82, 291–318 (2007). 20. Steimer, T., La Fleur, S. & Schulz, P. E. Neuroendocrine correlates of emotional reactivity and coping in male rats from the Roman high (RHA/Verh)- and low (RLA/Verh)-avoidance lines. Behav. Genet. 27, 503–512 (1997). 21. Benus, R. F., Koolhaas, J. M. & van Oortmerssen, G. A. Aggression and adaptation to the light–dark cycle: role of intrinsic and extrinsic control. Physiol. Behav. 43, 131–137 (1988). 22. Bolhuis, J. E., Schouten, W. G., de Leeuw, J. A., Schrama, J. W. & Wiegant, V. M. Individual coping characteristics, rearing conditions and behavioural flexibility in pigs. Behav. Brain. Res. 152, 351–360 (2004). 23. West, G. B., Brown, J. H. & Enquist, B. J. A general model for ontogenetic growth. Nature 413, 628–631 (2001). 24. Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Metabolic theory of ecology. Ecology 85, 1771–1789 (2004). 25. Speakman, J. R. Doubly Labelled Water: Theory and Practice (Chapman & Hall, 1997). 26. Lighton, J. R. B. Measuring Metabolic Rate: a Manual for Scientists (Oxford University Press, 2008), 27. Careau, V., Thomas, D., Humphries, M. M. & Réale, D. Energy metabolism and animal personality. Oikos 117, 641–653 (2008). 28. Krams, I. A. et al. Metabolic rate associates with, but does not generate covariation between, behaviours in western-trilling crickets, Gryllus integer. Proc. Biol. Sci. 284, e20162481 (2017). 29. Koolhaas, J. M., de Boer, S. F., Coppens, C. M. & Buwalda, B. Neuroendocrinology of coping styles: Towards understanding the biology of individual variation. Front. Neuroendocrinol. 31, 307–321 (2010). 30. Seltmann, M. W. et al. Stress responsiveness, age and body condition interactively ae ff ct flight initiation distance in breeding female eiders. Anim. Behav. 84, 889–896 (2012). 31. Koolhaas., J. M. et al. Coping styles in animals: current status in behavior and stress-physiology. Neurosci. Biobehav. Rev. 23, 925–935 (1999). 32. Gelder, M., Mayou, R. & Geddes, J. Psychiatry. (Oxford University Press, 2005). 33. Kravitz, E. A. & Huber, R. Aggression in invertebrates. Curr. Opin. Neurobiol. 13, 736–743 (2003). 34. Anstey, M. L., Rogers, S. M., Ott, S. R., Burrows, M. & Simpson, S. J. Serotonin mediates behavioral gregarization underlying swarm formation in desert locusts. Science 323, 627–630 (2009). 35. Donovan, M. H. & Tecott, L. H. Serotonin and the regulation of mammalian energy balance. Front. Neurosci. 7, e36 (2013). 36. De Boer, S. F. & Newman-Tancredi, A. Anti-aggressive effects of the selective high-efficacy ‘biased’ 5-HT1A receptor agonists F15599 and F13714 in male WTG rats. Psychopharmacology 233, 937–947 (2016). 37. Dolan, R. J., Morris, J. S. & de Gelder, B. Crossmodal binding of fear in voice and face. Proc. Natl. Acad. Sci. USA 98, 10006–10010 (2001). 38. Murphy, D. L. & Lesch, K. P. Targeting the murine serotonin transporter: insights into human neurobiology. Nat. Rev. Neurosci. 9, 85–96 (2008). 39. Caramaschi, C., De Boer, S. F. & Koolhaas, J. M. Differential role of the 5-HT1A receptor in aggressive and non-aggressive mice: an across-strain comparison. Physiol. Behav. 90, 590–601 (2007). 40. Steimer, T. & Driscoll, P. Divergent stress responses and coping styles in psychogenetically selected Roman high-(RHA) and low- (RLA) avoidance rats: behavioural, neuroendocrine and developmental aspects. Stress 6, 87–100 (2003). 41. Niemelä, P. T., Vainikka, A., Lahdenperä, S. & Kortet, R. Nymphal density, behavioural development, and life history in a field cricket. Behav. Ecol. Sociobiol. 66, 645–652 (2012). 42. Niemelä, P. T., Vainikka, A., Hedrick, A. V. & Kortet, R. Integrating behaviour with life history: boldness of the field cricket, Gryllus integer, during ontogeny. Funct. Ecol. 26, 450–456 (2012). 43. Niemelä, P. T., Dingemanse, N., Alioravainen, N., Vainikka, A. & Kortet, R. Personality pace-of-life hypothesis: testing genetic associations among personality and life-history. Behav. Ecol. 24, 935–941 (2013). 44. Rothwell, N. J. & Stock, M. J. Effect of diet and fenfluramine on thermogenesis in the rat: possible involvement of serotonergic mechanisms. Int. J. Obes. 11, 319–324 (1987). 45. Le Feuvre, R. A., Aisenthal, L. & Rothwell, N. J. Involvement of corticotrophin releasing factor (CRF) in the thermogenic and anorexic actions of serotonin (5-HT) and related compounds. Brain Res. 555, 245–250 (1991). 46. Laporta, J. et al. Increasing serotonin concentrations alter calcium and energy metabolism in dairy cows. J. Endocrinol. 226, 43–55 (2015). SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 10 www.nature.com/scientificreports/ 47. Lin, M. T., Pang, I. H., Chern, S. I. & Chia, W. Y. Changes in serotonin contents in brain ae ff ct metabolic heat production of rabbits in cold. Am. J. Physiol. 235, R41–47 (1978). 48. Nonogaki, K., Abdallah, L., Goulding, E. H., Bonasera, S. J. & Tecott, L. H. Hyperactivity and reduced energy cost of physical activity in serotonin 5-HT2C receptor mutant mice. Diabetes 52, 315–320 (2003). 49. Weng, R. et al. Metabolomics Approach Reveals Integrated Metabolic Network Associated with Serotonin Deficiency. Sci. Rep. 5, e11864 (2015). 50. Klein, D. F. Delineation of two drug-responsive anxiety syndromes. Psychopharmacologia 5, 397–408 (1964). 51. Graeff, F. G. & Del-Ben, C. M. Neurobiology of panic disorder: from animal models to brain neuroimaging. Neurosci. Biobehav. Rev. 32, 1326–1335 (2008). 52. De Boer, S. F. & Koolhaas, J. M. 5-HT1A and 5-HT1B receptor agonists and aggression: a pharmacological challenge of the serotonin deficiency hypothesis. Eur. J. Pharmacol. 526, 125–139 (2005). 53. Gutiérrez-Adán, A. et al. Effect of speed of development on mRNA expression pattern in early bovine embryos cultured in vivo or in vitro. Mol. Reprod. Dev. 68, 441–448 (2004). 54. Chung, J. Y. Y. & Evans, M. A. Shyness and symptoms of illness in young children. Can. J. Behav. Sci. 32, 49–57 (2000). 55. Rettew, D. C. & McKee, L. Temperament and its role in developmental psychopathology. Harv. Rev. Psychiatry 13, 14–27 (2005). 56. Rantala, M. J., Luoto, S., Krams, I. & Karlsson, H. Depression subtyping based on evolutionary psychiatry: Proximate mechanisms and ultimate functions. Brain Behav. Immun. 69, 603–617 (2018). 57. Koren-Schwartzer, N., Chen-Zion, M., Ben-Porat, H. & Beitner, R. Serotonin-induced decrease in brain ATP, stimulation of brain anaerobic glycolysis and elevation of plasma hemoglobin; the protective action of calmodulin antagonists. Gen. Pharmacol. 25, 1257–1262 (1994). 58. Gardner, A. & Boles, R. G. Beyond the serotonin hypothesis: mitochondria, inflammation and neurodegeneration in major depression and ae ff ctive spectrum disorders. Prog. Neuropsychopharmacol. Biol. Psychiatry 35, 730–743 (2011). 59. Volpi-Abadie, J., Kaye, A. M. & Kaye, A. D. Serotonin syndrome. Ochsner J. 13, 533–540 (2013). 60. Richmond, E. K., Rosi-Marshall, E. J., Lee, S. S., Thompson, R. M. & Grace, M. R. Antidepressants in stream ecosystems: influence of selective serotonin reuptake inhibitors (SSRIs) on algal production and insect emergence. Freshw. Sci. 35, 845–855 (2016). 61. Niemelä, P. T., Vainikka, A., Forsman, J. T., Loukola, O. J. & Kortet, R. How does variation in the environment and individual cognition explain the existence of consistent behavioural differences? Ecol. Evol. 3, 457–464 (2013). 62. Kaplan, D. D., Zimmermann, G., Suyama, K., Meyer, T. & Scott, M. P. A nucleostemin family GTPase, NS3, acts in serotonergic neurons to regulate insulin signaling and control body size. Genes Dev. 22, 1877–1893 (2008). 63. Veenema, A. H. et al. Differences in the effects of 5-HT(1A) receptor agonists on forced swimming behavior and brain 5-HT metabolism between low and high aggressive mice. Psychopharmacology 178, 151–160 (2005). 64. Cryan, J. F., Redmond, A. M., Kelly, J. P. & Leonard, B. E. The effects of the 5-HT1A agonist flesinoxan, in three paradigms for assessing antidepressant potential in the rat. Eur. Neuropsychopharmacol. 7, 109–114 (1997). 65. Neckameyer, W. S. A trophic role for serotonin in the development of a simple feeding circuit. Dev. Neurosci. 32, 217–237 (2010). 66. Dacks, A. M., Nickel, T. & Mitchell, B. K. An examination of serotonin and feeding in the flesh fly Neobellieria bullata (Sarcophagidae: Diptera). J. Insect. Behav. 16, 1–21 (2003). 67. Pener, M. & Simpson, S. Locust phase polyphenism: an update. Adv. Insect. Physiol. 36, 1–272 (2009). 68. Uvarov, B. P. Grasshoppers and Locusts (Cambridge University Press, 1966). 69. Saller, C. F. & Stricker, E. M. Hyperphagia and increased growth in rats aer in ft traventricular injection of 5,7-dihydroxytryptamine. Science 192, 385–387 (1976). 70. Ruaud, A. F. & Thummel, C. S. Serotonin and insulin signaling team up to control growth in Drosophila. Genes Dev. 22, 1851–1855 (2008). 71. Schultz, M. M. & Furlong, E. T. Trace analysis of antidepressant pharmaceuticals and their select degradates in aquatic matrixes by LC/ESI/MS/MS. Anal. Chem. 80, 1756–1762 (2008). 72. Metcalfe, C. D. et al. Antidepressants and their metabolites in municipal wastewater, and downstream exposure in an urban watershed. Environ. Toxicol. Chem. 29, 79–89 (2010). 73. Cardoso, O., Porcher, J.-M. & Sanchez, W. Factory-discharged pharmaceuticals could be a relevant source of aquatic environment contamination: review of evidence and need for knowledge. Chemosphere 115, 20–30 (2014). 74. Larsson, D. G. J. Pollution from drug manufacturing: review and perspectives. Philos. Trans. R. Soc. Lond. B Biol. Sci. 369, e20130571 (2014). 75. Briscoe, A. D. & Chittka, L. The evolution of color vision in insects. Ann. Rev. Entomol. 46, 471–510 (2001). 76. Burghardt, G. M. et al. Perspectives – minimizing observer bias in behavioral studies: a review and recommendations. Ethology 118, 511–517 (2012). 77. Chelini, M.-C., Willemart, R. H. & Hebets, E. A. Costs and benefits of freezing behavior in the harvestman Eumesosoma roeweri (Arachnida, Opiliones). Behav. Processes. 82, 153–159 (2009). 78. Krams, I. et al. High repeatability of anti-predator responses and resting metabolic rate in a beetle. J. Insect. Behav. 27, 57–66 (2014). 79. Mattila, A. L. K. & Hanski, I. Heritability of flight and resting metabolic rates in the Glanville fritillary butterfly. J. Evol. Biol. 27, 1733–1743 (2014). 80. McGraw, K. O. & Wong, S. P. Forming inferences about some intraclass correlation coefficients. Psychol. Methods 1, 30–46 (1996). Acknowledgements We thank Ilkka Sirkka, Sanni Hietaketo, Leena Koponen, Leena Pääkkönen, Julia Keronen, Sirpa Kaunisto, Anssi Vainikka, Inese Kivleniece, Jolanta Vrublevska and Ljubova Sivačova for assistance, and Raivo Mänd for comments on the manuscript. Author Contributions I.K., G.T., Ra.K. and T.K. conceived the study. I.K., A.K., R.K., Ra.K., T.K., R.K., M.M., M.J.R., J.K., P.J. and S.K. contributed to the design of the study and carried out the experiments. D.E., G.T. and I.K. analyzed the data. I.K., S.L. and G.T. wrote the manuscript. All authors participated in results interpretation and they approved the final manuscript. Additional Information Competing Interests: The authors declare no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 11 www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2018 SCIENTIFIC REPo R TS | (2018) 8:8599 | DOI:10.1038/s41598-018-26722-9 12

Journal

Scientific ReportsSpringer Journals

Published: Jun 5, 2018

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off