- Open Access
SalmoSim: the development of a three-compartment in vitro simulator of the Atlantic salmon GI tract and associated microbial communities
Microbiome volume 9, Article number: 179 (2021)
The aquaculture sector now accounts for almost 50% of all fish for human consumption and is anticipated to provide 62% by 2030. Innovative strategies are being sought to improve fish feeds and feed additives to enhance fish performance, welfare, and the environmental sustainability of the aquaculture industry. There is still a lack of knowledge surrounding the importance and functionality of the teleost gut microbiome in fish nutrition. In vitro gut model systems might prove a valuable tool to study the effect of feed, and additives, on the host’s microbial communities. Several in vitro gut models targeted at monogastric vertebrates are now in operation. Here, we report the development of an Atlantic salmon gut model, SalmoSim, to simulate three gut compartments (stomach, pyloric caecum, and midgut) and associated microbial communities.
The gut model was established in a series of linked bioreactors seeded with biological material derived from farmed adult marine-phase salmon. We first aimed to achieve a stable microbiome composition representative of founding microbial communities derived from Atlantic salmon. Then, in biological triplicate, the response of the in vitro system to two distinct dietary formulations (fishmeal and fishmeal free) was compared to a parallel in vivo trial over 40 days. Metabarcoding based on 16S rDNA sequencing qPCR, ammoniacal nitrogen, and volatile fatty acid measurements were undertaken to survey the microbial community dynamics and function. SalmoSim microbiomes were indistinguishable (p = 0.230) from their founding inocula at 20 days and the most abundant genera (e.g., Psycrobacter, Staphylococcus, Pseudomonas) proliferated within SalmoSim (OTUs accounting for 98% of all reads shared with founding communities). Real salmon and SalmoSim responded similarly to the introduction of novel feed, with majority of the taxa (96% Salmon, 97% SalmoSim) unaffected, while a subset of taxa (e.g., a small fraction of Psychrobacter) was differentially affected across both systems. Consistent with a low impact of the novel feed on microbial fermentative activity, volatile fatty acid profiles were not significantly different in SalmoSim pre- and post-feed switch.
By establishing stable and representative salmon gut communities, this study represents an important step in the development of an in vitro gut system as a tool for the improvement of fish nutrition and welfare. The steps of the system development described in this paper can be used as guidelines to develop various other systems representing other fish species. These systems, including SalmoSim, aim to be utilised as a prescreening tool for new feed ingredients and additives, as well as being used to study antimicrobial resistance and transfer and fundamental ecological processes that underpin microbiome dynamics and assembly.
In the last 50 years, per capita fish consumption has almost doubled from 10 kg in the 1960s to over 19 kg in 2012 . This increase in the demand for fish protein has put wild fish stocks under pressure. The aquaculture sector now accounts for almost 50% of all fish for human consumption and is anticipated to provide 62% by 2030 . Freshwater fish species currently dominate aquaculture fish production, such as freshwater carps and cyprinids, which constitute over 53% of total fish production by mass, with tilapia and other cichlids (11.0%) . However, when reported by value, Atlantic salmon (Salmo salar) ranked highest . While husbandry and disease control have advanced to improve fish health and welfare, improved feeds and their additives are the fundamental means of enhancing fish performance in aquaculture . To improve aquaculture environmental and financial sustainability, the ratio of the marine origin components (Fishmeal (FM) and Fish Oil (FO)) within feeds has reduced considerably. For example, in Norway the ratio of the marine origin components within farmed salmon feed reduced from around 90% in 1990 to 30% in 2013 .
Recent studies suggest that teleost microbiome plays a vital role in fish’s health and performance [7, 8], and that alternative feeds, such as ones containing non-marine dietary ingredients, can result in poor fish growth, altered gut health alongside a modified fish gut microbial community composition and activity [9,10,11]. For instance, Atlantic salmon feed supplementation with dietary soybean protein concentrate can induce intestinal disorder . Concomitant alterations in gut microbiota can result in the undesirable fermentation of various feed components [12, 13]. In view of all this, considerable interest lies around the development of novel ingredients and additives to enhance the performance of many species of farmed fish and their associated microbes.
To study the impact of novel feed ingredients on gut microbial communities (e.g., Gajardo et al., 2017), as well as the supplements (e.g., prebiotics, probiotics) tailored to modify microbial community diversity and function (e.g., Gupta et al. ), in vivo trials are widely deployed in aquaculture. Although physiologically relevant, in vivo trials have several scientific, ethical, and practical disadvantages. In salmonids, for example, gut sampling is terminal, preventing the generation of time series data from individual animals/microbial communities. Furthermore, microbial impacts on feed ingredients cannot be subtractively isolated from host enzymatic/cellular activity. From an ethical perspective, in vitro models offer the opportunity to reduce harm via the replacement of in vivo models . Practically, in vivo, testing of novel feed ingredients is both time consuming and costly. A reliable in vitro model could offer advantages in this respect. To the best of our knowledge, there is only one other gut system in place simulating a generalised teleost gut, (‘fish-gut-on-chip’ ). The ‘fish-gut-on-chip’ system exploits microfluidic technology and is based on the reconstruction of the rainbow trout's intestinal barrier by culturing only intestinal cell lines in an artificial microenvironment and currently does not involve microbial communities isolated from the fish’s gut.
Prior to deploying an in vitro gut microbiome simulator to perform biological experiments, several criteria must be met. Firstly, steady-state microbial communities need to be established prior to the experimental procedure to ensure that results due to experimental treatments are not confounded with bacterial adaptation to the in vitro environment . Secondly, physicochemical conditions within the artificial gut simulator and the gut of the target species should be similar. Thirdly, the bacterial communities need to be gut compartment–specific and representative of (if not identical to) the in vivo situation . Finally, the in vitro gut simulator should be validated against a parallel in vivo experiment, to establish the degree to which the results from the experimental protocol within the artificial gut are generalisable to the in vivo situation . Towards this end, several molecular techniques can then be deployed to analyse microbial populations within the gut. Multiplex quantitative PCR (qPCR) coupled with taxon-specific primers can rapidly detect and quantify the bacterial consortia within a large population . Whilst shotgun metagenomics and amplicon sequencing approaches can provide a detailed taxonomic assessment of the microbial composition of the gut, they may be less useful for day-to-day monitoring of specific taxa .
In view of the above criteria, the aim of the current study is to develop a synthetic, continuous teleost gut microbial fermentation simulator. Teleosts are the largest and most diverse group of living vertebrates . Due to the wide variety of habitats teleosts exploit and foods that they consume (ranging from bottom-living seaweeds and plankton organisms to actively swimming animals), fish have a wide array of gut morphologies . Teleost stomachs can be classified into four general configurations: a straight stomach with enlarged lumen (e.g., northern pike (Esox lucius)), Y-shaped stomach (e.g., Nile tilapia (Oreochromis niloticus)), the absence of stomach (e.g., common carp (Cyprinus carpio)) and a U-shaped sack-like stomach with enlarged lumen, such as is found in Atlantic salmon, receiving food from the fish via the oesophagus [24, 25]. If an acidic stomach is present, it is responsible for initial unspecific digestion of incoming food by secretion of hydrochloric acid and endopeptidase pepsin from gastric mucosa glands in the stomach lining . From the stomach, the food transport to the midgut is controlled by a muscular sphincter, the pylorus, which is developed to varying degrees in different species for unknown reasons, in some species even being absent (e.g., common carp in which the midgut attaches directly to the oesophagus) . The pylorus has proven to be important in enzymatic breakdown of ingested macromolecules via the secretion of proteases (e.g., trypsinogen), glucosidases (e.g., α-amylase), and lipases (e.g., bile salt-dependent lipase))  and to utilize the counter-current multiplication in generating osmoregulatory mechanisms for absorption of glucose, amino acids, dipeptides, and medium chain fatty acids [29,30,31,32]. The gastrointestinal tract ends with mid and distal intestines that further digest and absorb nutrients . Thus, considering that Atlantic salmon has all three gut compartments, and is a leading species in worldwide aquaculture by value, it and its gut system, were chosen as a model for simulation in vitro in this study. Nonetheless, the principles and methodology described in this study could be applied to the development and construction of other systems representing different fish species as well.
Our experimental gut system simulates the stomach, the pyloric caeca, and the midgut regions of the gastrointestinal tract of the generalised marine lifecycle stage of the farmed Atlantic salmon. In this context, we first aimed to establish a stable gut community, representative of the salmon gut communities used to found it. Secondly, we validated the system as a potential means of testing the impact of feeds on salmon gut microbial communities by comparing the performance and response of the in vitro simulator during a feed trial with a parallel in vivo modulation of the gut microbial community in a cohort of marine-phase Atlantic salmon.
Experimental set-up and sample collection in an aquaculture setting
The Atlantic salmon (Salmo salar) in vivo feed trial was performed by MOWI ASA at their research site in Averøy, Norway. Prior to commencement of the feed trial, salmon were fed on a Fishmeal diet (FMD) until they reached circa 750 g in mass. Fish were separated into 5 × 5 m marine pens (150 randomly distributed fish per pen) in a 4 × 4 modular design. Four pens were randomly assigned to each of the trial diets. This study focused on eight pens housing fish fed on FMD and Fishmeal-free diet (FMF) (Table 1, Fig. 1D). The feed trial was conducted over 5 months (November 2017–March 2018). For in vivo samples recovered from actual salmon, three randomly selected fish were collected at the end of the feed trial for two different feeds (N = 3 fish/feed: Fish 1, 2, and 3 for FMD and Fish 4, 5, and 6 for FMF) and sacrificed by MOWI employees (Fig. 1E). After, 1 cm in length samples from three salmon gut compartments were collected (stomach (N = 3/feed) pyloric caeca (N = 3/feed) and midgut (N = 3/feed) (approximately 20 cm from the vent)), placed into 1.5-mL cryovials, and kept on ice before long-term storage in – 80 °C conditions. For in vitro initial inoculum samples (the founding community for SalmoSim runs), a further three fish fed on FMD were sacrificed (Fish 7, 8, and 9), and 5 cm in length samples from three distinct gut compartments were collected (stomach (N = 3), pyloric caecum (N = 3) and midgut (N = 3)), transferred to 15-mL Falcon tubes and kept on ice before long-term storage in – 80 °C conditions (Fig. 1E). Details of samples collected from farmed Atlantic salmon have been described previously .
Establishment of stable, representative gut communities in SalmoSim and subsequent feed trial
Physicochemical conditions within Atlantic salmon gastrointestinal tract and microbiome sampling
Physicochemical conditions (temperature, pH, dissolved oxygen) were directly measured in adult Atlantic salmon from MOWI salmon farm in Loch Linnhe, Scotland (Supplementary Figure 1A-C). Bacterial inoculums were prepared for the in vitro trial from the different gut compartments sampled from individual fish (three biological replicates, three gut compartments per fish—stomach, pyloric caecum, and midgut) collected at the start of the in vivo feed trial in Averøy, Norway. Prior to SalmoSim inoculation, inoculums that were stored in 15-mL Falcon tubes in − 80°C freezer were dissolved in 1 mL of autoclaved 35 g/L Instant Ocean® Sea Salt solution. Distinct individual fish collected in Averøy formed the founder community for each distinct replicate of the in vitro trial (i.e., a true biological replicate (Fig. 1E)).
In vitro system ‘feed’ preparation
In vitro system feed media was prepared by combining the following for a total of 2 L: 35 g/L of Instant Ocean® Sea Salt, 10 g/L of the FMD or FMF used in the MOWI feed trial (Table 1: the concentration of feed added was optimized to run through the system without clogging the narrow tubing), 1 g/L freeze-dried mucous collected from the pyloric caecum of Scottish marine-phase Atlantic salmon and 2 L of deionised water. This feed was then autoclave-sterilised, followed by sieving of the bulky flocculate, and finally subjected to the second round of autoclaving. This subjection of the feed to two rounds of autoclave should have degraded the vast majority of DNA, thus ensuring that feed microbiome will not influence the microbiome within the SalmoSim system.
In vitro system preparation
Three 500-mL Applikon Mini Bioreactors (one for each gut compartment: stomach, pyloric caecum, and midgut) were filled with four 1 cm3 cubes made from sterile aquarium sponge filters used as a surface for biofilm formation, assembled by attaching appropriate tubing and probes (redox, temperature, and dissolved oxygen; Fig. 1A), and autoclaved. Bioreactor preparation was followed by attachment of reactor vessels to the Applikon electronic control module, connection of feed and acid and base bottles (0.01 M hydrochloric acid and 0.01 M sodium hydroxide solutions filtered through a 0.22-μm polyethersulfone membrane filter unit (Millipore, USA)). Nitrogen gas was periodically bubbled through each vessel to maintain anaerobic conditions and dissolved oxygen continually monitored. The bioreactors were then allowed to fill with 400 mL of feed media. Once the system was set up, media transfer, gas flow, and acid/base addition occurred for 24 h in sterile conditions (without microorganisms present in the system) in order to stabilise the temperature, pH, and oxygen concentration with respect to levels measured from adult salmon.
Initial pregrowth period during in vitro trial
In order to allow bacterial communities to proliferate in the in vitro environment without washing through the system, the microbial populations within the inoculum from real salmon were pre-grown inside the SalmoSim system for 4 days. During this phase, the system was filled with FMD media preparation and inoculum, and no media transfer occurred.
Performing feed trial within SalmoSim system
After the initial pregrowth period, each validation experiment was run for 20 days while supplying the SalmoSim system with FMD. After 20 days, SalmoSim was run for 20 additional days while supplying FMF food. During the full 44-day experiment (4-day pregrowth period, 20-day system fed on FMD, and 20-day system fed on FMF) physiochemical conditions within three SalmoSim gut compartments were kept similar to the values measured in real salmon: temperature inside the reactor vessels were maintained at 12 °C, dissolved oxygen contents were kept at 0% by daily flushing with N2 gas for 20 min, and pH was kept stable in each bioreactor by the addition of 0.01 M NaOH and 0.01 M HCl (stomach pH 4.0, pyloric caecum pH 7.0, and mid intestine pH 7.6). During this experiment (apart from the pregrowth period), the transfer rate of slurry between reactor vessels was 238 mL per day as described in vivo [34, 35]. Finally, 1 mL of filtered salmon bile and 0.5 mL of sterile 5% mucous solution (both collected from Scottish marine-phase Atlantic salmon) were added daily to the reactor, simulating the pyloric caecum compartment. The schematic representation of the SalmoSim system is visualised in Fig. 1A, and full feed trial within SalmoSim is visually summarised in Fig. 1B and C.
Sampling was performed in several steps. First, samples from initial inoculums were collected for each gut compartment. Once SalmoSim main experiment started, the sampling from each bioreactor vessel was performed every second day throughout the 40-day run period (20 samplings in total). The SalmoSim samplings entailed collecting 30 mL of the bioreactor contents (10% of the total bioreactor volume—the maximum volume of sampling without disturbing microbial community structure) into 50 mL Falcon tubes, centrifuging them for 10 min at 5000 rpm speed, and freezing the pellets at – 20 °C for storage. The pellets were frozen to perform DNA extraction all together with the same kit (to prevent batch effect), and supernatant was used for VFA analysis.
Nitrogen metabolism within the SalmoSim system
At each sampling point, the protein concentration in each chamber of the system was measured using Thermo Scientific™ Pierce™ BCA Protein Assay Kit (Thermo Fisher Scientific, USA) and the ammonia concentration using Sigma-Aldrich® Ammonia Assay Kit (Sigma-Aldrich, USA) to assay the bacterial community activity. Both methods were performed according to manufacturer’s protocol by using the Jenway 6305 UV/Visible Spectrophotometer (Jenway, USA). The same samples were used for both of these analyses immediately after sampling (no freezing or intermediate steps required).
Volatile fatty acid (VFA) production in SalmoSim
The last two time points for each diet were selected from the SalmoSim system (for all three gut compartments) for VFA analysis—18 and 20 for FMD and time points 38 and 40 for FMF, respectively—to ensure that the bacterial communities had enough time to adapt to the SalmoSim system (for FMD time points) and the feed change (for FMF time points). During runs, 1 mL of supernatant from SalmoSim bioreactors was frozen in – 80 °C which was then used for VFA extraction. The protocol involved combining 1 mL of supernatant with 400 μL of sterile phosphate-buffered saline (PBS) solution (Sigma Aldrich, USA) and vortexing the mixture for 1 min. The sample was then centrifuged at 16,000 g for 30 min, followed by two rounds of supernatant removal, before centrifuging at 16,000 g for 30 min. Finally, the supernatant was then filtered through 0.2 μm Costar SpinX centrifuge tube filters (Corning, USA) at 15,000 g for 2 min until clear. The extracted VFAs were sent for gas chromatographic analysis at MS-Omics (Denmark). In order to determine if the VFA concentrations were statistically different between SalmoSim fed on FMD and FMF diets, measured VFA values dataset were subjected to statistical analysis using linear mixed effect models (See Supplementary methods 2). Results are shown in supplementary Figure 9.
In vivo phenotypic fish performance fed on two different feeds
Phenotypic performance data (fish length, weight, gutted weight, carcass yield, gonad, and liver weights) was collected and provided at the end of the in vivo feed experiment by MOWI. The differences between each feed group (n = 32 fish per feed) for each phenotype were visualised, and statistical analysis was applied (independent two-sample t test) to identify statistically significant differences between the two feed groups.
Measuring bacterial population dynamics in SalmoSim
Genomic DNA extraction
The DNA extraction protocol followed was previously described by . In short, samples were subjected to a bead-beating step for 60 s by combining samples with MP Biomedicals™ 1/4" CERAMIC SPHERE (Thermo Fisher Scientific, USA) and Lysing Matrix A Bulk (MP Biomedicals, USA). Later, DNA was extracted by using the QIAamp® DNA Stool kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s protocol .
NGS library preparation and sequencing
In the first instance, microbial population dynamics in SalmoSim were measured in near real-time using a set of qPCR primers including published and custom sequences to enable the stability of the system to be monitored (See supplementary Methods 1 and data Supplementary Figure 4). Subsequently, 16S rRNA sequencing was deployed to provide a fuller picture community dynamics. The commonly used 16S ribosomal hypervariable region 4 primers were shown to cross-amplify Salmo salar 12S ribosomal gene sequences [33, 37] and hence were not used in this study. Rather, amplification of the 16S V1 hypervariable region was adopted as an alternate approach . Amplification of 16S V1 hypervariable region from diluted DNA samples was achieved using redundant tagged barcode 27F and 338R at final concentration of 1 pM of each primer. Primer sequences are summarised in Supplementary Table 2. First-round PCR was performed in triplicate on each sample, and reaction conditions were 95 °C for 10 min, followed by 25 cycles at 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s, followed by a final elongation step of 72 °C for 10 min. After the triplicate reactions were pulled together into one, their concentration was measured using a Qubit® fluorometer (Thermo Fisher Scientific, USA), and all of them were diluted to 5 ng/μL using Microbial DNA-Free Water (Qiagen, Valencia, CA, USA). The second-round PCR, which enabled the addition of the external multiplex identifiers (barcodes), involved only six cycles with otherwise identical reaction conditions to the first. The detailed composition of second-round PCR primers is summarised in Supplementary Table 3. This was followed by the DNA clean-up using Agencourt AMPure XP beads (Beckman Coulter, USA) according to the manufacturers' protocol. The cleaned DNA was then gel-purified by using the QIAquick Gel Extraction Kit (Qiagen, Valencia, CA, USA) and then quantified by using Qubit® (Thermo Fisher Scientific, USA). All the PCR products were pulled together at 10 nM concentration and sent for sequencing using HiSeq 2500.
Bioinformatic analysis of 16S rRNA gene amplicon sequencing data
Sequence analysis was performed with our bioinformatic pipeline as described previously with slight modifications . First, quality filtering and trimming (> Q30 quality score) was performed on all the reads of the 16S rRNA V1 hypervariable region using Sickle version V1.2 software . Second, read error correction was performed using the BayesHammer module within SPAdes V2.5.0 software to obtain high-quality assemblies . Third, paired-end reads were merged (overlap length 50 bp) using PANDAseq v2.11 software with a simple_Bayesian read merging algorithm [41, 42]. After overlapping paired-end reads, merged reads were dereplicated, sorted, and filtered chimaeras using the GOLD SILVA reference dataset , and singletons were removed by using the VSEARCH version 2.3.4 tool . Merged paired-end filtered reads were clustered in operational taxonomic units (OTUs) using the VSEARCH software at 97% identity followed by a decontamination step by mapping OTUs against the host (Salmo salar) reference genome (available on NCBI) DNA using the BWA aligner implemented in the DeconSeq v0.4.3 tool . Taxonomic assignment of OTUs was achieved using the Naïve Bayesian Classifier  implemented in the QIIME 2 platform using the SILVA 132 database [47, 48]. Phylogenetic trees of OTUs were generated using FastTree software after using MAFFT for multiple sequence alignment . The resultant OTU table was converted to a biological observation matrix (BIOM) format for the post-OTU statistical analysis .
Post-OTU statistical analysis: diversity metrics and community structure and composition analysis
All statistical analysis of the OTU table was performed by using R v 4.0.1 and RStudio v 1.3.959 . Alpha diversity analysis was performed using Rhea pipeline , supplemented by microbiomeSeq , and PhyloSeq  for ANOVA and visualisation steps. Two alpha diversity metrics were calculated: microbial richness (estimated number of observed OTUs) and Shannon diversity (an estimate of balance of the community using the effective Shannon index [55, 56]. Before calculating effective microbial richness, proportional filtering was performed at a relative abundance of 0.25% in each community to minimise the inflation in microbial richness caused by the very low abundant OTUs. Afterwards, a one-way ANOVA of diversity between groups was calculated with the p value threshold for significance (p value < 0.05) represented using boxplots.
To investigate the effect of time on the bacterial community stability, beta diversity analysis was performed using different phylogenetical distance metrics to assay phylogenetic similarities between samples (weighted, generalised, and unweighted UniFrac). To compare communities isolated from various sources (SalmoSim, inoculum, and real salmon), only samples fed on FMD were included as initial inoculum were collected from fish fed on FMD alone. Furthermore, only SalmoSim samples from the last 3 time points fed on FDM were selected as they are considered stable time points (once bacterial communities had over 2 weeks to adapt and grow within the SalmoSim system). In short, the resulting dataset contained real salmon samples fed on FMD, all inoculum samples, and stable SalmoSim time points fed on FMD (days 16, 18, and 20). This dataset was then subdivided into several smaller datasets that included OTUs, shared by various percentages of samples (60%, 50%, 40%, and 30% of samples), with the aim of minimising the impact of rare OTUs (low prevalence) on comparisons and only focusing on prevalent OTUs between samples (see details in Supplementary Table 4).
To analyse the response of microbes to the diet change (see Table 1 for feed formulation) in real salmon and SalmoSim, in addition to the full dataset (in vivo and in vitro samples); three different full dataset subsets were used to perform beta diversity analysis: samples from in vivo study, all samples from SalmoSim (all data points), and samples only from SalmoSim once it had achieved stability (the last 3 time points fed on FMD (days: 16, 18, and 20) and FMFD (days 36, 38, and 40). These datasets were used to compute ecological (Bray–Curtis and Jaccard) and phylogenetic (unweighted, weighted, and generalised UniFrac) distances with vegdist function from the vegan v2.4-2 package and GUniFrac function (generalised UniFrac) from the Rhea package [52, 57]. Both ecological and phylogenetical distances were then visualised in two dimensions by Multi-Dimensional Scaling (MDS) and nonmetric MDS (NMDS) . Finally, a permutational multivariate analysis of variance (PERMANOVA) by using both calculated distances was performed using adonis function to determine if the separation of selected groups was significant as a whole and in pairs . The full beta diversity workflow is summarised in Supplementary Methods 3.
To provide an overall visualisation of microbial composition across all samples, a Principal Coordinates Analysis (PCoA) was performed using the microbiomeSeq  package based on the phyloseq package  with Bray–Curtis dissimilarity measures calculated and visualised for four different subset datasets: the full dataset (real salmon, inoculum, and all SalmoSim samples) and three different subsets each containing only one of the free biological replicate samples from SalmoSim (Fish 1, 2, or 3), along with all real salmon and inoculum samples.
Differential abundance was calculated by using microbiomeSeq based on the DESeq2 package [53, 59]. The BIOM-generated OTU table was used as an input to calculate differentially abundant OTUs between selected groups based on the Negative Binomial (Gamma-Poisson) distribution.
Stabilisation of representative microbial communities within the SalmoSim system
Effective richness (Fig. 2A) indicates that within the stomach and midgut compartments, the initial inoculum contained the highest number of OTUs compared with later sampling time points from the SalmoSim system: in the stomach compartment, effective richness was statistically different between time point 0 (initial inoculum derived from salmon guts) and time points 18, 30, 36, and 38, and within the midgut compartment, the number of OTUs within the initial inoculum (time point 0) was statistically different from time points 2, 4, 6, 16, 34, 36, 38, and 40. However, within the pyloric caeca compartment, only one time point (time point 34) had a significantly different number of OTUs from the initial inoculum (time point 0).
Our results reveal that within the stomach compartment over time (including initial inoculum), the effective Shannon diversity was stable with a downwards but nonsignificant trend over the course of the experiment (Fig. 2B). A similar downwards trend was observed in the pyloric caecum, with significant differences between later time points, but no significant differences were noticed between the inoculum and SalmoSim. Within the midgut compartment, Shannon diversity was statistically lower than the inoculum (time point 0) over most time points (sampling days 2–40).
Taken together, diversity and richness metrics suggest some loss of microbial taxa as a result of transfer of salmon gut communities into SalmoSim in the pyloric caecum and midgut, but not in the stomach. Subsequently, richness and evenness were then stable over the time course of the experiment in the stomach and midgut compartments (some instabilities seen only between the initial inoculum and later time points), whilst much more instability within the alpha diversity metrics were detected in the pyloric caecum compartment.
To assess the compositional stability of the system, comparisons over time were undertaken with reference to pairwise beta diversity metrics. Significant differences in composition between time points represent instability in the system. Figure 3 visually summarises between–time point comparison of beta diversity metrics within the SalmoSim system across all replicates using generalised UniFrac (visual representations of the results using unweighted and weighted UniFrac are summarised in Supplementary Figure 3). Irrespective of the metric used, microbial community composition appeared to stabilise in all gut compartments by approximately 20 days, with little-observed impact of introducing the different feed at day 20. This trend was supported by our qPCR data, suggesting increasing stability over the course of the 40-day experiment (Supplementary Figure 4). Prevailing stability was also observed when each biological replicate’s individual gut compartment was examined separately (stomach, pyloric caecum, and midgut) (Supplementary Figure 5). Importantly, stabilisation over 20 days was a feature of two previous pilot runs of the system using unrelated marine salmon gut communities (Supplementary Figure 6).
Microbial identity and diversity compared between SalmoSim and salmon
To compare microbial identity and microbiome composition between salmon and SalmoSim sample sizes were first balanced by examining a reduced dataset that contained real salmon samples (three fish individuals used to inoculate and three others, all fed FMD) and stable SalmoSim time points fed on FMD (days 16, 18, and 20). Alpha diversity comparisons between inoculum, real salmon, and SalmoSim are shown in Fig. 2.
Compositional comparisons between different sample types (inoculum salmon, three other individuals and SalmoSim) were made using several pairwise beta diversity metrics (phylogenetic (unweighted, generalised, and weighted UniFrac and ecological distances (Bray–Curtis and Jaccard)) (Table 2). Ecological metrics could not distinguish between SalmoSim (days 16, 18, and 20) from the individual fish used to inoculate the system, suggesting that microbiome composition was very similar between salmon and SalmoSim (Table 2). In contrast, metrics that incorporate phylogenetic differences between taxa (i.e., Unifrac) did identify significant differences, indicating that there is variability between the inoculum and SalmoSim, but the taxa involved are closely related. To explore the impact of rare OTUs when accounting for observed differences between sample types, the dataset was partitioned and analysed. Partitioned datasets indicated that progressive removal of rare OTUs increased the compositional similarity between SalmoSim and the fish gut communities used to inoculate the system (Table 2). Inoculum samples (18 samples) generated 388 OTUs, and SalmoSim stable samples (54 samples in total: SalmoSim time points fed on FMD (days 16, 18, and 20) and SalmoSim time points fed on FMF (days 36, 38, and 40)) generated 508 OTUs. A total of 291 OTUs were present in both sample types. Importantly, the shared 291 OTUs accounted for > 97% of the total reads in inoculum samples and > 98% in stable SalmoSim samples (Table 2), again suggesting that the majority of abundant microbes in real salmon were maintained in the SalmoSim system.
Between real salmon that were not the direct source of inocula and SalmoSim, and between salmon used as inocula and other individual salmon, however, statistically significant differences were found in using all metrics regardless of inclusion of rare OTUs. These observations are consistent with inter-individual variation—SalmoSim and inoculum samples originated from the same individuals, while other salmon samples were, by necessity, collected from different individuals during the in vivo trial. Furthermore, while the number of OTUs between salmon not used to inoculate (192 OTUs) and inoculum salmon samples (388 OTUs) are different (Table 2), these non-inoculum salmon share 131 OTUs out of 192 OTUs with inoculums, and these 131 OTUs account for around 98% of the total reads. Thus, extra OTUs found only in inoculum salmon and not in others are relatively rare in abundance terms. Differences in OTU numbers and composition are not unexpected as a slightly larger amount of inoculum sample was collected (5 cm of intestine length vs 1 cm for other salmon).
Effect of changing diet on the microbiome of real salmon (in vivo) and SalmoSim (in vitro)
The impact of diet on the abundance of individual taxa
In response to the change of diet, the relative abundances of individual taxa in salmon vs SalmoSim also revealed some differences, as well as multiple similarities in response of the two systems (Fig. 4). In this respect, the abundance of the vast majority of OTUs (SalmoSim: 97%; Salmon: 95%; Fig. 4C) were unaffected by the change in feed; these included 161 OTUs shared by SalmoSim and the real salmon assayed. For OTUs whose individual abundance was impacted by feed across the two systems, only a single common OTU changed in the same way in both Salmon and SalmoSim (Fig. 4A). qPCR-based estimates of taxon abundance variation in response to diet (Supplementary Table 5) and corresponding data for the same taxa from 16S OTU profiles (Fig. 4D) show several similarities and differences between SalmoSim and real salmon. Again, however, the overall pattern is that of limited change in both in vivo and in vitro systems in response to the change in diet. Invariance observed in the microbiome in response to feed was reflected in estimates of physical attributes of fish in response to the change in feed formulation. As such, no statistically significant differences in various phenotypic measurements (fish length, weight, gutted weight, carcass yield, gonad, and liver weights) were noted in salmon fed on the two different diets used in the experiment (see Supplementary Figure 8). VFA measurements were undertaken to assay any differences between the microbial fermentation profiles of SalmoSim microbes fed on the different diets. Consistent with a limited impact of the two different feeds observed in both in vivo and in vitro datasets, invariance was also observed in in vitro VFA production data, in which no significant differences were observed in SalmoSim between the FMD and FMF diets (see Supplementary Figure 9). Finally, while ammonia production levels were largely static throughout the experiment, protein levels in the system did fluctuate (Figure S5), however, not apparently in response to the change in feed.
Microbial composition in SalmoSim and real salmon fed different feeds
Most gut compartments for real salmon, SalmoSim, and the salmon used to inoculate SalmoSim were abundant in Pseudomonas, Psychrobacter, and Staphylococcus genera, suggesting that genera present in the marine phase-salmon are generally maintained in SalmoSim (Fig. 5), with these three genera accounting for 42%, 39%, and 34% of all OTUs in non-inocula salmon, inoculum salmon, and SalmoSim samples, respectively. In terms of change in alpha diversity, the only statistically significant difference in response to the switch in feed was observed in the pyloric caeca compartment within SalmoSim based on the Shannon diversity metric (Supplementary Figure 6), where a slight decrease alongside the FMF occurred. Otherwise, the change in feed formulation did not impact alpha diversity in any gut compartment, either in real salmon or in SalmoSim. Furthermore, no differences were indicated between real salmon and SalmoSim within each gut compartment (Supplementary Figure 7).
To provide an overview of microbial composition and variation in the experiment, a PCoA (Principal Coordinates Analysis) based on Bray–Curtis distance was performed and plotted (Fig. 6A–D). The biological replicate (the fish providing the founding inoculum of each SalmoSim run) appears to be a major driver of community composition in the experiment (Fig. 6A). Taxonomic composition represented in stacked bar plots in Fig. 5 also supports this observation. Once individual SalmoSim runs (biological replicates) are visualised separately, changes to microbial communities in response to the feed become apparent (Fig. 6B–D). Statistical comparisons based on PERMANOVA show there is an effect of feed on microbial composition in both salmon and SalmoSim (Table 3); however, based on OTU differential abundance data (above), the effect seems to be small (Fig. 4). Samples from the real salmon fed on the different diets also diverge from one and other (supported by Table 3, Fig. 5), however, not necessarily along the same axes as each SalmoSim replicate. This divergence is potentially indicative of an effect of the biological replicate (i.e., inter-individual variation). Consistent with Fig. 5, inoculum for the respective SalmoSim replicates cluster among SalmoSim samples for the FMD in each case.
Our findings suggest a loss of microbial taxa diversity and richness as a result of transferring initial inoculums from real salmon into the SalmoSim system in the pyloric caeca and midgut compartments. Several lines of evidence from our core OTU analysis suggest that low prevalence (rare) OTUs make up most of the taxa lost, and progressive removal of rare OTUs increased the compositional similarity between inoculum and SalmoSim samples using both phylogenetic and ecological distances. Furthermore, shared OTUs between inoculum salmon and SalmoSim samples accounted for around 97% and 98% of total reads respectively, and as such, the microbiota in the model are highly representative of those from the fish that founded them. A general trend was observed, in which all gut compartments became increasingly stable throughout the 40-day experiment, with little-observed impact of introducing the different feed at day 20. Comparison of non-inoculum salmon and SalmoSim samples at the microbial level showed significant differences using both ecological and phylogenetic metrics as well as a different number of OTUs (SalmoSim 508 OTUs and real salmon samples 192 OTUs). These differences may be explained by the fact that samples used for non-inoculum salmon and SalmoSim originated from different individuals, whereas initial inoculum and SalmoSim samples for a given run originated from the same fish. Furthermore, the non-inoculum salmon and inoculum samples were derived slightly differently—with a longer section of the gut sampled for the inoculum samples, which could explain the higher diversity of bacteria (number of OTUs) which in turn can affect beta diversity output. However, despite this, shared OTUs between non-inoculum salmon and inoculum samples accounted for around 98% of total reads generated in non-inoculum real salmon and around 70% within inoculum real salmon samples. Correspondingly, we observed that the biological replicate (the founding inoculum of each SalmoSim run that originated from different individuals) was the major driver of community composition in the experiment. Once the individual runs were separated, phylogenetic and ecological distances suggested that changing feed was a statistically significant driver of community composition in both real salmon and SalmoSim. However, the vast majority of OTUs remained unchanged by the switch in feed in both systems, and no changes were noticeable in the bacterial activity (VFA production) within the system after the introduction of the plant-based feed, nor in phenotypic performance of Atlantic salmon fed on two different feeds (fish length, weight, gutted weight, carcass yield, gonad, and liver weights) in in vivo trial.
Many of the microbes we detected and cultured from the salmon gut microbiome have been reported previously in this species. For example, gram-negative Pseudomonas and Psychrobacter, the most abundant genera we observed, are among the core bacterial taxa known to reside within the real salmon gut [38, 60, 61]. Staphylococcus genera have also been reported widely in freshwater- and marine-farmed salmon . SalmoSim was able to maintain these species in culture throughout the experimental run, and although some diversity was lost, no statistical differences could be detected between the composition of SalmoSim and that of the fish gut communities used to found the different biological replicates via ecological metrics. Notable by their scarcity were mycoplasma OTUs, which occurred at relatively low abundance in both the in vivo and in vitro systems in this study. Mycoplasma OTUs were recovered from most SalmoSim gut compartments at low abundances (see supplementary Table 5), suggesting that these fastidious microbes can survive in the bioreactors. Our group and several others have widely reported Mycoplasma species from marine and freshwater stage of wild and farmed Atlantic salmon (higher abundance in farmed salmon), where many proliferate intracellularly in the gut epithelial lining [33, 63,64,65]. Establishing whether mycoplasma can actively proliferate in SalmoSim would require the use of founding communities rich in these organisms in a future experiment. One means of achieving this could be by using mock communities to better understand the temporal development of the observed microbial compositions from the inocula to the SalmoSim community .
We identified that a change in feed resulted in an overall shift in microbial community structure in both real salmon and SalmoSim system, as has been found to be the case in many previous studies [10, 67, 68], but the vast majority of OTUs within both real salmon and SalmoSim were not affected by the feed change. The direction of this shift, and the microbial taxa involved, were not equivalent in SalmoSim and real salmon, although no overall trend was observed at higher taxonomic levels in either system. Importantly, it is also the case that the vast majority of OTUs within both real salmon and SalmoSim were not affected by the switch in feed. Furthermore, it was found that change in feed did not affect VFA production in the SalmoSim system. As such, it is not clear whether any relevant functional shifts occurred in the microbiome of SalmoSim or real salmon as a result of the treatment. Furthermore, we did not identify any phenotypic changes (fish length, weight, gutted weight, carcass yield, gonad, and liver weights) within in vivo trial of Atlantic salmon fed on two different feeds. This lack of change is not unexpected, considering the plant-based feed was developed to have similar macronutrient composition to a Fishmeal-based feed. One difference is a slightly higher crude fibre (fermentable substrate) proportion in Fishmeal-free diet, which could explain higher microbial diversity in in vivo samples fed on plant-based feed .
The use of in vitro systems to study and model the microbial communities of monogastric vertebrates is becoming increasingly widespread, with systems simulating Sus scrofa (pig) , Gallus gallus (chicken) , Canis lupus (dog) , and other vertebrate guts. Using in vitro gut simulators is also a widely accepted approach to study the human gut microbiome [18, 72, 73]. One of the most established systems is the Simulator of the Human Intestinal Microbial Ecosystem (SHIME) that mimics the entire gastrointestinal tract incorporating the stomach, small intestine, and different colon regions . This system was used to study the effects of many different dietary additives on the human microbiome [74, 75]. The value of in vitro simulators in providing genuine insights is limited only by the research question and the corresponding level of sophistication required. The host component of the system, for example, is often poorly modelled, although cell lines, artificial mucosae, and digestion / absorbance systems can be included, which can provide specific physiological and metabolic insights [72, 76]. For instance, the growth medium in vessels representing stomach and small intestine of the SHIME system is enriched with pancreatic enzymes and bile, while a further upgrade of SHIME incorporates the mucosal environment allowing a portion of the microbiota present in the vessels to adhere to a gut mucus layer [18, 76]. Furthermore, more sophisticated in vitro gut models, such as TIM-2 (the TNO computer-controlled, dynamic in vitro gastroIntestinal Model of the colon), incorporates a dialysis system, which helps to maintain a physiological concentration of metabolites and prevents inhibition of the microbial growth by microbial metabolite accumulates [77, 78]. Finally, microfluid-based in vitro systems, such as HuMiX, allow coculturing of the host gut cells and microbial cells under conditions and processes mimicking gastrointestinal host-microbiome interface .
Whilst studies have been performed using a simple in vitro batch method to study fermentation of various feeds using gut inocula from Nile tilapia (Oreochromis niloticus) and European sea bass (Dicentrarchus labrax), these systems, to date, lack complexity . In this paper, we identified that in order to design and build an artificial gut simulator, parameters including physiochemical conditions (pH, temperature, dissolved oxygen), feed media composition, and flow rate between different gut compartments need to be tightly controlled. If all these parameters can be determined and maintained, and there is access to crude enzyme extracts from relevant species, adaptation of the SalmoSim model to other fish species (e.g., trout, carp) should be possible alongside suitable validation. As we found, interindividual variability may be an important consideration, and adequate biological replication is necessary to enable reliable interpretation of results, a consideration that can be overlooked by even the most sophisticated systems. Interindividual variation in gut microbial communities is widely observed in human studies that demonstrate more between-person variation than within-person variation with adults having an average unique microbial signature that is largely stable over time [81,82,83,84,85]. This is also true in Atlantic salmon—our previous work clearly shows high levels of interindividual variability in farmed  and wild  fish, as does the work of all others. It was reported, for example, that a single Lactobacillales OTU represented 96% of the microbiome of one fish which compared with a mean of only 3.5% relative abundance in the other fish from the same shoal in an aquaculture setting . Well-established gut microbiome systems, such as SHIME, use inoculums from only one donor  or in recent studies 2 donors in which differences were found in microbiota distribution even when using control diets . Furthermore, some artificial gut systems pool biological replicates together to produce a “representative microbiome inoculum”, such as in a recent in vitro chicken gut model, and even in these systems, microbiomes still show variability . To our knowledge in this study, we are the first to run a gut microbiome model in biological triplicate and to highlight the importance of accounting for interindividual differences before drawing conclusions. Prior to the current study, only one other attempt was made to study the effect of diet on Atlantic salmon gut microbial composition in vitro . In this preliminary study, a simple in vitro system was used to assess the impact of different feed formulations on the microbial communities of faecal slurries prepared from live salmon. However, no direct comparison was made with a true in vivo trial, nor were the different gut compartments present in salmon modelled in any detail and the predictive value for such simple in vitro systems in not immediately clear. Nonetheless, the work provided an important catalyst for the development of more sophisticated systems.
In the future, SalmoSim could be improved even further by introducing the ability to assess the digestibility of various feeds. In the aquaculture industry, determination of the digestibility of nutrients in various feeds provides the indication of their nutritional value and is often considered as the first step in feed quality evaluation [91,92,93,94]. While some in vitro systems used to assess digestibility in Atlantic salmon use, enzymes sourced from pigs, and bovines, such as porcine pepsin/porcine trypsin, bovine chymotrypsin, and porcine peptidase , enzyme extracts from different parts of the digestive system from salmon should be chosen to ensure accurate simulation of the in vivo digestion process. A modified SalmoSim could include the addition of enzyme extracts from the stomach and pyloric caecum compartments of real salmon in order to catabolise ingredients within the feed being tested , as well as suitable system for continuously removing the small molecular products of digestion [97,98,99].
Our results indicate that SalmoSim can not only stably maintain the most abundant microbial communities from real salmon, but also demonstrate similar responses to experimental feed treatments as those seen in real salmon. These results are encouraging, however, the nature of the treatment applied in this study—a switch between two similar feeds that had little effect on the gut microbiota in vivo—suggests that further experimentation with SalmoSim would be beneficial. For example, the survival and influence of probiotics within the system or the influence of known prebiotics could also be assessed, as has been previously studied in other in vitro gut systems . Gut models such as SalmoSim could provide a powerful tool for aquaculture, where there is considerable interest associated with the development of feed and feed additives [100,101,102], but where the capacity for in vivo trials is limited. The steps of the system development described in this paper can be used as guidelines to develop various other systems representing other fish species. The aim of such systems could be to provide prescreening tool for new feed ingredients and additives with the aim of reducing the cost and scale of in vivo testing. In parallel, an in vitro gut model for salmon could also be exploited to understand questions of public health importance (e.g., antimicrobial resistance and transfer ), as well as the fundamental ecological processes that underpin microbiome dynamics and assembly.
FAO. The State of Fisheries and Aquaculture in the world 2018: FAO; 2018.
Moffitt CM, Cajas-Cano L. Blue growth: the 2014 FAO State of World Fisheries and Aquaculture: Fisheries; 2014.
Food and Agriculture Organization. The State of World Fisheries and Aquaculture (SOFIA) | FAO | Food and Agriculture Organization of the United Nations: United Nations; 2016.
Junning C, Zhou X, Xue Y, Lucentea D, Laganaa C. Top 10 species groups in global aquaculture 2017. Food Agric Organ United Nations. 2019;June:7 Available from: www.fao.org/fishery/statistics/software/wapi/en.
Merrifield DL, Dimitroglou A, Foey A, Davies SJ, Baker RTM, Bøgwald J, et al. The current status and future focus of probiotic and prebiotic applications for salmonids. Aquaculture. 2010.
Ytrestøyl T, Aas TS, Åsgård T. Utilisation of feed resources in production of Atlantic salmon (Salmo salar) in Norway. Aquaculture. 2015;448:365–74.
Egerton S, Culloty S, Whooley J, Stanton C, Ross RP. The gut microbiota of marine fish. Front Microbiol Frontiers. 2018;9:873.
Perry WB, Lindsay E, Payne CJ, Brodie C, Kazlauskaite R. The role of the gut microbiome in sustainable teleost aquaculture. Proc R Soc B Biol Sci. 2020;287.
Ingerslev HC, Strube ML, von Gersdorff Jørgensen L, Dalsgaard I, Boye M, Madsen L. Diet type dictates the gut microbiota and the immune response against Yersinia ruckeri in rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol. 2014;40(2):624–33.
Gajardo K, Jaramillo-Torres A, Kortner TM, Merrifield DL, Tinsley J, Bakke AM, et al. Alternative protein sources in the diet modulate microbiota and functionality in the distal intestine of Atlantic salmon (Salmo salar). Appl Environ Microbiol. 2017;83(5):e02615–6.
Beheshti Foroutani M, Parrish CC, Wells J, Taylor RG, Rise ML, Shahidi F. Minimizing marine ingredients in diets of farmed Atlantic salmon (Salmo salar): effects on growth performance and muscle lipid and fatty acid composition. Soengas JL, editor. PLoS One. 2018 [cited 2020 Feb 3];13:e0198538. Available from: http://dx.plos.org/10.1371/journal.pone.0198538
Green TJ, Smullen R, Barnes AC. Dietary soybean protein concentrate-induced intestinal disorder in marine farmed Atlantic salmon, Salmo salar is associated with alterations in gut microbiota. Vet Microbiol. 2013;166:286–92.
Michl SC, Beyer M, Ratten JM, Hasler M, LaRoche J, Schulz C. A diet-change modulates the previously established bacterial gut community in juvenile brown trout (Salmo trutta). Sci Rep. 2019;9(1):2339.
Gupta S, Fečkaninová A, Lokesh J, Koščová J, Sørensen M, Fernandes J, et al. Lactobacillus dominate in the intestine of atlantic salmon fed dietary probiotics. Front Microbiol. 2019;9:3247.
Payne AN, Zihler A, Chassard C, Lacroix C. Advances and perspectives in in vitro human gut fermentation modeling. Trends Biotechnol. 2012;30:17–25.
Drieschner C, Könemann S, Renaud P, Schirmer K. Fish-gut-on-chip: development of a microfluidic bioreactor to study the role of the fish intestine: in vitro. Lab Chip. 2019;19(19):3268–76.
Possemiers S, Verthé K, Uyttendaele S, Verstraete W. PCR-DGGE-based quantification of stability of the microbial community in a simulator of the human intestinal microbial ecosystem. FEMS Microbiol Ecol. 2004;49:495–507.
Van Den Abbeele P, Grootaert C, Marzorati M, Possemiers S, Verstraete W, Gérard P, et al. Microbial community development in a dynamic gut model is reproducible, colon region specific, and selective for bacteroidetes and Clostridium cluster IX. Appl Environ Microbiol. 2010;76:5237–46.
Molly K, Woestyne MV, Smet ID, Verstraete W. Validation of the simulator of the human intestinal microbial ecosystem (SHIME) reactor using microorganism-associated activities. Microb Ecol Health Dis. 1994;7(4):191–200.
Postollec F, Falentin H, Pavan S, Combrisson J, Sohier D. Recent advances in quantitative PCR (qPCR) applications in food microbiology. Food Microbiol. 2011;28:848–61.
Malla MA, Dubey A, Kumar A, Yadav S, Hashem A, Allah EFA. Exploring the human microbiome: the potential future role of next-generation sequencing in disease diagnosis and treatment. Front Immunol. 2019;9:2868.
Ravi V, Venkatesh B. Rapidly evolving fish genomes and teleost diversity. Curr Opin Genet Dev Elsevier Current Trends. 2008;18:544–50.
Horn MH, Gawlicka A. Digestive system of fish. Encycl Life Sci. 2001.
Smith LS. Digestion in teleost fishes: Lect. Present. FAO/UNDP Train. course fish Feed. Technol. ACDP/REP/80/11; 1980. p. 3–17.
Mello GCG, Santos ML, Arantes FP, Pessali TC, Brito MFG, Santos JE. Morphological characterisation of the digestive tract of the catfish Lophiosilurus alexandri Steindachner, 1876 (Siluriformes, Pseudopimelodidae). Acta Zool. 2019;100.
Krogdahl Å, Sundby A, Holm H. Characteristics of digestive processes in Atlantic salmon (Salmo salar). Enzyme pH optima, chyme pH, and enzyme activities. Aquaculture. 2015.
Smith LS. Chapter 1. Digestion in teleost fishes. Fish Feed Technol. 1980.
Sahlmann C, Gu J, Kortner TM, Lein I, Krogdahl Å, Bakke AM. Ontogeny of the digestive system of Atlantic salmon (Salmo salar L.) and effects of soybean meal from start-feeding. PLoS One. 2015;10.
Buddington RK, Diamond JM. Aristotle revisited: the function of pyloric caeca in fish. Proc Natl Acad Sci U S A. 1986;83.
Buddington RK, Diamond JM. Pyloric ceca of fish: a “new” absorptive organ. Am J Physiol Gastrointest Liver Physiol. 1987;252.
Bakke-McKellep AM, Nordrum S, Krogdahl Å, Buddington RK. Absorption of glucose, amino acids, and dipeptides by the intestines of Atlantic salmon (Salmo salar L.). Fish Physiol Biochem. 2000;22.
Denstadli V, Vegusdal A, Krogdahl Å, Bakke-Mckellep AM, Berge GM, Holm H, et al. Lipid absorption in different segments of the gastrointestinal tract of Atlantic salmon (Salmo salar L.). Aquaculture. 2004;240.
Heys C, Cheaib B, Busetti A, Kazlauskaite R, Maier L, Sloan WT, et al. Neutral processes dominate microbial community assembly in Atlantic salmon, <em>Salmo salar</em>. Appl Environ Microbiol. 2020:AEM.02283-19 Available from: http://aem.asm.org/content/early/2020/02/03/AEM.02283-19.abstract.
Waagbø R, Jørgensen SM, Timmerhaus G, Breck O, Olsvik PA. Short-term starvation at low temperature prior to harvest does not impact the health and acute stress response of adult Atlantic salmon. PeerJ. 2017; [cited 2018 Mar 5];5:e3273. Available from: https://peerj.com/articles/3273.
Aas TS, Sixten HJ, Hillestad M, Sveier H, Ytrestøyl T, Hatlen B, et al. Measurement of gastrointestinal passage rate in Atlantic salmon (Salmo salar) fed dry or soaked feed. Aquac Rep. 2017.
Claassen S, du Toit E, Kaba M, Moodley C, Zar HJ, Nicol MP. A comparison of the efficiency of five different commercial DNA extraction kits for extraction of DNA from faecal samples. J Microbiol Methods. 2013.
Werner JJ, Koren O, Hugenholtz P, Desantis TZ, Walters WA, Caporaso JG, et al. Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys. ISME J. 2012;6(1):94–103.
Gajardo K, Rodiles A, Kortner TM, Krogdahl Å, Bakke AM, Merrifield DL, et al. A high-resolution map of the gut microbiota in Atlantic salmon (Salmo salar): a basis for comparative gut microbial research. Sci Rep. 2016; [cited 2020 Aug 27];6:1–10. Available from: www.nature.com/scientificreports/.
Joshi N, Fass J. Sickle - a windowed adaptive trimming tool for FASTQ files using quality. (Version 133); 2011.
Nikolenko SI, Korobeynikov AI, Alekseyev MA. BayesHammer: Bayesian clustering for error correction in single-cell sequencing. BMC Genomics. 2013;14 Suppl 1(Suppl 1):S7.
Masella AP, Bartram AK, Truszkowski JM, Brown DG, Neufeld JD. PANDAseq: paired-end assembler for illumina sequences. BMC Bioinformatics. 2012;13:31.
Schirmer M, D’Amore R, Ijaz UZ, Hall N, Quince C. Illumina error profiles: resolving fine-scale variation in metagenomic sequencing data. BMC Bioinformatics. 2016;17:125.
Mukherjee S, Stamatis D, Bertsch J, Ovchinnikova G, Katta HY, Mojica A, et al. Genomes OnLine database (GOLD) v.7: updates and new features. Nucleic Acids Res. 2019;47(D1):D649–59.
Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584.
Schmieder R, Edwards R. Fast identification and removal of sequence contamination from genomic and metagenomic datasets. PLoS One. 2011;6(3):e17288.
Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: machine learning in Python. J Mach Learn Res. 2011;12:2825–30.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.
Douglas GM, Maffei VJ, Zaneveld J, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2: an improved and extensible approach for metagenome inference. bioRxiv. 2019.
Rstudio Team. RStudio: integrated development for R. Boston MA: RStudio, Inc.; 2019.
Lagkouvardos I, Fischer S, Kumar N, Clavel T. Rhea: a transparent and modular R pipeline for microbial profiling based on 16S rRNA gene amplicons. PeerJ. 2017;5:e2836.
Ssekagiri A, Sloan WT, Zeeshan Ijaz U. microbiomeSeq: an R package for analysis of microbial communities in an environmental context. ISCB Africa ASBCB Conf. 2017.
McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217.
Jost L. Entropy and diversity. Oikos. 2006.
Jost L. Partitioning diversity into independent alpha and beta components. Ecology. 2007;88(10):2427–39.
Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB. Package vegan. R Pack Ver. 2013.
Anderson MJ. A new method for non-parametric multivariate analysis of variance. Aust Ecol. 2001;26(1):32–46.
Love M, Anders S, Huber W. Analyzing RNA-seq data with DESeq2. Bioconductor. 2017.
Navarrete P, Espejo RT, Romero J. Molecular analysis of microbiota along the digestive tract of juvenile Atlantic salmon (Salmo salar L.). Microb Ecol. 2009;57:550–61.
Webster TMU, Consuegra S, Hitchings M, de Leaniz CG. Interpopulation variation in the Atlantic salmon microbiome reflects environmental and genetic diversity. Appl Environ Microbiol. 2018;84(16):e00691–18.
Dehler CE, Secombes CJ, Martin SAM. Environmental and physiological factors shape the gut microbiota of Atlantic salmon parr (Salmo salar L.). Aquaculture. 2017;467:149–57.
Llewellyn MS, Boutin S, Hoseinifar SH, Derome N. Teleost microbiomes: the state of the art in their characterization, manipulation and importance in aquaculture and fisheries. Front Microbiol. 2014;5:207.
Cheaib B, Seghouani H, Ijaz UZ, Derome N. Community recovery dynamics in yellow perch microbiome after gradual and constant metallic perturbations. Microbiome. 2020;8(1):14.
Holben WE, Williams P, Saarinen M, Särkilahti LK, JHA A. Phylogenetic analysis of intestinal microflora indicates a novel Mycoplasma phylotype in farmed and wild salmon. Microb Ecol. 2002;44(2):175–85.
Amos GCA, Logan A, Anwar S, Fritzsche M, Mate R, Bleazard T, et al. Developing standards for the microbiome field. Microbiome. 2020;8(1):98.
Egerton S, Wan A, Murphy K, Collins F, Ahern G, Sugrue I, et al. Replacing fishmeal with plant protein in Atlantic salmon (Salmo salar) diets by supplementation with fish protein hydrolysate. Sci Rep. 2020;10(1):4194.
Michl SC, Ratten J-M, Beyer M, Hasler M, La Roche J, Schulz C. The malleable gut microbiome of juvenile rainbow trout (Oncorhynchus mykiss): diet-dependent shifts of bacterial community structures. Prunet P, editor. PLoS One. 2017 [cited 2020 Sep 4];12:e0177735. Available from: https://dx.plos.org/10.1371/journal.pone.0177735
Tanner SA, Berner AZ, Rigozzi E, Grattepanche F, Chassard C, Lacroix C. In vitro continuous fermentation model (PolyFermS) of the swine proximal colon for simultaneous testing on the same gut microbiota. PLoS One. 2014;9(4):e94123.
Card RM, Cawthraw SA, Nunez-Garcia J, Ellis RJ, Kay G, Pallen MJ, et al. An in vitro chicken gut model demonstrates transfer of a multidrug resistance plasmid from Salmonella to commensal Escherichia coli. MBio. 2017;8(4):e00777–17.
Duysburgh C, Ossieur WP, De Paepe K, Van Den Abbeele P, Vichez-Vargas R, Vital M, et al. Development and validation of the Simulator of the Canine Intestinal Microbial Ecosystem (SCIME). J Anim Sci. 2020;98(1):skz357.
Déat E, Blanquet-Diot S, Jarrige JF, Denis S, Beyssac E, Alric M. Combining the dynamic TNO-gastrointestinal tract system with a Caco-2 cell culture model: application to the assessment of lycopene and α-tocopherol bioavailability from a whole food. J Agric Food Chem. 2009;57:11314–20.
Kim HJ, Li H, Collins JJ, Ingber DE. Contributions of microbiome and mechanical deformation to intestinal bacterial overgrowth and inflammation in a human gut-on-a-chip. Proc Natl Acad Sci. 2016;113:E7–15 Available from: http://www.pnas.org/lookup/doi/10.1073/pnas.1522193112.
Sánchez-Patán F, Barroso E, Van De Wiele T, Jiménez-Girón A, Martín-Alvarez PJ, Moreno-Arribas MV, et al. Comparative in vitro fermentations of cranberry and grape seed polyphenols with colonic microbiota. Food Chem. 2015;183:273–82.
Giuliani C, Marzorati M, Innocenti M, Vilchez-Vargas R, Vital M, Pieper DH, et al. Dietary supplement based on stilbenes: a focus on gut microbial metabolism by the in vitro simulator M-SHIME®. Food Funct. 2016;7:4564–75 Available from: http://xlink.rsc.org/?DOI=C6FO00784H.
Van den Abbeele P, Roos S, Eeckhaut V, Mackenzie DA, Derde M, Verstraete W, et al. Incorporating a mucosal environment in a dynamic gut model results in a more representative colonization by lactobacilli. Microb Biotechnol. 2012;5:106–15.
Lefebvre DE, Venema K, Gombau L, Valerio LG, Raju J, Bondy GS, et al. Utility of models of the gastrointestinal tract for assessment of the digestion and absorption of engineered nanomaterials released from food matrices. Nanotoxicology. 2015;9:523–42 Available from: https://www.tandfonline.com/doi/full/10.3109/17435390.2014.948091.
Minekus M, Smeets-Peeters M, Bernalier A, Marol-Bonnin S, Havenaar R, Marteau P, et al. A computer-controlled system to simulate conditions of the large intestine with peristaltic mixing, water absorption and absorption of fermentation products. Appl Microbiol Biotechnol. 1999;53:108–14.
Shah P, Fritz JV, Glaab E, Desai MS, Greenhalgh K, Frachet A, et al. A microfluidics-based in vitro model of the gastrointestinal human-microbe interface. Nat Commun. 2016;7:11535. Available from: https://doi.org/10.1038/ncomms11535.
Leenhouwers JI, Pellikaan WF, Huizing HFA, Coolen ROM, Verreth JAJ, Schrama JW. Fermentability of carbohydrates in an in vitro batch culture method using inocula from Nile tilapia (Oreochromis niloticus) and European sea bass (Dicentrarchus labrax). Aquac Nutr. 2008;14.
Stearns JC, Lynch MDJ, Senadheera DB, Tenenbaum HC, Goldberg MB, Cvitkovitch DG, et al. Bacterial biogeography of the human digestive tract. Sci Rep. 2011;1.
Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R. Bacterial community variation in human body habitats across space and time. Science (80- ). 2009;326.
Rajilić-Stojanović M, Heilig HGHJ, Tims S, Zoetendal EG, De Vos WM. Long-term monitoring of the human intestinal microbiota composition. Environ Microbiol. 2013;15.
Human Microbiome Project Consortium T. A framework for human microbiome research The Human Microbiome Project Consortium*. Nature. 2012.
Huttenhower C, Gevers D, Knight R, Abubucker S, Badger JH, Chinwalla AT, et al. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486.
Heys C, Cheaib B, Busetti A, Kazlauskaite R, Maier L, Sloan WT, et al. Neutral processes dominate microbial community assembly in Atlantic Salmon, Salmo salar. Appl Environ Microbiol. 2020;86.
Llewellyn MS, McGinnity P, Dionne M, Letourneau J, Thonier F, Carvalho GR, et al. The biogeography of the Atlantic salmon (Salmo salar) gut microbiome. ISME J. 2016;10:1280–4.
Schmidt V, Amaral-Zettler L, Davidson J, Summerfelt S, Good C. Influence of fishmeal-free diets on microbial communities in Atlantic salmon (Salmo Salar) recirculation aquaculture systems. Appl Environ Microbiol. 2016;82.
Rovalino-Córdova AM, Fogliano V, Capuano E. Effect of bean structure on microbiota utilization of plant nutrients: an in-vitro study using the simulator of the human intestinal microbial ecosystem (SHIME®). J Funct Foods. 2020;73.
Zarkasi KZ, Taylor RS, Glencross BD, Abell GCJ, Tamplin ML, Bowman JP. In vitro characteristics of an Atlantic salmon (Salmo salar L.) hind gut microbial community in relation to different dietary treatments. Res Microbiol. 2017;168(8):751–9.
Allan GL, Parkinson S, Booth MA, Stone DAJ, Rowland SJ, Frances J, et al. Replacement of fish meal in diets for Australian silver perch, Bidyanus bidyanus: I. Digestibility of alternative ingredients. Aquaculture. 2000;186.
Glencross BD, Booth M, Allan GL. A feed is only as good as its ingredients - a review of ingredient evaluation strategies for aquaculture feeds. Aquac Nutr. 2007;13(1):17–34.
Luo Z, Li XD, Gong SY, Xi WQ. Apparent digestibility coefficients of four feed ingredients for Synechogobius hasta. Aquac Res. 2009;40.
Liu H, Wu X, Zhao W, Xue M, Guo L, Zheng Y, et al. Nutrients apparent digestibility coefficients of selected protein sources for juvenile Siberian sturgeon (Acipenser baerii Brandt), compared by two chromic oxide analyses methods. Aquac Nutr. 2009;15.
Stewart Anderson J, Lall SP, Anderson DM, McNiven MA. Evaluation of protein quality in fish meals by chemical and biological assays. Aquaculture. 1993;115.
Moyano FJ, Saénz de Rodrigáñez MA, Díaz M, Tacon AGJ. Application of in vitro digestibility methods in aquaculture: constraints and perspectives. Rev Aquac. 2015;7.
Yúfera M, Romero MJ, Pujante IM, Astola A, Mancera JM, Sánchez-Vázquez FJ, et al. Effect of feeding frequency on the daily rhythms of acidic digestion in a teleost fish (gilthead seabream). Chronobiol Int. 2014;31.
Montoya A, López-Olmeda JF, Yúfera M, Sánchez-Muros MJ, Sánchez-Vázquez FJ. Feeding time synchronises daily rhythms of behaviour and digestive physiology in gilthead seabream (Sparus aurata). Aquaculture. 2010;306.
Martínez-Llorens S, Peruzzi S, Falk-Petersen IB, Godoy-Olmos S, Ulleberg LO, Tomás-Vidal A, et al. Digestive tract morphology and enzyme activities of juvenile diploid and triploid Atlantic salmon (Salmo salar) fed fishmeal-based diets with or without fish protein hydrolysates. PLoS One. 2021;16.
Encarnação P. Functional feed additives in aquaculture feeds. Aquafeed Formul. 2016.
Kristiansen M, Merrifield DL, Vecino JLG, Myklebust R, Ringø E. Evaluation of prebiotic and probiotic effects on the intestinal gut microbiota and histology of Atlantic salmon (Salmo salar L.). J Aquac Res Dev. 2011.
Hartviksen M, Vecino JLG, Ringø E, Bakke AM, Wadsworth S, Krogdahl Å, et al. Alternative dietary protein sources for Atlantic salmon (Salmo salar L.) effect on intestinal microbiota, intestinal and liver histology and growth. Aquac Nutr. 2014.
We thank the Llewellyn Environmental Biotechnology Laboratories teams for their help in sampling. Big thanks to the MOWI team, especially Jarred Lee Knapp and Matthew Watkins-Baker, in Averøy Norway, for sampling and providing physiological data for real salmon.
Availability of data and material
Sequence data have been deposited alongside metadata to the NCBI Short Read Archive.
This research was supported in part by research grants from the BBSRC (grant number BB/P001203/1 & BB/N024028/1), by Science Foundation Ireland, the Marine Institute, and the Department for the Economy, Northern Ireland, under the Investigators Program grant number SFI/15/IA/3028, and by the Scottish Aquaculture Innovation Centre. U.Z.I. is supported by a NERC independent research fellowship (NERC NE/L011956/1) as well as a Lord Kelvin Adam Smith Leadership Fellowship (Glasgow). RK is supported by an Alltech PhD Studentship award to the University of Glasgow.
Ethics approval and consent to participate
Animals sampled in the study were euthanised by authorised MOWI employees under Home Officer Schedule 1 of the Animals (Scientific Procedures) Act 1986.
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.
16S rRNA gene-targeted group-specific primers. Table lists primer sets that already published and validated in the literature. All primers were used on mouse faeces samples apart from Alphaproteobacterial specific primers that were used on marine biofilm samples.
First round PCR primers used for the first round of NGS library preparation.
Table S3 Second round PCR primers used for the first round of NGS library preparation.
OTUs prevalence analysis by sub-setting full dataset into multiple core OTUs. The table summarises the number OTUs within each subset dataset (subset by the % of samples that share OTUs) and percentage of the total number of OTUs within the full dataset (100%). This table also shows the number of reads and the percentage of total reads (100%) within each of the subset datasets. Note: in the 60% subset three samples were lost as they did not retain any under that criteria OTUs: Id-val1-PC1, Id-Val2-MG4 and Id-Val2-PC1.
Bacterial group responses to feed change within different gut compartments in real salmon and SalmoSim based on qPCR data. The table summarises the Estimated Marginal Means output for each mixed-effect linear model run with different qPCR measured relative abundance values identifying the difference between real salmon and SalmoSim response to feed change (Fish meal to Fish meal free diet) within different gut compartments (S – stomach, PC – pyloric caeca, and MG – mid gut). P>0.05 values identify no change in the bacterial group, p<0.05 identifies decrease (Est is negative), and p<0.05 identifies increase (Est is positive) in the relative abundance of target group after the feed change. Bold values identify similarities between SalmoSim and real salmon samples. The SalmoSim values used for this test involves stable SalmoSim time points: days 16, 18 and 20 (Fish meal diet – once bacterial communities adapted to the SalmoSim environment), and days 36, 38 and 40 (Fish meal free diet – once bacterial communities adapted to feed change).
Physiochemical conditions measured within different real Atlantic salmon gut compartments. 1A-1C measured physicochemical conditions within real salmon (n=3) gut compartments: pH (1A), temperature (°C, 1B), dissolved oxygen (mg/L, 1C).
Specificity of the primers that target Lactobacillus and Mycoplasma genus. The results in figure summarise bacterial genus targeted by Lactobacillus (Fig. 1 A, B, C) and Mycoplasma (Fig. 1 D, E, F) specific primer set. It shows that of all genus captured by Lactobacillus primer pair 98% were Lactobacillus in fish 1, 78% in fish 2, and 65% in fish 3. While of all genus captured by Mycoplasma primer pair 95% were Mycoplasma in fish 1, 79% in fish 2, and 56% in fish 3.
Stability within SalmoSim system calculated by using unweighted and weighted UniFrac values for pairwise beta diversity analysis. The figure represents microbial stability within the SalmoSim system (data from all gut compartments combined) as the pairwise beta diversity comparison between different sampling time points (days), calculated by using A unweighted (0%) and B weighted (100%) UniFrac as a distance measure. A small p-value indicates that the two time points are statistically different, and p>0.05 indicates that two time points are not statistically different. The colour key illustrates the p-value: red end of spectrum denoting low p values (distinct compositions between time points) and dark green indicating high p values (similar compositions between timepoints).
Measured value (qPCR, ammonia and protein concentrations) stability within different SalmoSim compartments fed on Fish meal and Fish meal free diets. The figure summarises the Estimated Marginal Means output for each mixed-effect linear model (Model 1) run with different values measured in different SalmoSim compartments (qPCR measurements, ammonia and protein concentrations) identifying the difference between different time points during the first (system fed on Fish meal diet) and last 20 days (system fed on Fish meal free diet) of validation experiment. A small p-value indicates that the two time points are statistically different, and p>0.05 indicates that two time points are not statistically different. The colour key illustrates the p-value: red end of spectrum denoting low p values (low correlation between time points) and dark green indicating high p values (no differences between timepoints).
Stability within SalmoSim system, within different biological replicates and different gut compartments, calculated by using generalised UniFrac values for pairwise beta diversity analysis. The figure represents microbial stability within the SalmoSim system (data separated by different biological replicates and gut compartments) as the pairwise beta diversity comparison between different sampling time points (days), calculated by using generalised (50%) UniFrac as a distance measure. A small p-value indicates that the two time points are statistically different, and p>0.05 indicates that two time points are not statistically different. The colour key illustrates the p-value: red end of spectrum denoting low p values (distinct compositions between time points) and dark green indicating high p values (similar compositions between timepoints).
Stability within SalmoSim system calculated by using different UniFrac values for pairwise beta diversity analysis. The figure represents microbial stability within the SalmoSim system (data from all gut compartments and two different technical replicate runs combined) as the pairwise beta diversity comparison between different sampling time points (days), calculated by using unweighted (0%), generalised (50%) and weighted (100%) UniFrac as a distance measure. A small p-value indicates that the two time points are statistically different, and p>0.05 indicates that two time points are not statistically different. The colour key illustrates the p-value: red end of spectrum denoting low p values (distinct compositions between time points) and dark green indicating high p values (similar compositions between timepoints).
Calculated alpha-diversity metrics within different gut compartments of real salmon and SalmoSim fed on Fish meal and Fish meal free diets. Figure visually represents different alpha diversity outputs within different gut compartments of real salmon in red and SalmoSim in yellow (stable time points: 16, 18 and 20 fed on Fish meal, and 36, 38 and 40 fed on Fish meal free diet) fed on Fish meal and Fish meal free diets. A visually represents effective richness (number of OTUs), B represents effective Shannon diversity. The lines above bar plots represent statistically significant differences after feed change. The stars flag the levels of significance: one star (*) for p-values between 0.05 and 0.01, two stars (**) for p-values between 0.01 and 0.001, and three stars (***) for p-values below 0.001.
In vivo phenotypic fish performance fed on two different feeds. Figure visually represents different phenotypic performance data of fish (n=32 per feed) fed on two different feed. A Atlantic salmon length in centimetres; B Atlantic salmon length in weight in kilograms; C Atlantic salmon percentage carcass yield; D Atlantic salmon gonad weight in grams; E Atlantic salmon gutted weight in kilograms; F Atlantic salmon liver weight in grams. Blue box plots represent data from salmon (n=32) fed on Fish meal free diet, and red represents Atlantic salmon fed on Fish meal diet (n=32).
VFA production within different SalmoSim compartments fed on different feeds. Figure represents 11 volatile fatty acid production within SalmoSim system fed on Fish meal and Fish meal free diets within different gut compartments. Y axis represents the concentration of specific volatile fatty acid (mM) while the X axis represents each gut compartment (stomach, pyloric caeca, midgut). Red colour denoted Fish meal and blue – Fish meal free diets. The lines above bar plots represent statistically significant differences between different feeds and gut compartments. The stars flag the levels of significance: one star (*) for p-values between 0.05 and 0.01, two stars (**) for p-values between 0.01 and 0.001, and three stars (***) for p-values below 0.001.
About this article
Cite this article
Kazlauskaite, R., Cheaib, B., Heys, C. et al. SalmoSim: the development of a three-compartment in vitro simulator of the Atlantic salmon GI tract and associated microbial communities. Microbiome 9, 179 (2021). https://doi.org/10.1186/s40168-021-01134-6
- Gut microbiota
- In vitro gut model system
- Atlantic salmon