- Open Access
Effects of predation stress and food ration on perch gut microbiota
Microbiome volume 6, Article number: 28 (2018)
Gut microbiota provide functions of importance to influence hosts’ food digestion, metabolism, and protection against pathogens. Factors that affect the composition and functions of gut microbial communities are well studied in humans and other animals; however, we have limited knowledge of how natural food web factors such as stress from predators and food resource rations could affect hosts’ gut microbiota and how it interacts with host sex. In this study, we designed a two-factorial experiment exposing perch (Perca fluviatilis) to a predator (pike, Esox lucius), and different food ratios, to examine the compositional and functional changes of perch gut microbiota based on 16S rRNA amplicon sequencing. We also investigated if those changes are host sex dependent.
We showed that overall gut microbiota composition among individual perch significantly responded to food ration and predator presence. We found that species richness decreased with predator presence, and we identified 23 taxa from a diverse set of phyla that were over-represented when a predator was present. For example, Fusobacteria increased both at the lowest food ration and at predation stress conditions, suggesting that Fusobacteria are favored by stressful situations for the host. In concordance, both food ration and predation stress seemed to influence the metabolic repertoire of the gut microbiota, such as biosynthesis of other secondary metabolites, metabolism of cofactors, and vitamins. In addition, the identified interaction between food ration and sex emphasizes sex-specific responses to diet quantity in gut microbiota.
Collectively, our findings emphasize an alternative state in gut microbiota with responses to changes in natural food webs depending on host sex. The obtained knowledge from this study provided us with an important perspective on gut microbiota in a food web context.
Animal hosts provide associated microorganisms with suitable ecological niches in their intestines . The intestine of most animals such as human and fish develops from an initial sterile environment, followed by a subsequent microbial colonization leading to a matured intestine inhabited by a diverse microbial community . These microbiota help the hosts digest food, protect against pathogens, and influence the host’s metabolisms . Consequently, the mutual benefits between host and their gut microbiota may contribute to the host fitness through metabolites .
Gut microbiota have been shown to be affected by host genetics . For example, Bolnick et al.  have shown that the variation in gut microbiota in three-spine stickleback was affected by the major histocompatibility class II (MHC) genotypes. In addition, host sex, another genetic trait, has been linked to gut microbial composition [7, 8], which furthermore, can interact with other environmental factors in affecting gut microbial composition, e.g., sex-dependent effects of diet on gut microbiota . Besides host genetics, host diet choice is also an important determinant of gut microbiota composition. Bacteria differ in their substrate use; thus, niche specialization in gut microbiota causes changes in bacterial taxa as a consequence of diet choice of the host .
In natural animal populations, there is a relationship between gut microbial composition and their host trophic level in the food web . One important aspect of trophic level for gut microbiota is that it is associated with shifts in diet quality and quantity [11, 12]. Ley et al. , Sullam et al. , and Liu et al.  have shown that gut microbial composition could change along trophic levels resulting from diet change in mammals and fish. Furthermore, trophic interactions can influence the abundance of organisms along trophic levels through cascading effects, which are called trophic cascades , which in turn indirectly could influence the gut microbiota. Increasing predator density will decrease prey densities and indirectly increase the food resources availability for the prey, thus lowering intraspecific competition among the prey . On the other hand, when predator density is low, competition among the prey will increase, resulting in less food availability [17, 18]. Low food ration, i.e., reduced food intake, especially in the case of starvation and hibernation has been shown to affect gut microbiota [19,20,21]. For example, hibernating brown bear and squirrels show reduced gut microbiota diversity and reduced levels of certain phyla [21, 22]. Carey et al.  have found that the remaining microbiota phyla during hibernation mainly consist of taxa that can use host-derived substrates as a food source.
Besides affecting resource levels for its prey, predation is also an important food web factor that may cause stress for individual prey in nature. Chemical cues from predators have been shown to reduce activities in fish and subsequently induce morphological changes that could decrease predation risk [23,24,25]. Moreover, predation stress also influences the physiological status of prey, for example, hormones released from stress could mediate immunological and behavioral responses in vertebrates . It has been shown that nerve and immune system can play important roles in regulating gut microbiota communities [27,28,29]. O’Mahony et al.  showed that early life stress from maternal separation in rats could alter the gut microbiota in the offspring. Thus, the association between trophic level and gut microbiota may not only depend on diet shifts along trophic levels but also on the risk of predation, i.e., stress. However, this association of trophic levels and gut microbiota is less studied despite the prevalence of predation and its effects on prey.
In this study, we use Eurasian perch (Perca fluviatilis) as prey and Northern pike (Esox lucius) as predators to investigate the relative importance of predator presence and food ration to perch gut microbial communities. Previous studies have shown that predator cues  and food availability  can affect perch behavior and morphology. As stress has been shown to affect gut microbiota, we predict that (1) pike predation stress could alter perch gut microbial community composition and consequently affect their functions; (2) increasing the amount of food fed to perch will also influence gut microbial communities, as it will change competition within the microbial communities. Furthermore, previous studies of perch have shown that sex affects both composition and diversity of gut microbiota as well as interacts with diet in affecting perch gut microbiota [9, 33]. Similarly, (3) if microbes that respond to stress or food ration are also influenced by sex hormones, we might expect to see stress and food ration effects that differ between fish sexes.
We collected 1-year-old perch from Lake Mälaren (N59° 20′, E17° 52′) in Sweden in May 2013 using cast net. We also collected pike (341.6 ± 49.2 mm, 207.3 ± 76.7 g, mean fish length and weight ± SD) from Lake Messormen and Hersjön between May and July in 2013. We acclimated all perch and pike to lab conditions for 6 weeks before starting the experiment. During the acclimation, we fed perch with frozen chironomids (Imazo AB, Sweden) daily, and pike were fed with perch (from the same pool of perch for the experiment) two times a week.
To examine how predation stress, food ration, and fish sex affect gut microbiota of perch, four perch were put into one aquarium with or without the presence of pike. The 105 L aquaria (75 × 40 × 35 cm, length × width × height) were divided into two parts by a transparent plastic board to separate pike and perch. The plastic boards had holes drilled into them to allow for predator kairomones to freely pass to the other side of the aquarium containing perch. This setup allowed perch to be affected by predator cues both visually and by olfactory in each aquarium, but all aquaria were visually isolated from each other to avoid pikes influencing perch in predator-free treatments. All aquaria were kept at temperatures ranging between 19 and 20 °C with a thermostat heater in each aquarium under a photoperiod of 12-h light and 12-h dark. Each treatment with four perch was replicated six times resulting in 36 aquaria.
All perch were fed once a day with frozen chironomids at three quantity rations, 5, 10, and 15% of the average perch weight, in which 15% ration is close to the maximum food conversion at the specific size and temperature in perch . All pikes were hand-fed with one juvenile perch two times a week during the experiment. We ran the experiment for 10 weeks to observe perch growth and morphological changes for an accompanying paper from this study  as well as to allow the perch gut microbiome to adapt to the treatments. Due to occasional death, 91 perch remained. All perch were killed with an overdose of benzocaine (ethylene glycol monophenyl ether, Merck). We recorded final weight, length, and sex (41 females and 48 males, 2 undetermined) and calculated relative intestine length (intestine length/fish length) for all fish (Additional file 1: Table S1). The entire intestine including both intestine tissue and the gut content from each fish was immediately frozen and stored at − 80 °C until later analysis of bacterial composition. To assess bacterial community composition in the surrounding water, we filtered 50 ml water through 0.2 μm Supor 200 filters (Pall Corporation, Port Washington, NY, USA) from each aquarium at the end of the experiment and stored the filter at − 80 °C. In addition, we took samples from the chironomids to check diet-associated bacteria.
DNA extraction and bacterial 16S rRNA genes Illumina sequencing
The entire intestine from perch, water filter samples, and 0.25 g chironomids were processed to extract bacterial DNA using PowerSoil® DNA Isolation Kit (MO BIO Laboratories, Inc., Carlsbad, CA, USA) with a modification from the manufacture protocol in which we incubated the samples at 65 °C for 10 min after adding the C1 solution.
The variable region V4 of the 16S rRNA gene was amplified by using bacterial primers (515F and 806R) . Polymerase chain reaction was done with two steps . Triplicates of 20 μl reaction for each sample were done in the first step PCR with 515F (5’-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). Each reaction contained 10 μM of forward and reverse primers, 1× reaction buffer, 200 μM of dNTPs and 0.02 U/μl Q5 HF DNA polymerase, and 1 μl of DNA template. The reaction started with initial denaturation at 98 °C for 30 s, and then 30 cycles of denaturation at 98 °C for 10 s, annealing at 58 °C for 30 s, and extension at 72 °C for 30 s. It was finished with a final extension at 72 °C for 2 min. Triplicate PCRs for each sample were pooled and purified with Agencourt® AMPure® XP (Beckman Coulter). The purified sample was used as the template for the second step PCR to attach the Illumina handles and index primers. Triplicates of PCR product for each sample were prepared. Each reaction contained 1.25 μM of forward and reverse primers, 5 × reaction buffer, 2 mM of dNTPs and 2 U/μl Q5 HF DNA polymerase, and 1 μl of DNA template. Each reaction started with initial denaturation at 98 °C for 30 s and then continued with 20 cycles of denaturation at 98 °C for 10 s, annealing at 65 °C for 30 s and extension at 72 °C for 30 s. A final extension was done at 72 °C for 2 min. Each sample was purified with Agencourt® AMPure® XP (Beckman Coulter) and quantified with Quant-iT™ PicoGreen® dsDNAReagent Kit (Invitrogen). Equal amounts of DNA were mixed in one pool with a final concentration of 2.7 ng/μl. Samples were sent for Illumina Miseq sequencing at National Genomics Infrastructure by ScilifeLab, Uppsala, in Sweden.
Sequencing data analysis
The raw amplicon sequencing data was demultiplexed, and sequence pairs were assembled using a pipeline developed by Sinclair et al. . The pipeline further removed sequences with missing primers and unassigned base pairs. Resulting quality-filtered assembled reads were clustered into operational taxonomical units (OTUs) using UPARSE (cutoff of 3% sequence dissimilarity) . Taxonomy was assigned using CREST  and the ribosomal sequence database SilvaMod.
We used PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States, version 1.1.1)  to obtain the relative abundance of gene families (gene ontology categories or GO)  within individual perch. A closed reference OTUs were prepared at 97% level against the gg_13_5_otus.tar.gz from Greengenes using Macqiime (1.9.1 20150604) before using PICRUSt. The newly picked OTUs were then used as input for PICRUSt following the workflow suggested by the developers, including normalization by dividing each OTU by the known/predicted 16S copy number abundance, and then calculated the final metagenome functional predictions. The predicted functions were then categorized with KEGG pathways on level 2. Quality control steps for PICRUSt were also performed, which gave the percentage of successful reads that were mapped to Greengenes when using the closed reference OTU picking (Additional file 2: Table S2), and the calculations of the reference genome coverage for each fish sample presented as NSTI scores (Additional file 3: Table S3). While PICRUSt provide some first predictions on metabolic pathway direct methodologies such as shotgun metagenomics can be used to confirm the metabolic pathways predictions by PICRUSt.
The 97% OTU table was rarified to 8000 reads per sample before the statistical analysis. All statistical analyses were performed in R (version 3.2.2). Alpha diversities were calculated using observed richness (S. Obs) and Chao1 by package phyloseq (version 1.12.2). Faith’s phylogenetic diversity (PD) was estimated using the picante package (version 1.6-2).
We used generalized linear models (GLMs) with quasibinomial link function to test whether food ration, pike presence, and sex could affect the relative abundance of all the functional categories obtained from PICRUSt and each of the top 10 phyla, which represented more than 90% of the total relative abundance in the whole microbial community. For the phyla with significant treatment effect(s), we re-analyzed to microbiome shifts at the genus level, also using GLMs with quasibinomial link function.
To observe the overall pattern of microbial community composition across all the treatments, we used non-metric multidimensional scaling (NMDS) with Bray-Curtis (vegan 2.3-5), weighted (w), and unweighted (uw) UniFrac distance matrices calculated by phyloseq. Next, we used PERMANOVA to test the effects of all treatments on Bray-Curtis, weighted and unweighted UniFrac distance matrix with 10,000 permutations using vegan. We used the r2 value from PERMANOVA to estimate the relative effect size (% of variation explained) of our treatments (food ration, pike, sex, and the interactions food ration × pike, food ration × sex, and pike × sex) on the gut microbiota.
We identified the OTUs that were over- or under-represented in each treatment using EdgeR (version 3.10.5) . EdgeR was originally designed for differential expression analysis of RNA-seq expression profiles but can be applied to any technology that produces read counts for genomic features. In our analysis, we adapted EdgeR to test the OTUs that were significantly over- or under-represented in response to the three factors (pike presence, food ration, and sex). We analyzed over- and under-representative OTUs both for each factor alone and the interactions between two factors (food ration × pike, pike × sex, and food ration × sex). We then assigned each OTU to their minimum genus level and used Fisher’s exact test to check at the phyla level which phylum can be representative for the treatment based on the contingency table including number of “success OTUs” (number of significant representative OTU in one phylum from EdgeR analysis), “failure OTUs” (total number of OTU within a phylum subtracted with number of “success” OTUs), and number of “success” and “failure” OTUs in all phyla. Finally, we calculated the average relative abundance at both OTU and phyla levels.
Abundant phyla in perch gut microbiota communities
The relative abundance of the ten most abundant phyla varied among treatments (Fig. 1a, Additional file 4: Table S4). Tenericutes was the most abundant phylum, and it increased with food ration whereas the relative abundances of Fusobacteria and Proteobacteria seemed to decrease with increasing food ration. However, the effects of food ration on Tenericutes and Proteobacteria were non-significant (p > 0.05), but significantly affected Fusobacteria (F2, 78 = 5.420, p = 0.006). The relative abundance of Fusobacteria was also significantly affected by pike presence (F1, 78 = 6.114, p = 0.016). Pike presence also significantly affected the relative abundance of Proteobacteria (F1,78 = 4.833, p = 0.031). The relative abundance of Bacteroidetes was, however, significantly influenced by the interaction of food ration and sex (Fig. 1b, ANOVA, F2, 78 = 3.592, p = 0.032). The relative abundance of Cyanobacteria decreased significantly with pike presence (Fig. 1c, ANOVA, F1, 78 = 11.614, p = 0.001) and was marginally insignificantly influenced by the interaction between food ration and sex (ANOVA, F2, 78 = 2.730, p = 0.071).
At the genus level, the interaction of food ration and sex significantly influenced the relative abundance of Myroides (F2, 82 = 5.852, p = 0.004) and an unknown Flavobacteriaceae genus (F2, 82 = 4.122, p = 0.020) from the phylum Bacteriodetes (Additional file 5: Figure S1a). Pike presence also significantly affected the relative abundance of Prochiorococcus (F1, 88 = 11.98, p = 0.001) and Anabaenopsis (F1, 88 = 4.948, p = 0.028) from the phylum Cyanobacteria (Additional file 5: Figure S1b), and Cetobacterium (F1, 88 = 5.457, p = 0.021) and an unknown Fusobacteriaceae genus (F1, 88 = 15.99, p = 0.0001) (Additional file 5: Figure S1c).
Alpha diversity in perch gut microbiota
Alpha diversity of perch gut microbiota, such as Chao1 (Fig. 2a), PD, and S. Obs (Additional file 6: Figure S2, Additional file 7: Table S5), significantly decreased in the presence of pike (Additional file 7: Table S5). When adding food ration and sex in the model, their interaction had a significant effect on Chao1 (Fig. 2b, t = − 2.03, p = 0.046), showing a sex-dependent effect of increasing food ration. However, this interaction between sex and food ration could not be observed in PD and S. Obs. We also added perch relative intestine length as a covariate to our test, but we did not find significant effects from it on microbiota alpha diversity (p = 0.765).
Gut microbiota composition among perch individuals
Perch gut microbiota communities were distinctively separated from water and food samples in terms of NMDS using Bray-Curtis (Fig. 3a) distance matrix. Unweighted and weighted UniFrac distance matrix-based NMDS also showed similar patterns (Additional file 8: Figure S3). NMDS for only gut microbial communities showed that predation stress was the main factor to determine the microbial compositional variations (Fig. 3b). PERMANOVA on Bray-Curtis, unweighted UniFrac distances showed that both pike presence and food ration had significant effects on perch gut microbiota (Table 1). In contrast, sex, food ration × sex, pike × sex, and food ration × pike were all non-significant. PERMANOVA on the weighted UniFrac distances corroborated this since food ration significantly affected perch gut microbiota, while pike presence had a marginally insignificant effect (Table 1). However, sex, food ration × sex, pike × sex, and food ration × pike were all found non-significant.
There were large overlaps in OTUs of gut microbiota within predation treatments, but perch with pike presence had fewer unique OTUs compared to perch with no pike (Fig. 4a). The overlaps were smaller and the number of unique OTUs was greater among the perch with different food ration (Fig. 4b). By detailed profiling of the OTU dynamics in each treatment, we found 23 OTUs belonged to, for example, Cetobacterium and Fusobacteriaceae in Fusobacteria, Lactococcus and Clostridium in Firmicutes, that were significantly over-represented in the pike treatments (Additional file 9: Table S6). At phylum level, Firmicutes and Fusobacteria were tested to have more odds to respond to pike presence/absence than the other phyla (Firmicutes: odds ratio = 7.64, p < 0.001, Fusobacteria: odds ratio = 25.19, p = 0.004). When considering the interaction of factors, we also found that Firmicutes and Fusobacteria had significantly more odds to respond than other phyla in the 10% food ration and pike presence interactions (Firmicutes: odds ratio = 6.27, p < 0.001, Fusobacteria: odds ratio = 20.68, p = 0.005). However, the sum of all responding OTUs only made up a small proportion of the total relative abundance of all phyla (Additional file 10: Figure S4).
PICRUSt gave 44 predicted functional categories that represented 7 pathway maps in the KEGG level 2 functional modules. Each functional category showed variations in their average relative abundance (Additional file 11: Table S7). We did not find any treatment effects in an overall model including all 2-way interactions. However, when analyzing the effects of food ration and pike presence separately, we found that food ration and pike presence had significant effects on several functional categories, corresponding to 5 pathway maps (Additional file 11: Table S7). Food ration had a significant effect on metabolic pathways, involved in biosynthesis of other secondary metabolites, metabolism of cofactors and vitamins, and digestive system (Fig. 5a). Pike presence influenced the functional categories membrane transport, signaling molecules and interaction as well as environmental adaptation (Fig. 5b).
Perch intestine length
Besides the effect on microbial communities, we found that pike presence and food ration had a significant impact on perch relative intestine length (ANOVA: pike, F1, 85 = 16.18, p = 0.0001; food ration, F2, 85 = 11.09, p < 0.0001). Perch intestine length decreased in the presence of pike and increased with increased food ration (Fig. 6). Despite strong effects of treatments on both intestine length and microbial diversity, we found no significant correlations between perch relative intestine length and gut microbiota diversity (Chao1: t = 1.016, p = 0.313; PD: t = 0.825, p = 0.411; S. Obs: t = 0.608, p = 0.545).
Gut microbiota serve the host with crucial roles in immune function and metabolism. The coupling of gut microbiota and their host has been attributed to the production of vitamins and other metabolites by microbiota in the gut of aquatic vertebrates, which is similar to what has been found in terrestrial mammals [42, 43]. For example, enzymes from gut microbes are important sources besides the enzymes produced by the fish gut for food digestion . In this study, we tested the effects of predation stress, food ration, and host sex on gut microbiota in perch. We found that gut microbial diversity, as well as metabolic potential predicted from 16S rRNA, responded to predation stress and food ration. Furthermore, Bacteroidetes communities showed sex-dependent responses to food ration. These community responses to the treatments coincided with decreased intestine length in the presence of pike and increased intestine length with larger food ration. Hence, it can be suggested that predation stress and food ration have consequences for fish body condition through inducing changes in gut microbial communities, and these changes can be sex-dependent.
Like previous studies on human and laboratory animals, we have shown that stress can influence and change the gut microbiota community. Knowles et al.  found significantly lower fecal lactic acid bacterial levels when students were facing academic stress. Similarly, Bailey et al.  showed that social disruption stressor could impact the gut microbiota community in mice. In our perch experiment, stress responses were observed in the relative abundance of Proteobacteria and Fusobacteria. A well-studied response to predator presence is that prey will change habitat and diet to reduce the risk of being captured by the predator [47,48,49]. To eliminate the effect of diet type which has been shown to shape gut microbiota of multiple vertebrate species [50,51,52,53,54], we only fed perch with one type of food (chironomids). Hence, our experimental setup allowed us to show that direct physiological responses of perch to predation-stress could modify gut microbiota communities, and thus, in our case, diet shifts could be excluded. In an accompanying paper from this experiment, we found that perch reduced their foraging activity and space use in the presence of predators . Such responses in behavior emphasize that the perch actually experienced predation stress in our experiment. Furthermore, predation stress also changed the relative abundance of microbial functional abilities as predicted by PICRUSt. Predator presence was shown to affect the presence of signaling molecules including the cytokine-cytokine receptor interactions, where microbes are suggested to be necessary for the stressor-induced increases in circulating cytokine .
The stress responses in teleost fish are similar to those found in terrestrial vertebrates. The primary response is to release stress hormones into the circulations, for example, corticosteroids that could combine with glucocorticoid to restore the hydromineral homeostasis . The main purpose of our experiment was not aimed at measuring hormones. Nevertheless, a primary hormone response would also start an endocrine response that has been shown to have a close interactions with fish immune systems . Both innate and adaptive immune system were previously suggested to influence gut microbiota composition and diversity [56, 57]. A recent study also gave evidence that guppies living in high-predation and low-predation areas expressed different gut enzymatic profiles . These differences in gut enzymatic profiles can result in changes in abundance and composition of the corresponding microbiota. Sympathetic nervous system (SNS) activities induced by stress could be another significant influence on gut secretion and motility and has been shown to influence the stability of microbiota communities in mice . These multiple molecular studies together with our experimental results thus suggest that the physiological changes from stress can result in an environmental disturbance in the gut environment, which can alter the composition and functioning of the microbial communities .
In an accompanying paper from the experiment of this study, we found that perch had lower body condition in the presence of pike even though they were fed a similar amount of food , suggesting that the microbial change in response to predator presence could have consequences on perch body condition and vice versa. In our study, the perch relative intestine length was significantly shorter with pike presence. Such a decrease in intestine length when facing predation could be the result of reduced food intake . However, during the experiment period, we detected no food residuals after feeding the fish. Thus, other factors than food limitation probably caused the intestine shortening when pike was present. Dealing with stress is an energy demanding process for animals and can have consequences on metabolism. Animals might need to re-distribute metabolic substrates to other tissues so that they can handle the increasing energy needs when facing stress, such as to stimulate oxygen in gills , instead of spending them on a high-energy intestine . Interestingly, we did not find a significant correlation between relative intestine length and the microbial community diversity. This is surprising, as predation-stress decreased microbial diversity as well as intestine length. Our results thus suggest that the host physiological changes from stress and/or the food substrates on microbe use in intestine have bigger effects on gut microbial communities than intestine length per se.
In natural animals, predators strongly influence their prey, including changes in prey’s food intake in terms of both quality and quantity. Food intake from the host is provided as substrates for gut microbes to use for their own growth. Meanwhile, gut microbiota greatly contribute to regulate host energy harvesting , as shown by responses to different diet, and also in response to calorie intake . This is reflected by our experiments with one type of food, where food ration influenced the gut microbiota community in perch. One possible explanation for the effect of food ration on gut microbiota is that a high ration of food favors bacteria that are quick colonizers and fast growers, as food is not limiting. At lower food ratios on the other hand, bacteria that are good competitors will be favored . An extreme situation of lower food ration is starvation, which has been shown to change physiological state in fish to meet energy requirements , and several studies have shown that starvation is a stress factor for fish [20, 64, 65].
When comparing the relative abundance of bacterial phyla among treatments, we found that Fusobacteria increased both at the lowest food ration and at predation stress conditions. This suggests that Fusobacteria, especially Cetobacterium could be used as indicators for fish experiencing stress. Furthermore, Fusobacteria have also been suggested to be associated with many human infections, such as colonic mucosa inflammation [66, 67], inflammatory bowel disease (IBD), where IBD often increase with host psychological stress . In addition, we suggest that changes of other secondary metabolites as an effect of food ration could also be important and possibly contribute to the nutrition absorption, energy obtaining, and weight gaining of the host .
We also found that microbial diversity was affected by the interaction between food ration and sex. For example, gut microbiota in perch male and female responded differently when fed with different amounts of food. The mechanisms behind this pattern are not clear. One explanation can be that males and females differ in intestinal tract physiology. For example in humans, females have a longer transit time in their intestine compared to males , which could give gut microbial community longer time to use the substrate. Another explanation can be related to sex hormones, which have been shown to play important roles in regulating bacterial metabolism and growth . Studies in mice have also shown sex-specific differences in gut microbiota composition at puberty . Similar to our study, Bolnick et al.  found that gut microbiota responded differently to diet in male and female stickleback. The role of sex hormones rests on the assumption that male and female sex hormones affect bacterial growth differently at high and low food rations, which requires further investigations.
In this study, we showed that perch gut microbiota communities react to predation stress and food ration, with parts of the gut community also showing host sex-specific responses. The observed effects of predation and food ration call for an assessment of the role of gut microbiota in food web dynamics and trophic energy transfer.
Cahill MM. Bacterial flora of fishes: a review. Microb. Ecol. 1990;19:21–41. Available from: http://www.springerlink.com/index/10.1007/BF02015051.
Fraune S, Bosch TCG. Why bacteria matter in animal development and evolution. BioEssays. 2010;32:571–80. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20544735.
Dutton RJ, Turnbaugh PJ. Taking a metagenomic view of human nutrition. Curr. Opin. Clin. Nutr. Metab. Care. 2012;15:448–54. Available from: http://www.ncbi.nlm.nih.gov/pubmed/22878238.
Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, et al. Host-gut microbiota metabolic interactions. Science (80-. ). 2012;336:1262–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/22674330.
Spor A, Koren O, Ley R. Unravelling the effects of the environment and host genotype on the gut microbiome. Nat. Rev. Microbiol. 2011;9:279–90. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21407244.
Bolnick DI, Snowberg LK, Caporaso JG, Lauber C, Knight R, Stutz WE. Major histocompatibility complex class IIb polymorphism influences gut microbiota composition and diversity. Mol Ecol. 2014;23:4831–45.
García-Gómez E, González-Pedrajo B, Camacho-Arroyo I. Role of sex steroid hormones in bacterial-host interactions. Biomed Res. Int. 2013 Available from: https://doi.org/10.1155/2013/928290.
Markle JGM, Frank DN, Mortin-toth S, Robertson CE, Feazel LM, Rolle-kampczyk U, et al. Sex differences in the gut microbiome drive hormone-dependent regulation of autoimmunity. Science (80-. ). 2013;339:1084–1088.
Bolnick DI, Snowberg LK, Hirsch PE, Lauber CL, Org E, Parks B, et al. Individual diet has sex-dependent effects on vertebrate gut microbiota. Nat. Commun. 2014;5:4500. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25072318.
Sullam KE, Essinger SD, Lozupone CA, O’Connor MP, Rosen GL, Knight R, et al. Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. Mol. Ecol. 2012;21:3363–78. Available from: http://www.ncbi.nlm.nih.gov/pubmed/22486918.
Gerking SD. Trophic levels and optimal foraging theory. Feed. Ecol. Fish. 1994. p. 3–13. Available from: http://www.sciencedirect.com/science/article/pii/B9780122807800500062
Polis GA, Strong DR. Food web complexity and community dynamics. Am Nat. [University of Chicago Press, American Society of Naturalists]; 1996;147:813–846.
Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbiota. Science (80-. ). 2008;320:1647–1652.
Liu H, Guo X, Gooneratne R, Lai R, Zeng C, Zhan F, et al. The gut microbiome and degradation enzyme activity of wild freshwater fishes influenced by their trophic levels. Sci. Rep. 2016;6:24340. Available from: http://www.nature.com/articles/srep24340
Carpenter SR, Kitchell JF, Hodgson JR. Cascading trophic interactions and lake productivity. Bioscience. 1985;35:634–9.
Tonn WM, Holopainen IJ, Paszkowski CA. Density-dependent effects and the regulation of crucian carp populations in single-species ponds. Ecology. 1994;75:824–34.
Svanbäck R, Persson L. Individual diet specialization, niche width and population dynamics: implications for trophic polymorphisms. J. Anim. Ecol. 2004;73:973–82. Available from: http://www.jstor.org/stable/3505532
Svanbäck R, Eklöv P, Fransson R, Holmgren K. Intraspecific competition drives multiple species resource polymorphism in fish communities. Oikos. 2008;117:114–24.
Crawford PA, Crowley JR, Sambandam N, Muegge BD, Costello EK, Hamady M, et al. Regulation of myocardial ketone body metabolism by the gut microbiota during nutrient deprivation. Proc. Natl. Acad. Sci. U. S. A. 2009;106:11276–81. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19549860.
Baumgarner BL, Bharadwaj AS, Inerowicz D, Goodman AS, Brown PB. Proteomic analysis of rainbow trout ( Oncorhynchus mykiss ) intestinal epithelia: physiological acclimation to short-term starvation. Comp. Biochem. Physiol. - Part D Genomics Proteomics. Elsevier Inc.; 2013;8:58–64. Available from: https://doi.org/10.1016/j.cbd.2012.11.001
Carey HV, Walters WA, Knight R. Seasonal restructuring of the ground squirrel gut microbiota over the annual hibernation cycle. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2013;304:R33–42. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3543654&tool=pmcentrez&rendertype=abstract
Sommer F, Ståhlman M, Ilkayeva O, Arnemo JM, Kindberg J, Josefsson J, et al. The gut microbiota modulates energy metabolism in the hibernating brown bear Ursus arctos. Cell Rep. 2016;1655–61. Available from: http://linkinghub.elsevier.com/retrieve/pii/S2211124716000474
Johansson F, Andersson J. Scared fish get lazy, and lazy fish get fat. J. Anim. Ecol. 2009;78:772–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19302323
Brönmark C, Miner JG. Predator-induced phenotypical change in body morphology in crucian carp. Science (80-. ). 1992;258:1348–50. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17778362
Svanbäck R, Zha Y, Brönmark C, Johansson F. The interaction between predation risk and food ration on behavior and morphology of Eurasian perch. Ecol Evol. 2017;7:8567–77.
Fischer EK, Harris RM, Hofmann HA, Hoke KL. Predator exposure alters stress physiology in guppies across timescales. Horm. Behav. Elsevier Inc.; 2014;65:165–72. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24370688
Collins SM, Bercik P. The relationship between intestinal microbiota and the central nervous system in normal gastrointestinal function and disease. Gastroenterology. AGA Institute American Gastroenterological Association; 2009;136:2003–14.http://www.ncbi.nlm.nih.gov/pubmed/19457424
Candela M, Biagi E, Maccaferri S, Turroni S, Brigidi P. Intestinal microbiota is a plastic factor responding to environmental changes. Trends Microbiol. Elsevier Ltd; 2012;20:385–91. Available from: http://www.ncbi.nlm.nih.gov/pubmed/22672911
Forsythe P, Kunze W, Bienenstock J. On communication between gut microbes and the brain. Curr. Opin. Gastroenterol. 2012;28:557–62. Available from: http://www.ncbi.nlm.nih.gov/pubmed/23010679
O’Mahony SM, Marchesi JR, Scully P, Codling C, Ceolho A-M, Quigley EMM, et al. Early life stress alters behavior, immunity, and microbiota in rats: implications for irritable bowel syndrome and psychiatric illnesses. Biol. Psychiatry. Society of Biological Psychiatry; 2009;65:263–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/18723164
Eklöv P, Jonsson P. Pike predators induce morphological changes in young perch and roach. J Fish Biol. 2007;70:155–64.
Olsson J, Svanback R, Eklov P. Growth rate constrain morphological divergence when driven by competition. Oikos. 2006;115:15–22.
Bolnick DI, Snowberg LK, Hirsch PE, Lauber CL, Knight R, Caporaso JG, et al. Individuals’ diet diversity influences gut microbial diversity in two freshwater fish (threespine stickleback and Eurasian perch). Ecol. Lett. 2014; Available from: http://www.ncbi.nlm.nih.gov/pubmed/24847735
Lessmark O (Lund U (Sweden). LI. Competition between perch (Perca fluviatilis) and roach (Rutilus rutilus) in south Swedish lakes. Limnol. Institutionen, Lunds Univ. 1983.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U. S. A. 2011;108:4516–22. Available from: http://www.pnas.org/content/108/Supplement_1/4516.short
Sinclair L, Osman OA, Bertilsson S, Eiler A. Microbial community composition and diversity via 16S rRNA gene amplicons: evaluating the illumina platform. PLoS One. 2015;10:1–18. Available from: https://doi.org/10.1371/journal.pone.0116955
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods. 2013;10:996–8. Available from: https://doi.org/10.1038/nmeth.2604
Lanzén A, Jørgensen SL, Huson DH, Gorfer M, Grindhaug SH, Jonassen I, et al. CREST––classification resources for environmental sequence tags. PLoS One. 2012;7:e49334.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. Nature Publishing Group; 2013;31:814–21. Available from: http://www.ncbi.nlm.nih.gov/pubmed/23975157
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. The gene ontology consortium. Gene ontology: tool for the unification of biology. Nat. Genet. 2011;25:25–9.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26:139–40.
Skrodenyte-arbaciauskiene V. Enzymatic activity of intestinal bacteria in roachRutilus rutilusL. Fish. Sci. 2007;73:964–6. Available from: http://link.springer.com/10.1111/j.1444-2906.2007.01421.x
Rekecki A, Dierckens K, Laureau S, Boon N, Bossier P, Van den Broeck W. Effect of germ-free rearing environment on gut development of larval sea bass (Dicentrarchus labrax L.). Aquaculture. Elsevier B.V.; 2009;293:8–15. Available from:Dicentrarchus labrax). Aquaculture. Elsevier B.V.; 2009;293:8–15. Available from:https://doi.org/10.1016/j.aquaculture.2009.04.001
Bairagi A, Ghosh KS, Kumar S, Ray AK. Enzyme producing bacterial flora isolated from fish. Aquac Int. 2002;10:109–21.
Knowles SR, Nelson EA, Palombo EA. Investigating the role of perceived stress on bacterial flora activity and salivary cortisol secretion: a possible mechanism underlying susceptibility to illness. Biol. Psychol. 2008;77:132–7. Available from:http://www.ncbi.nlm.nih.gov/pubmed/18023961
Bailey MT, Dowd SE, Galley JD, Hufnagle AR, Allen RG, Lyte M. Exposure to a social stressor alters the structure of the intestinal microbiota: implications for stressor-induced immunomodulation. Brain. Behav. Immun. 2011;25:397–407. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3039072&tool=pmcentrez&rendertype=abstract
Abrahams MV, Dill LM. A determination of the energetic equivalence of the risk of predation. Ecology. 1989;70:999–1007.
Heithaus MR, Dill LM. Food availability and tiger shark predation risk influence bottlenose dolphin habitat use. Ecology. 2002;83:480–91.
Eklöv P, Svanbäck R. Predation risk influences adaptive morphological variation in fish populations. Am. Nat. 2006;167:440–52. Available from:http://www.ncbi.nlm.nih.gov/pubmed/16673351.
de Vries EJ, Jacobs G, Sabelis MW, Menken SBJ, Breeuwer JAJ. Diet-dependent effects of gut bacteria on their insect host: the symbiosis of Erwinia sp. and western flower thrips. Proc. Biol. Sci. 2004;271:2171–8. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1691834&tool=pmcentrez&rendertype=abstract
Ley RE, Bäckhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. Proc. Natl. Acad. Sci. U. S. A. 2005;102:11070–5. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1176910&tool=pmcentrez&rendertype=abstract
Nielsen S, Nielsen DS, Lauritzen L, Jakobsen M, Michaelsen KF. Impact of diet on the intestinal microbiota in 10-month-old infants. J. Pediatr. Gastroenterol. Nutr. 2007;44:613–8. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17460496
De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc. Natl. Acad. Sci. U. S. A. 2010;107:14691–6. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2930426&tool=pmcentrez&rendertype=abstract
Desai AR, Links MG, Collins SA, Mansfield GS, Drew MD, Van Kessel AG, et al. Effects of plant-based diets on the distal gut microbiome of rainbow trout (Oncorhynchus mykiss). Aquaculture. 2012;350–353:134–42. Available from: http://linkinghub.elsevier.com/retrieve/pii/S0044848612002098
Wendelaar ES. The stress response in fish. Physilogical Rev. 1997;77:592–616.
Vijay-kumar M, Aitken JD, Carvalho FA, Cullender TC, Mwangi S, Srinivasan S, et al. Metabolic syndrome and altered gut microbiota in mice lacking toll-like receptor 5. Science (80-. ). 2010;328:228–231.
Zhang H, Sparks JB, Karyala S V, Settlage R, Luo XM. Host adaptive immunity alters gut microbiota. ISME J. Nature Publishing Group; 2015;9:770–81. Available from: http://www.nature.com/doifinder/10.1038/ismej.2014.165
Sullam KE, Dalton CM, Russell JA, Kilham SS, El-Sabaawi R, German DP, et al. Changes in digestive traits and body nutritional composition accommodate a trophic niche shift in Trinidadian guppies. Oecologia. 2014;177:245–57. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25430044
Bailey MT, Dowd SE, Parry NMA, Galley JD, Schauer DB, Lyte M. Stressor exposure disrupts commensal microbial populations in the intestines and leads to increased colonization by Citrobacter rodentium. Infect. Immun. 2010;78:1509–19. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2849416&tool=pmcentrez&rendertype=abstract
Werner EE, Gilliam JF, Hall DJ, Mittelbach GG, Gilliam2 JF, Mittelbach3 GG. An experimental test of the effects of predation risk on habitat use in fish. Ecology. 1983;64:1540–8. Available from: http://www.jstor.org/stable/1937508?seq=1#page_scan_tab_contents
Cant JP, McBride BW, Croom WJ. The regulation of intestinal metabolism and its impact on whole animal energetics. J Anim Sci. 1996;74:2541–53.
Bäckhed F, Ding H, Wang T, Hooper L V, Koh GY, Nagy A, et al. The gut microbiota as an environmental factor that regulates fat storage. Proc. Natl. Acad. Sci. U. S. A. 2004;101:15718–23. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=524219&tool=pmcentrez&rendertype=abstract
Johnston IA. Quantitative studies of muscle breakdown during starvation in the marine flat fish Pleuronectes platessa. Cell Tissue Res. 1981;214:369–86.
Montgomery WL, Pollak PE. Gut anatomy and pH in a Red Sea surgeonfish, Acanthurus nigrofuscus. Mar Ecol - Prog Ser. 1988;44:7–13.
Hall K, Bellwood D. Histological effects of cyanide, stress and starvation on the intestinal mucosa of Pomacentrus coelestis, a marine aquarium fish species. J Fish Biol. 1995;47:438–54.
Bennett KW, Eley A. Fusobacteria: new taxonomy and related diseases. J Med Microbiol. 1993;39:246–54.
McCoy AN, Araújo-Pérez F, Azcárate-Peril A, Yeh JJ, Sandler RS, Keku TO. Fusobacterium is associated with colorectal adenomas, e53653. PLoS One. 2013;8(1)
Mawdsley JE, Rampton DS. Psychological stress in IBD: new insights into pathogenic and therapeutic implications. Gut. 2005;54:1481–91.
Ridaura VK, Faith JJ, Rey FE, Cheng J, Duncan AE, Kau AL, et al. Cultured gut microbiota from twins discordant for obesity modulate adiposity and metabolic phenotypes in mice. Science (80-. ). 2013;341:1241214–1–1241214–10. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3829625&tool=pmcentrez&rendertype=abstract
Freire AC, Basit AW, Choudhary R, Piong CW, Merchant HA. Does sex matter? The influence of gender on gastrointestinal physiology and drug delivery. Int. J. Pharm. 2011;415:15–28. Available from: https://doi.org/10.1016/j.ijpharm.2011.04.069
We thank Konrad Karlsson and Dandan Shen for their help in the field and in the laboratory. We thank Peter Eklöv, Eva Lindström, and Raven Bier for their constructive comments on earlier versions of this manuscript and Brendan Bohannan for his nice discussions and suggestions on this study.
This work was funded by the Swedish Research Council (VR – 2011-05646) to RS, the Swedish Research Council (VR - 2012-4592), and the Swedish Foundation for strategic research (grant ICA10-0015) to AE.
Availability of data and materials
The data of raw sequences for microbial sequence data were deposited in the GenBank sequence read archive (SRA) under accession number SRP090413.
The experiments were carried out in accordance with the recommendations of the ethical committee of Uppsala Djurförsöketiska Nämnd (permit number C80/13).
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Measurements of weight, length, intestine length, and sex of all perch used. (XLSX 45 kb)
Quality control results for PICRUSt showing the percentage of successful reads that were mapped to Greengenes when using the closed reference OTU picking. (XLSX 46 kb)
Quality control results for PICRUSt showing NSTI scores of the reference genome coverage for each fish sample. (XLSX 48 kb)
ANOVA test of effects of food ration, predation stress, sex, and their two-way interactions on the ten most abundant phyla in intestinal microbiota community. Significant treatment effects are highlighted in bold text. (DOCX 69 kb)
Relative abundance changes of the top genus from phyla Bacteroidetes, Cyanobacteria, Fusobacteria, and Proteobacteria across different factors. (a) Changes of the relative abundance of the top genus in Bacteroidetes affected by the interaction of food ration (5, 10, and 15%) and host sex. (b) Changes of the relative abundance of the top genus in Cyanobacteria affected by pike presence (yes pike) and pike absence (no pike). (c) Changes of the relative abundance of the top genus in Fusobacteria affected by pike presence (yes pike) and pike absence (no pike). (d) Changes of the relative abundance of the top genus in Proteobacteria affected by pike presence (yes pike) and pike absence (no pike). (PDF 615 kb)
Alpha diversity changes across the factors. (a) Phylogenetic diversity (PD) affected by pike absence (N) and pike presence (Y). (b) Observed species richness affected by predation absence (N) and pike presence (Y). p values were obtained from TukeyHSD test with ANOVA model. (PDF 119 kb)
Results of ANOVA testing effects of food ration, predation stress, sex, and their two-way interactions on intestinal microbiota diversity indices Chao1, phylogenetic diversity (PD), and observed species richness (S. Obs). Significant treatment effects are highlighted in bold text. (DOCX 54 kb)
Two-dimensional non-metric multidimensional scaling (NMDS) plot of bacterial communities. Point patterns denote bacterial communities from food (round), perch intestine (triangle), and water (square). (a) NMDS plot generated by using unweighted UniFrac distance matrix. (b) NMDS plot generated by using weighted UniFrac distance matrix. (PDF 184 kb)
Results of EdgeR showing the representative OTUs in treatments. (XLSX 39 kb)
Sum of relative abundance of the representative OTUs clustered into phyla level. (a) Relative abundance of representative phyla across the interaction of 10% food ration and pike predation. (b) Relative abundance of representative phyla across the pike predation treatment. (PDF 232 kb)
ANOVA test of food ration and predation stress effect on functional categories. Significant treatment effects are highlighted in bold text. (DOCX 119 kb)