- Research article
- Open Access
Modulation of posterior intestinal mucosal proteome in rainbow trout (Oncorhynchus mykiss) after Yersinia ruckeri infection
Veterinary Researchvolume 50, Article number: 54 (2019)
Yersinia ruckeri is the causative agent of enteric redmouth disease in salmonids. In fish, the intestine represents an important site of nutrient uptake, host–pathogen interactions, and defense. The posterior intestine can be inflamed, reddened, and filled with an opaque, yellowish fluid during Y. ruckeri infection. Herein, we report an investigation on the proteome alteration in the posterior intestinal mucosa of rainbow trout (Oncorhynchus mykiss) after exposure to Y. ruckeri. The intestinal mucosal proteins were identified and quantified by a shotgun proteomic approach by applying data-independent quantification with sequential windowed acquisition of all theoretical mass spectra (SWATH). A total of 437 proteins were found to be differentially up- or downregulated in the posterior intestine. Gene ontology of upregulated proteins pointed to their involvement into exopeptidase, endopeptidase, and hydrolase activities, while the downregulated proteins were involved in lipid metabolism, actin binding, and translation processes. Additionally, upregulated proteins were predicted to be involved in lysosome, oxidative phosphorylation, and metabolic pathways, while downregulated proteins were implicated in focal adhesion, regulation of actin cytoskeleton, protein digestion and absorption pathways. This study showed that Y. ruckeri infection can alter protein abundance involved in serine-type carboxypeptidase, cysteine and aspartic-type endopeptidases, metallopeptidases, antioxidant defense, calcium ion binding, glycolytic and carbohydrate metabolic processes in the proteome of the intestinal mucosa of rainbow trout.
Enteric redmouth disease (ERM) is one of the major diseases of mainly salmonid fish. It was first described in the 1950s in rainbow trout (Oncorhynchus mykiss) in the Hagerman Valley of Idaho, USA. The disease is caused by Yersinia ruckeri, a Gram-negative rod shaped enterobacterium, and prevalent in Europe, North and South America, Australia, South Africa, the Middle East, and China [1, 2]. The signs of the disease include exophthalmia, subcutaneous hemorrhages, and splenomegaly . In most cases, the intestine is inflamed, reddened, and filled with an opaque, yellowish fluid [1, 3]. In rainbow trout, up to 25% of infected fish could carry Y. ruckeri in their posterior intestines . Y. ruckeri appeared to adhere or invade gills and gut tissue of rainbow trout [3, 5]. Y. ruckeri can also be seen in the lumen of the intestine, adhering to mucus and invading villi within 30 min after exposure and can be detected 7 days post-exposure (dpe) in the liver, spleen, brain, and heart .
The strains of Y. ruckeri have been classified on the basis of biotypes, serotypes (O1, O2, O3, and O4), and outer-membrane protein types . The majority of epizootics in salmonids are caused by the motile serotype O1a . There are two biotypes of Y. ruckeri; biotype 1 is positive for motility and lipase secretion, whereas biotype 2 is negative for both tests . The outbreaks of ERM were recorded in biotype 1-vaccinated salmonids and the isolates associated with these outbreaks were identified as Y. ruckeri biotype 2 [9,10,11,12].
In fish, the gastrointestinal (GI) tract is a multifunctional organ, which serves a diverse range of functions from nutrient absorption to ionic and osmotic regulation and even air breathing [13,14,15,16,17,18,19]. Gut epithelial cells are protected by a mucus layer, which creates a physical and chemical barrier against an intruder and acts as an important mechanism of innate defense that maintains tissue homeostasis . The GI tract is continuously challenged with food antigens as well as pathogens entering the body via feed and water intake, and acts as the first line of defense against pathogen attachment and invasion [20,21,22]. The posterior portion of the GI tract is a site where uptake of macromolecules and foreign antigens occurs, which leads to antigen uptake and processing as well as the initiation of a systemic immune response [20, 23]. The disease has been of a major concern in salmonids, which influences the physiological response of fish. Rainbow trout is highly susceptible to ERM . The effect of natural herbal products on growth performance, hematological and biochemical values and resistance to Y. ruckeri infection was studied as a feed additives for rainbow trout [24, 25].
In ERM, hyperemia and hemorrhages throughout the intestinal mucosa in rainbow trout have been observed . The intestinal response of rainbow trout has been mainly studied at the transcription level, especially regarding immune gene expression [26,27,28,29]. It is still largely unknown how the gut mucosal proteome of the rainbow trout responds to alternations in the luminal environment triggered by Y. ruckeri. Therefore, in the current study, we aimed to evaluate the possible effects of Y. ruckeri exposure on the modulation of the proteome profile in the posterior intestine of rainbow trout. In order to meet this aim, we applied a label-free shotgun proteomic approach.
Materials and methods
Collecting posterior intestine samples
The details of the experimental setup have been described in our previous study . Briefly, specific pathogen free rainbow trout (15 ± 1 cm) were allocated to 9 aquaria, 18 fish per aquarium. In order to monitor the effect of Y. ruckeri exposure, the following three treatments were given in parallel—Treatment 1: an exposure with biotype 1 (CSF007–82), Treatment 2: an exposure with biotype 2 (A7959–11), and Treatment 3: a mock exposure to sterile tryptic soya broth. In both the test groups, the fish were challenged by bath exposing to 2 × 106 colony forming units (CFU) of Y. ruckeri biotype 1 (CSF007-82) and Y. ruckeri biotype 2 (A7959-11) strains for 2 h. The fish were maintained in a flow-through system supplied with UV-treated ground water at 19 ± 1 °C, monitored daily, and morbid, moribund, and dead fish were taken out immediately from the aquaria. The mortality of fish was only considered to be caused by Y. ruckeri if the bacteria were recovered from the head kidney and confirmed using the MONO-Yr kit (Bionor, Skien, Norway) or PCR . No morbid and moribund fish were sampled for the study. The fish were maintained in starvation 48 h before the sampling. Nine fish from each group were sampled at 3, 9, and 28 dpe. The posterior intestine of each sample was dissected and placed in cold sterile PBS (Sigma-Aldrich, Neustadt, Germany). The intestine was opened along the mesenteric border and washed three times with cold sterile PBS containing a cocktail of mammalian protease inhibitors (Sigma-Aldrich) to remove digesta and fecal matter. Intestinal mucosa was scraped with a sterile large scalpel blade. Each sample was divided into two parts, one immediately snap-frozen in liquid nitrogen for proteomic analysis, and one fixed in RNALater (Sigma-Aldrich) for molecular analysis and stored at −80 °C.
Equal amounts (30 mg) of nine individual intestinal mucosal samples from each group were pooled randomly as three pools of three samples each in order to minimize the effects of individual variation. It refers to N = 3 for the biological replicates per time point (3, 9 and 28 dpe) and per exposed group (control, Y. ruckeri biotype 1 and biotype 2). Each pool was ground using a sterilized mortar and pestle in the presence of liquid nitrogen to a fine powder, which was then mixed with 800 μL of precooled denaturing lysis buffer (7 M urea, 2 M thiourea, 4% CHAPS, and 1% DTT) containing mammalian protease inhibitor cocktail (Sigma-Aldrich). The samples were subjected to ultrasonic disruption for 10 cycles of 10 s pulse-on and 30 s pulse-off. The samples were then centrifuged at 12 000 rpm for 20 min at 4 °C to remove any cellular debris. The supernatant was collected. The protein content was measured using the Pierce 660 nm Protein Assay (Thermo Scientific, Vienna, Austria).
The protein digestion was performed following a standard enhanced, filter-aided sample preparation protocol (FASP) . Briefly, after washing, proteins were reduced and alkylated. On-filter digestion was performed with 1.2 µg Trypsin/Lys-C mix (Promega, Madison, USA) for 14 h at 37 °C. Digested peptides were recovered, dried, and redissolved in 0.1% aqueous trifluoroacetic acid. After desalting and sample cleanup by C18 ZipTips (Sigma-Aldrich), all samples were spiked with standardized indexed retention time reference peptides for facilitation of retention time alignment (iRT kit, Biognosys, Switzerland). For each liquid chromatography (LC)–mass spectrometer (MS) analysis, 3 µg of digested protein was loaded onto the LC column.
Microflow high-performance liquid chromatography electrospray-ionisation quadrupole time-of-flight tandem mass spectrometry (micro LC ESI QTOF MS/MS)
Micro LC ESI QTOF MS/MS was carried out as described earlier . Briefly, peptides were separated on an Eksigent NanoLC 425 system using a microflow pump module (Sciex, Concord, Canada). The samples were pre-concentrated and desalted on a 5 mm YMC-Triart C18 precolumn using ultra-pure LC–MS grade water with 0.1% formic acid as a mobile phase and a flow rate of 10 μL/min. Desalted peptides were separated on a 15 cm YMC-Triart C18 column at a flow rate of 5 µL/min. Mobile phase A was ultra-pure water with 0.1% formic acid, whereas mobile phase B consisted of acetonitrile with 0.1% formic acid. The gradient started with 3% B and increased in two steps to 25% B (68 min) and 35% (73 min) followed by a washing step with 80% B. Total acquisition time was 87 min. For mass spectrometric analysis, LC was directly coupled via a DuoSpray ion source in electrospray mode to a high resolution quadrupole time-of-flight mass spectrometer (Triple TOF 5600+, Sciex). For information dependent data acquisition (IDA runs), the MS1 spectra were collected in a range of 400–1250 m/z for 250 ms. The 40 most intense precursors with charge state 2–4, which exceeded 150 counts per second, were selected for fragmentation. MS2 spectra were collected in the range of 200–1500 m/z for 50 ms. The precursor ions were dynamically excluded from reselection for 13 s. Based on the data-dependant acquisition spectra, an ion library was established for the next step of data-independent acquisition of SWATH spectra (Sequential Windowed Acquisition of All Theoretical Mass Spectra). Seventy-three variable windows were created in a mass range of 400–1250 Da depending on the precursor ion density. All precursors of each window were fragmented. MS2 spectra were acquired for 50 ms (SWATH runs).
The acquired raw data were processed with ProteinPilot Software version 5.0 (Sciex) for re-calibration and database searches as described by Kumar et al. . The database contained entries of following taxonomies: Oncorhynchus mykiss (NCBI, Refseq: 71 285 entries), and Y. ruckeri (UniProt, taxonomy id 29486: 4493 entries) as well as cRAP (common Repository of Adventitious Proteins). The database search parameters applied were trypsin digestion, cysteine alkylation set to iodoacetamide, search effort set to rapid ID. The mass tolerance in MS mode was set by program default with 0.05 Da in MS and 0.1 Da in MSMS mode for the rapid recalibration search, and 0.0011 Da in MS and 0.01 Da in MSMS mode for the final search. False discovery rate (FDR) analysis was performed using the integrated tools in ProteinPilot with < 1% on peptide as well as on protein level.
For the quantification of proteins by SWATH, the identified proteins of a combined search of all 27 IDA runs identifying 3425 proteins served as a basis for the creation of the SWATH ion library with the MS/MS (ALL) with SWATH Acquisition MicroApp 2.0 in PeakView 2.2 (both Sciex) (dataset PRIDE PXD011087, IDA FASP intestine Control/CSF/7959). Shared peptides were excluded. So the resulting ion library (Additional file 1) contained a total of 3372 proteins, 25 469 peptides and 295 168 transitions to be quantified potentially. The key criteria for processing of the SWATH samples were to use only proteins with an FDR rate below 1%, furthermore up to 6 peptides per protein and up to 6 transitions per peptide were chosen, false discovery rate threshold was 5%, peptide confidence threshold 98%, XIC extraction window 5 min, XIC width 100 ppm and modified peptides were excluded. Retention time alignment based on iRT peptides (iRT kit, Biognosys) and processing of SWATH raw samples for calculation of raw peak areas were performed in PeakView 2.2. These were normalized in the software MarkerView 1.2.1 (Sciex) with the integrated tool based on the total area sums, which assumes that the abundance of most of the proteins within the label free quantification approach remains unchanged and only a few of them are differentially regulated.
All statistical analyses were performed in R programming language . The protein abundance after retention time alignment and normalization on total area sums derived from MarkerView were first transformed into a logarithmic scale. Technical replicates compensating for instabilities during LC–MS acquisition (N = 2) were averaged for the mean values. Afterward, the differences in abundance of posterior intestinal proteins were assessed using one-way ANOVA for each protein comparing groups [(control vs. biotype 1 (CSF007-82) exposed vs. biotype 2 (A7959-11) exposed posterior intestine samples] for each time point (day 3, 9 and 28). The method of Benjamini and Hochberg  was applied to control the FDR. The differences between the biotypes-exposed and control posterior intestine samples were considered significant if FDR-adjusted p-values were smaller than the significance level of α = 0.05. In order to assess the significance of the pairwise comparisons, Tukey’s honest significant difference method was applied as a post hoc test. Protein level changes were considered differential if the adjusted p value was below α and the absolute fold change was at least two (fold change < −2 or > +2).
PANTHER classification system was used to categorize the Gene Ontology of all the differentially up- or downregulated proteins in biological process, cellular component, and protein class . Molecular function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway of differentially up- or downregulated proteins were further determined on the Cytoscape software (version 3.7.1)  plugin ClueGO (version 2.5.4) . We tested the significance of GO term and KEGG pathway using Fisher’s exact test and FDR-correction at< 0.05 significant level.
Quantitative real time PCR
In order to validate the abundance of identified proteins, eight kinds of proteins (probable serine carboxypeptidase, cathepsin D, caspase 6, lysozyme C II, precerebellin, protein S100, and tubulin alpha) were chosen for qPCR analysis carried out on a CFX96 Touch Real-Time PCR detection system (Bio-Rad, Hercules, USA). The sequences of the primers used in this study are listed in Additional file 2. Total RNA was extracted from posterior intestinal mucosa samples using an RNeasy Mini Kit (Qiagen, Hilden, Germany) and included an on-column DNase digestion step. cDNA was synthesized using an iScript cDNA Synthesis Kit (Bio-Rad) with 500 ng total RNA according to the user’s manual.
All reactions were done in a final volume of 20 μL, which contained 4 μL of 1:10-fold diluted cDNA, 0.5 μL of each primer (10 pmol/μL), 10 μL of 2× SsoAdvanced™ Universal SYBR Green Supermix (Bio-Rad), and 5 μL of RNase-free water. After 5 min of cDNA denaturation at 95 °C, 38 cycles were performed at 95 °C for 30 s, 57 °C for 30 s and 72 °C for 30 s. Each qPCR was performed for all biological replicates. Relative gene expression was assessed by the CFX Manager Maestro Software in normalized expression mode (∆∆Cq), using elongation factor alpha 1 as a Ref. , and was compared to the control sample. Additionally, standard curve was generated for qPCR using plasmid containing Y. ruckeri 16S ribosomal RNA fragment . Y. ruckeri load was measured in the exposed and unexposed control posterior intestine samples with a Y. ruckeri 16S rRNA specific primers . The copy number of Y. ruckeri 16S rRNA was calculated for each sample. The differences between unexposed control and exposed group and between exposed groups were tested for significant differences using a t-test. The sequential differences within exposed groups were tested for using a one-way ANOVA and significant differences revealed with the Tukey’s post hoc test. The Pearson correlation coefficient (r) was calculated to estimate the relationship between qPCR and proteomic quantifications. For all statistical tests, a p-value of < 0.05 was regarded as significant and all the data were analyzed in IBM SPSS software version 24.
Clinical signs and Y. ruckeri load
The Y. ruckeri-exposed fish were prepared as described in our previous study , where 30% fish showed severe clinical signs of the disease and the maximum morbidity and mortality were observed between 8 and 10 dpe. The fish exhibited external hemorrhages in the caudal and anal fins (Figure 1A) and internal signs, such as enlarged spleen and reddened intestine (Figure 1B), at 9 dpe. However, at 28 dpe, these organs were in normal form and did not show any clinical signs of the disease (Figure 1C). The Y. ruckeri 16S transcript increased and peaked in the fish exposed to both biotypes at 3 dpe and then decreased at 9 dpe and was almost undetected at 28 dpe (Figure 2). No Y. ruckeri was detected in unexposed control samples.
A total of 3006 proteins in the posterior intestine of rainbow trout were identified at 1% FDR and minimum of two matching peptides (Additional file 3) over all 27 samples including controls as well as the exposed groups. In total, 1880 proteins with at least two unique peptides could be quantified. Sixty-two percent of the total 3006 proteins identified amongst all samples could be quantified, after the exclusion of shared and modified peptides and the exclusion of so-called one-peptide wonders. The statistical analysis (exposed versus control) revealed a total of 437 differentially abundant proteins in the posterior intestine at 3 dpe (Additional file 4). Of these, 205 proteins were upregulated and 232 proteins were downregulated in the posterior intestine at 3 dpe. We did not find statistically different changes in the protein abundance at 9 and 28 dpe. As shown in Table 1, upregulated proteins were related to serine -type carboxypeptidase, cysteine-type endopeptidase, aspartic-type endopeptidase, metallodipeptidase, metallopeptidase, hydrolase, peroxidase activity, signal transduction, and carbohydrate metabolic activities. On the other hand, the downregulated proteins were related to lipid metabolism, actin filament binding, stress response, and translation (Table 2).
Within the classification of biological process, the differentially up- or downregulated posterior intestine proteins were mainly associated with metabolic process, biological regulation, cellular processes, and cellular component organization (Figure 3A). In terms of molecular function, the majorities of up- or downregulated posterior intestine proteins were involved in cation-transporting ATPase, exopeptidase, lipase inhibitor, cysteine-type endopeptidase, oxidoreductase, and beta-N-acetylhexosaminidase activities (Figure 3B). Most of the proteins were localized in the cell, organelle, membrane, and extracellular region (Figure 3C). Classification based on protein class resulted in the identification of ten major categories: hydrolase, cytoskeletal protein, oxidoreductase, enzyme modulator, transferase, transporter, calcium binding protein, signaling molecule, nucleic acid binding, and cell adhesion molecule (Figure 3D). Additionally, upregulated proteins were predicted to be involved in lysosome, oxidative phosphorylation, and metabolic pathways, while downregulated proteins were involved in focal adhesion, regulation of actin cytoskeleton, protein digestion and absorption pathways.
Validation of differentially abundant proteins
The upregulation of six proteins (probable serine carboxypeptidase, cathepsin D, caspase 6, lysozyme C II, and precerebellin-like protein) and downregulation of two proteins (protein S100 and tubulin alpha) in the posterior intestine of the fish exposed to both biotypes were confirmed at the mRNA level by qPCR analysis. The transcript levels of these candidate genes were significantly (p < 0.05) either increased or decreased compared to the control posterior intestine samples (Figure 4). The results of qPCR were consistent with those obtained from the proteomic results mainly at 3 dpe (Figure 5), which confirmed reliability of our proteomic data (r = 0.406 or 0.877). Nevertheless, the transcript expression slightly varied in parallel to the corresponding protein abundance perhaps suggesting additional post-transcriptional regulation.
The gut immune system of teleost fish deserves special attention in the biology of the fish intestine. This provides important insights into the mechanisms induced by the progression of diseases . In the case of common carp (Cyprinus carpio), the response of intestinal mucosa against Aeromonas hydrophila was investigated at the proteome level , which improved the understanding of defense mechanisms of carp intestinal mucosa and associated molecular mechanisms. In rainbow trout, the intestinal response was monitored at transcription levels against pathogens like Y. ruckeri, Enteromyxum leei, Edwardsiella ictaluri, and Vibrio harveyi [26,27,28,29]. In the current study, the elements of the gut mucosal responses were examined by proteome profiling of the posterior intestine of rainbow trout following experimental exposure to Y. ruckeri. Additionally, proteins were significantly up- or downregulated in the posterior intestine mainly at 3 dpe, which correlated with increasing bacterial transcripts found in the posterior intestine at 3 dpe (Figure 2). This suggests that Y. ruckeri might affect the gut and modulate the intestinal proteome of fish whenever a high load of bacteria is present in the fish intestine.
In total, 30% mortality was found in rainbow trout following the exposure to Y. ruckeri (2 × 106 CFU) . This mortality rate is in contrast with the previous studies [41, 42], where mortality in rainbow trout was observed to be 20% and 66%, respectively, at lower (5 × 105 CFU) and higher (1.8 × 109 CFU) levels of the challenge dose of Y. ruckeri. These differences in the mortality data may be due to variations between the bacterial strains, challenge doses, fish size, and laboratory conditions. At 28 dpe, the intestine of rainbow trout exposed to Y. ruckeri strains was in normal form and also transcripts of Y. ruckeri 16S were almost undetectable in the posterior intestine. It seems that rainbow trout slowly cleared the Y. ruckeri during the analyzed period.
We found upregulation of exopeptidase and endopeptidase [probable serine carboxypeptidase (3.4- to 4.7-fold), dipeptidase I (2.1- to 2.7-fold), Xaa-Pro dipeptidase (2.4- to 3.1-fold), dipeptidyl peptidase 1 (2.3-fold), cathepsin D (2.4- to 3.2-fold), pro-cathepsin H (4.3- to 6.3-fold) and cathepsin L1 (4.3- to 5.9-fold)] in the posterior intestine of rainbow trout exposed to Y. ruckeri. These exo- and endopeptidases participate in several physiological and cellular processes in the guts of animals . The expression of these peptidase (dipeptidase I and cathepsin D) genes was observed in the gut of channel catfish (Ictalurus punctatus) and grass carp (Ctenopharyngodon idella) following Edwardsiella ictaluri and Aeromonas hydrophila challenges [44, 45], thereby suggesting that they participate in various degradation functions and digestive processes and are closely involved in the immune response of rainbow trout. Overall, our study suggests that proteolysis and peptide hydrolysis activities are upregulated in the GI tract during the bacterial infection, indicating that these identified peptidases might play some role in host defense.
Furthermore, we found upregulation of phagocytosis proteins such as lysozyme C II, macrophage mannose receptor 1, eosinophil peroxidase, and thioredoxin reductase 1 in the posterior intestine of rainbow trout. These phagocytosis proteins were also upregulated in head kidney and spleen of rainbow trout exposed to Y. ruckeri strains . Lysozymes are important defense proteins of the innate immune system of fish against bacterial pathogens [24, 46]. The upregulation of lysozyme C II (3.3-fold) in the posterior intestine, suggests that, besides its role in defense against bacterial pathogens, this protein may also play a significant role in the digestion during the infection in fish. In addition, the upregulation of macrophage mannose receptor 1 (3.2- to 4.7-fold) and thioredoxin reductase 1 (2.2-fold) in the posterior intestine of rainbow trout exposed to Y. ruckeri showed that the abundance of phagocytosis in the gut of rainbow trout results in enhanced intestinal defense against Y. ruckeri. Finally, these results suggest that Y. ruckeri triggers intestinal inflammation which enhances growth in the intestinal lumen and supports the invasion of the intestinal epithelium and mucosal macrophages. This macrophagic process promotes phagocytic activity in the gut to eliminate the bacteria and plays an important role in maintaining gut homeostasis and immune response .
In addition, caspase 6 and precerebellin-like protein were upregulated in the posterior intestine. Caspases are a family of proteases engaged in various important biological processes and play a significant role in the execution phase of the apoptotic death cascade . Caspase-6 was found to be involved in apoptosis and immune response in puffer fish (Takifugu obscurus) against Aeromonas hydrophila . The upregulation of caspase-6 (2.2-fold) in the gut mucosa of rainbow trout suggests that caspase-6 plays a role in the immune response and cell apoptosis against Y. ruckeri infection. However, another apoptosis related protein, programmed cell death protein 5 (PDCD5), was downregulated (−2.7- to 4.4-fold) in the posterior intestine as well. This indicates that PDCD5 may have an important role in the pathogenesis and development of the disease. However, the mechanisms underlying its apoptotic function are largely unknown in the fish. Precerebellin-like protein is a part of the acute phase response and has been characterized in rainbow trout . The expression of precerebellin-like protein gene was upregulated (sevenfold) in the liver of rainbow trout in response to Y. ruckeri at 3 dpe . This suggests that acute phase proteins act as a defensive agent against Y. ruckeri on the gut mucosal surface of rainbow trout. In support of this, the expression of some acute phase proteins such as serum amyloid protein and hepcidin was observed to be significantly increased in rainbow trout after Y. ruckeri challenge .
Evidently, calcium is engaged in signal transduction by acting as a second messenger . In our study, we found two calcium binding proteins (calretinin and regucalcin) and seven signal transduction proteins (Ras-related protein Rab-1B, Rab-5C, Rab7-like, Rab-32-like, Rab-25-like, Rab7-like and ORAB-1) being upregulated. Calretinin is a calcium-binding protein found broadly distributed in the central nervous system and regucalcin is crucial for the regulation of Ca2+ ion homeostasis [52, 53]. The distribution of calretinin immunoreactivity in the developing olfactory system of the rainbow trout was observed by using an immunocytochemistry technique , and the expression of regucalcin gene was upregulated in rainbow trout at different temperatures and Aeromonas salmonicida challenges . It could be speculated that calretinin and regucalcin execute their functions indirectly by controlling Ca2+ homeostasis in infected rainbow trout. Additionally, four calcium binding proteins [annexin A-13 (−4.1-fold), protein S100 (−6.6-fold), M-calpain (−4.0-fold), and calpain 2 catalytic subunit (−2.6-fold)] were downregulated. Annexins are a member of a multigene family of Ca2+ and phospholipid binding proteins . S100 proteins take part in many cellular processes such as modulation of protein kinases and signal transduction pathways, maintenance of cell shape, regulation of calcium homoeostasis  and interaction with the parasite, Tetracapsuloides bryosalmonae . However, calpains are calcium regulated proteases involved in many cellular functions and have been characterized in rainbow trout . In conclusion, the presented results establish the significance of calcium binding proteins and their downregulation in the gut of rainbow trout in response to Y. ruckeri.
Yersinia ruckeri causes profound alterations in the nutritional and metabolic status of the fish [18, 60]. The interaction of a pathogen with intestinal mucosa leads to a variety of physiological responses aimed at adjusting to the new condition and triggers different processes in the gut epithelial cells . Thus, we found upregulation of anabolic and catabolic proteins [fructose-1,6-bisphosphatase 1 (3.8-fold), beta-hexosaminidase (2.6-fold), alpha-N-acetylgalactosaminidase (4.1-fold), and long-chain specific acyl-CoA dehydrogenase (2.6-fold)], transport proteins [epididymal secretory protein E1 (3.9-fold), v-type proton ATPase subunit a (2.2-fold), and AP-1 complex subunit mu-2 (3.1-fold)]. The up-regulation of fructose-1,6-bisphosphatase and epididymal secretory protein E1 was observed in the spleen of rainbow trout in response to Y. ruckeri . This suggests that these proteins are involved in the gut’s immune response and might result in physiological adaptations that contribute to extended longevity of fish during bacterial infection. Additionally, we found some downregulated proteins in the posterior intestine that are involved in gluconeogenesis [l-lactate dehydrogenase (−7.4-fold)], lipid transport [apolipoprotein A-I-1 (−2.2-fold), apolipoprotein A-II (−4.0-fold), and apolipoprotein C-I (−4.6-fold)], and translation [eukaryotic translation initiation factor 5A (−2.4-fold), 40S ribosomal protein S29 (−3.4-fold), and 40S ribosomal protein S12 (−2.2-fold)]. This suggests that fish might adjust their metabolism to channel the energy to the process of defense by decreasing glycogenesis, lipid metabolism, and translation, and increasing carbohydrate and protein catabolism during yersiniosis.
In conclusion, this study provides the first evidence illustrating the proteomic alteration of the intestinal mucosa of rainbow trout in response to Y. ruckeri, suggesting Y. ruckeri exerts a profound impact on posterior intestine mainly at 3 dpe. Intestinal endopeptidase, exopeptidase, and the proteins involved in antioxidant defense processes were upregulated upon Y. ruckeri exposure, while those involved in lipid metabolism, actin filament, and translation processes were downregulated. It is expected that this new information will exemplify further how Y. ruckeri influences the intestinal mucosa proteome of fish. The predicted lysosomal, metabolic, and focal adhesion pathways might be useful in understanding the gut defense mechanisms of rainbow trout and further research work in this direction. Utilization of these results may improve approaches for selection of disease-resistant rainbow trout broodstock and evaluation of prevention opportunities. Finally, the relationship among the intestinal bacterial pathogen and host protein responses will be explored to improve nutrient uptake, fish performance, and vaccine efficiency. Further studies are needed to address the functions of these intestinal proteins in immune response and protection against bacterial infection for both disease control and improved performance.
Availability of data and materials
Shotgun proteomics data generated during the current study have been deposited in the ProteomeXchange Consortium via the PRIDE partner repository  with the dataset identifiers PXD011087.
enteric redmouth disease
false discovery rate
information dependent data acquisition
liquid chromatography–mass spectrometry
- Micro LC ESI QTOF MS/MS:
microflow high-performance liquid chromatography electrospray-ionisation quadrupole time-of-flight tandem mass spectrometry
principal component analysis
sequential window acquisition of all theoretical spectra
triple quadrupole time of flight
Horne MT, Barnes AC (1999) Enteric redmouth disease (Yersinia ruckeri). In: Woo PTK, Bruno DW (eds) Fish diseases and disorders. Viral, bacterial and fungal infections. CABI Publishing, Wallingford, pp 445–477
Kumar G, Menanteau-Ledouble S, Saleh M, El-Matbouli M (2015) Yersinia ruckeri, the causative agent of enteric redmouth disease in fish. Vet Res 46:103
Zorriehzahra MJ, Adel M, Torabi Delshad S (2017) Enteric redmouth disease: past, present and future: a review. Iran J Fish Sci 16:1135–1156
Busch RA, Lingg AJ (1975) Establishment of an asymptomatic carrier state infection of enteric redmouth disease in rainbow trout (Salmo gairdneri). J Fish Res Board Can 32:2429–2432
Tobback E, Hermans K, Decostere A, Van den Broeck W, Haesebrouck F, Chiers K (2010) Interactions of virulent and avirulent Yersinia ruckeri strains with isolated gill arches and intestinal explants of rainbow trout Oncorhynchus mykiss. Dis Aquat Org 90:175–179
Ohtani M, Villumsen KR, Strøm HK, Raida MK (2014) 3D Visualization of the initial Yersinia ruckeri infection route in rainbow trout (Oncorhynchus mykiss) by optical projection tomography. PLoS One 9:e89672
Davies RL (1991) Clonal analysis of Yersinia ruckeri based on biotypes, serotypes and outer membrane protein types. J Fish Dis 14:221–228
Romalde JL, Margariños B, Barja JL, Toranzo AE (1993) Antigenic and molecular characterization of Yersinia ruckeri. Proposal for a new intraspecies classification. Syst Appl Microbiol 16:411–419
Austin DA, Robertson PAW, Austin B (2003) Recovery of a new biogroup of Yersinia ruckeri from diseased rainbow trout (Oncorhynchus mykiss, Walbaum). Syst Appl Microbiol 26:127–131
Fouz B, Zarza C, Amaro C (2006) First description of non-motile Yersinia ruckeri serovar I strains causing disease in rainbow trout, Oncorhynchus mykiss (Walbaum), cultured in Spain. J Fish Dis 29:339–346
Arias CR, Olivares-Fuster O, Hayden K, Shoemaker CA, Grizzle JM, Klesius PH (2007) First report of Yersinia ruckeri biotype 2 in the U.S.A. J Aquat Anim Health 19:35–40
Strom-Bestor M, Mustamaki N, Heinikainen S, Hirvela-Koski V, Verner-Jeffreys D, Wiklund T (2010) Introduction of Yersinia ruckeri biotype 2 into Finnish fish farms. Aquaculture 308:1–5
Axelsson M, Fritsche R (1991) Effects of exercise, hypoxia and feeding on the gastrointestinal blood flow in the Atlantic cod Gadus morhua. J Exp Biol 158:181–198
Buddington RK, Krogdahl A, Bakke-Mckellep AM (1997) The intestines of carnivorous fish: structure and functions and the relations with diet. Acta Physiol Scand Suppl 638:67–80
Seth H, Axelsson M, Farrell AP (2011) The circulation and metabolism of the gastrointestinal tract. In: Grosell M, Farrell AP, Brauner CJ (eds) Fish physiology: the multifunctional gut of fish. The Multifunctional Gut of Fish Academic Press, Amsterdam, pp 351–393
Wulff T, Petersen J, Nørrelykke MR, Jessen F, Nielsen HH (2012) Proteome analysis of pyloric ceca: a methodology for fish feed development? J Agric Food Chem 60:8457–8464
Baumgarner BL, Bharadwaj AS, Inerowicz D, Goodman AS, Brown PB (2013) Proteomic analysis of rainbow trout (Oncorhynchus mykiss) intestinal epithelia: physiological acclimation to short-term starvation. Comp Biochem Physiol Part D Genomics Proteomics 8:58–64
Ingerslev HC, Strube ML, Lv Jørgensen, Dalsgaard I, Boye M, Madsen L (2014) Diet type dictates the gut microbiota and the immune response against Yersinia ruckeri in rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol 40:624–633
Na-Phatthalung P, Teles M, Voravuthikunchai SP, Tort L, Fierro-Castro C (2018) Immune-related gene expression and physiological responses in rainbow trout (Oncorhynchus mykiss) after intraperitoneal administration of Rhodomyrtus tomentosa leaf extract: a potent phytoimmunostimulant. Fish Shellfish Immunol 77:429–437
Salinas I, Parra D (2015) Fish mucosal immunity: intestine. In: Beck B, Peatman E (eds) Mucosal health in aquaculture. Academic Press, San Diego, pp 135–170
Tobback E, Decostere A, Hermans K, Ryckaert J, Duchateau L, Haesebrouck F, Chiers K (2009) Route of entry and tissue distribution of Yersinia ruckeri in experimentally infected rainbow trout Oncorhynchus mykiss. Dis Aquat Org 84:219–228
Wilson JM, Castro LFC (2011) Morphological diversity of the gastrointestinal tract in fishes. In: Grosell MA, Farrell P, Brauner CJ (eds) The multifunctional gut of fish. Academic Press, Amsterdam, pp 1–155
Martin SA, Dehler CE, Król E (2016) Transcriptomic responses in the fish intestine. Dev Comp Immunol 64:103–117
Acar Ü, Parrino V, Kesbiç OS, Lo Paro G, Saoca C, Abbate F, Yılmaz S, Fazio F (2018) Effects of different levels of pomegranate seed oil on some blood parameters and disease resistance against Yersinia ruckeri in rainbow trout. Front Physiol 9:596
Parrino V, Kesbiç OS, Acar Ü, Fazio F (2019) Hot pepper (Capsicum sp.) oil and its effects on growth performance and blood parameters in rainbow trout (Oncorhynchus mykiss). Nat Prod Res. https://doi.org/10.1080/14786419.2018.1550769 (in press)
Davey GC, Calduch-Giner JA, Houeix B, Talbot A, Sitja-Bobadilla A, Prunet P, Perez-Sanchez J, Cairns MT (2011) Molecular profiling of the gilthead sea bream (Sparus aurata L.) response to chronic exposure to the myxosporean parasite Enteromyxum leei. Mol Immunol 48:2102–2112
Evenhuis JP, Cleveland BM (2012) Modulation of rainbow trout (Oncorhynchus mykiss) intestinal immune gene expression following bacterial challenge. Fish Shellfish Immunol 146:8–17
Li C, Zhang Y, Wang R, Lu J, Nandi S, Mohanty S, Terhune J, Liu Z, Peatman E (2012) RNA-seq analysis of mucosal immune responses reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus. Fish Shellfish Immunol 32:816–827
Xia JH, Liu P, Liu F, Lin G, Sun F, Tu R, Yue GH (2013) Analysis of stress responsive transcriptome in the intestine of Asian seabass (Lates calcarifer) using RNA-seq. DNA Res 20:449
Kumar G, Hummel K, Noebauer K, Razzazi-Fazeli E, El-Matbouli M (2018) Proteome analysis reveals a role of rainbow trout lymphoid organs during Yersinia ruckeri infection process. Sci Rep 8:13998
Gibello A, Blanco MM, Moreno MA, Cutuli MT, Domenech A, Domínguez L, Fernández-Garayzábal JF (1999) Development of a PCR assay for detection of Yersinia ruckeri in tissues of inoculated and naturally infected trout. Appl Environ Microbiol 65:346–350
Erde J, Loo RR, Loo J (2014) Enhanced FASP (eFASP) to increase proteome coverage and sample recovery for quantitative proteomic experiments. J Proteome Res 13:1885–1895
R Core Team R (2015) A language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria. https://www.R-project.org/ Accessed 15 Mar 2019
Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B 57:289–300
Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, Thomas PD (2017) PANTHER version 11: expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements. Nucleic Acids Res 45:D183–D189
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13:2498–2504
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman WH, Pagès F, Trajanoski Z, Galon J (2009) ClueGO: a Cytoscape plug-into decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 25:1091–1093
Kumar G, Abd-Elfattah A, El-Matbouli M (2014) Differential modulation of host genes in the kidney of brown trout Salmo trutta during sporogenesis of Tetracapsuloides bryosalmonae (Myxozoa). Vet Res 45:101
Raida MK, Buchmann K (2008) Development of adaptive immunity in rainbow trout, Oncorhynchus mykiss (Walbaum) surviving an infection with Yersinia ruckeri. Fish Shellfish Immunol 5:533–541
Di G, Li H, Zhang C, Zhao Y, Zhou C, Naeem S, Li L, Kong X (2017) Label-free proteomic analysis of intestinal mucosa proteins in common carp (Cyprinus carpio) infected with Aeromonas hydrophila. Fish Shellfish Immunol 66:11–25
Wiens GD, Vallejo RL (2010) Temporal and pathogen-load dependent changes in rainbow trout (Oncorhynchus mykiss) immune response traits following challenge with biotype 2 Yersinia ruckeri. Fish Shellfish Immunol 29:639–647
Raida MK, Holten-Andersen L, Buchmann K (2011) Association between Yersinia ruckeri infection, cytokine expression and survival in rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol 30:1257–1264
Sternlicht MD, Werb Z (2001) How matrix metalloproteinases regulate cell behavior. Annu Rev Cell Dev Biol 17:463–516
Feng T, Zhang H, Liu H, Zhou Z, Niu D, Wong L, Kucuktas H, Liu X, Peatman E, Liu Z (2011) Molecular characterization and expression analysis of the channel catfish cathepsin D genes. Fish Shellfish Immunol 31:164–169
Dong ZD, Zhang J, Ji XS, Zhou FN, Fu Y, Chen W, Zeng YQ, Li TM, Wang H (2012) Molecular cloning, characterization and expression of cathepsin D from grass carp (Ctenopharyngodon idella). Fish Shellfish Immunol 33:1207–1214
Grinde B, Jollès J, Jollès P (1988) Purification and characterization of two lysozymes from rainbow trout (Salmo gairdneri). Eur J Biochem 173:269–273
Laing KJ, Holland J, Bonilla S, Cunningham C, Secombes CJ (2001) Cloning and sequencing of caspase 6 in rainbow trout, Oncorhynchus mykiss, and analysis of its expression under conditions known to induce apoptosis. Dev Comp Immunol 25:303–312
Fu S, Ding M, Yang Y, Kong J, Li Y, Guo Z, Wang A, Ye J (2018) Molecular cloning, characterization and expression analysis of caspase-6 in puffer fish (Takifugu obscurus). Aquaculture 490:311–320
Gerwick L, Reynolds WS, Bayne CJ (2000) A precerebellin-like protein is part of the acute phase response in rainbow trout, Oncorhynchus mykiss. Dev Comp Immunol 24:597–607
Raida MK, Buchmann K (2009) Innate immune response in rainbow trout (Oncorhynchus mykiss) against primary and secondary infections with Yersinia ruckeri O1. Dev Comp Immunol 33:35–45
Machaca K (2010) Ca2+ signaling, genes and the cell cycle. Cell Calcium 48:243–250
Takahashi H, Yamaguchi M (1995) Increase of (Ca(2+)–Mg2+)-ATPase activity in hepatic plasma membranes of rats administered orally calcium: the endogenous role of regucalcin. Mol Cell Biochem 144:1–6
Castro A, Becerra M, Manso MJ, Anadón R (2006) Calretinin immunoreactivity in the brain of the zebrafish, Danio rerio: distribution and comparison with some neuropeptides and neurotransmitter-synthesizing enzymes. II. Midbrain, hindbrain, and rostral spinal cord. J Comp Neurol 494:792–814
Porteros A, Arévalo R, Weruaga E, Crespo C, Briñón JG, Alonso JR, Aijon J (1997) Calretinin immunoreactivity in the developing olfactory system of the rainbow trout. Brain Res Dev Brain Res 100:101–109
Verleih M, Rebl A, Köllner B, Korytář T, Anders E, Wimmers K, Goldammer T (2012) Comparative molecular characterization of the regucalcin (RGN) gene in rainbow trout (Oncorhynchus mykiss) and maraena whitefish (Coregonus marena). Mol Biol Rep 39:4291–4300
Farber SA, De Rose RA, Olson ES, Halpern ME (2003) The zebrafish annexin gene family. Genome Res 13:1082–1096
Santamaria-Kisiel L, Rintala-Dempsey AC, Shaw GS (2006) Calcium-dependent and -independent interactions of the S100 protein family. Biochem J 396:201–214
Kumar G, Gotesman M, El-Matbouli M (2015) Interaction of Tetracapsuloides bryosalmonae, the causative agent of proliferative kidney disease, with host proteins in the kidney of Salmo trutta. Parasitol Res 114:1721–1727
Salem M, Nath J, Rexroad CE, Killefer J, Yao J (2005) Identification and molecular characterization of the rainbow trout calpains (Capn1 and Capn2): their expression in muscle wasting during starvation. Comp Biochem Physiol B: Biochem Mol Biol 140:63–71
Skov J, Kania PW, Holten-Andersen L, Fouz B, Buchmann K (2012) Immunomodulatory effects of dietary β-1,3-glucan from Euglena gracilis in rainbow trout (Oncorhynchus mykiss) immersion vaccinated against Yersinia ruckeri. Fish Shellfish Immunol 33:111–120
Vizcaíno JA, Csordas A, del Toro N, Dianes JA, Griss J, Lavidas I, Mayer G, Perez-Riverol Y, Reisinger F, Ternent T, Xu QW, Wang R, Hermjakob H (2016) 2016 update of the PRIDE database and related tools. Nucleic Acids Res 44:D447–D456
This research was supported using resources of the VetCore Facility (Proteomics) of the University of Veterinary Medicine Vienna. We are thankful to Dr. Maike Ahrens for helping in R statistical analysis. Yersinia ruckeri strain CSF007-82 was kindly provided by Dr. Timothy J. Welch, National Center for Cool and Cold Water Aquaculture, Kearneysville, USA.
This study was funded by the Austrian Science Fund (FWF) Project No. P 27489-B22 to GK.
Ethics approval and consent to participate and to publication
This study was approved by the institutional ethics committee of the University of Veterinary Medicine Vienna and the national authority, according to §26 of the Austrian Law for Animal Experiments, Tierversuchsgesetz 2012 under approval number BMWFW-68.205/0041-WF/V/3b/2015.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.