- Research article
- Open Access
Gastrointestinal microbial community changes in Atlantic cod (Gadus morhua) exposed to crude oil
BMC Microbiologyvolume 18, Article number: 25 (2018)
The expansion of offshore oil exploration increases the risk of marine species being exposed to oil pollution in currently pristine areas. The adverse effects of oil exposure through toxic properties of polycyclic aromatic hydrocarbons (PAHs) have been well studied in Atlantic cod (Gadus morhua). Nevertheless, the fate of conjugated metabolites in the intestinal tract and their effect on the diversity of intestinal microbial community in fish is less understood. Here, we investigated the intestinal microbial community composition of Atlantic cod after 28 days of exposure to crude oil (concentration range 0.0–0.1 mg/L).
Analysis of PAH metabolites in bile samples confirmed that uptake and biotransformation of oil compounds occurred as a result of the exposure. Various evidence for altered microbial communities was found in fish exposed to high (0.1 mg/L) and medium (0.05 mg/L) concentrations of oil when compared to fish exposed to low oil concentration (0.01 mg/L) or no oil (control). First, altered banding patterns were observed on denaturing gradient gel electrophoresis for samples pooled from each treatment group. Secondly, based on 16S rRNA sequences, higher levels of oil exposure were associated with a loss of overall diversity of the gut microbial communities. Furthermore, 8 operational taxonomic units (OTUs) were found to have significantly different relative abundances in samples from fishes exposed to high and medium oil concentrations when compared to samples from the control group and low oil concentration. Among these, only one OTU, a Deferribacterales, had increased relative abundance in samples from fish exposed to high oil concentration.
The results presented herein contribute to a better understanding of the effects of oil contamination on the gut microbial community changes in fish and highlight the importance of further studies into the area. Our findings suggest that increased relative abundance of bacteria belonging to the order Deferribacterales may be indicative of exposure to oil at concentrations higher than 0.05 mg/L.
Marine environments continue to be exposed to oil pollution events as oil and gas production activities continue to expand . Accidental and operational releases of crude oil remain a threat to marine biota, including fish, therefore monitoring activities play an important role. Atlantic cod (Gadus morhua) is a commonly used fish species in environmental monitoring in the North Atlantic . For this species, the adverse effects of oil exposure, in particular, the damage caused by petrogenic polycyclic aromatic hydrocarbons (PAHs) − a class of compounds found in crude oil and oil products − is well documented [3, 4]. In addition, several biomarkers of petrogenic PAH exposure are well-established for this species [3, 5]. PAHs found in crude oil are well-known carcinogens due to their ability to form DNA and protein adducts following their primary metabolism in the liver [6, 7]. Liver enzymes initiate the excretion process of PAHs by transforming them into more water soluble conjugated metabolites, which are sometimes more toxic than their corresponding parent PAHs [8, 9]. The conjugated metabolites are then excreted via the bile into the intestinal tract. The fate of conjugated metabolites in the intestinal tract and their effect on the gut microbiota is still under investigation. It is generally assumed that conjugated biliary PAH metabolites either (1) pass through the gut unchanged and are eliminated via feces, or (2) are hydrolyzed spontaneously, reabsorbed through the intestinal wall and transported to the liver via the portal vein (enterohepatic recirculation) . Considering the metabolic diversity of the gut microbiota, conjugated metabolites could also undergo biotransformation as a result of bacterial activity in the gastrointestinal (GI) tract . Moreover, there might be mechanisms that exacerbate the toxic effects of pollutants by for example reactivating conjugated biliary metabolites . Studies on the potential interactions between biliary metabolites and gut microbiota in fish are scarce, and most microbiota studies have so far been carried out on mammals, e.g., mice and humans [13, 14]. Their results indicate that gut microbiota can indeed be actively involved in biotransformation of environmental pollutants. Such biotransformation may influence the detoxification process and alter the fate of xenobiotic compounds. At the same time, these pollutants can in return influence the composition of the gut microbiota.
The species composition, richness and diversity of the fish gut microbiome has been studied extensively in fish species that are interesting for the aquaculture industry mainly with the interest of improving fish health through enhanced feed. Commonly used methods include cultivation and/or molecular profiling with denaturing gradient gel electrophoresis (DGGE) . Based on the increasing use of next-generation sequencing techniques a diverse gut microbiome composition has been detected in fish [16, 17]. While next-generation sequencing  indicated changes in the intestinal community of a freshwater fish (goldfish, Carassius auratus) upon exposure to an environmental pollutant (pentachlorophenol), studies investigating such effect in marine fish are lacking. Herein, we report the effect of oil exposure on the intestinal microflora in Atlantic cod. The aim of this study was to assess the potential shift in gut bacterial community composition of Atlantic cod due to oil exposure. The intestinal bacterial community composition, as analyzed by Illumina sequencing of V4 amplicons of bacterial 16S rRNA genes, was compared among gut microbiota of fish after a 28-day exposure to four levels of oil concentration (i.e., no oil, 0.01, 0.05 and 0.1 mg/L). Results from biological effect measurement, i.e. PAH metabolite analysis in fish bile, were used to confirm metabolic uptake of oil by the fish .
Exposure set up and sampling
Atlantic cod (weight range of 137–1455 g) were caught in March–April 2015 in the Stavanger region, transported to the International Research Institute of Stavanger facilities and placed in 1000 L tanks. After a 1–2 weeks’ acclimation, fish were exposed to three different concentrations of dispersed crude oil from the Troll C platform (North Sea, Norway): low (0.01 mg/L) (n of fish = 8), medium (0.05 mg/L) (n of fish = 9) and high (0.1 mg/L) (n of fish = 6). A group of fish were kept in a tank with clean seawater flow and used as control (n of fish = 7). A constant flow of 8 L/min seawater into each tank was ensured in both control and exposed groups. The crude oil was continuously distributed from a header tank (dispersed oil reservoir) into the different exposure tanks through a continuous flow system (CFS) as described earlier , in order to mimic natural environmental conditions, e.g., oil spill scenario . Illustration of the CFS system is shown in an additional figure [see Additional file 1]. The CFS was checked daily, oil droplet sizes were monitored three times a week by Multisizer measurements, and PAH concentrations were measured semi quantitatively by fixed wavelength fluorescence (FF) weekly (data not reported). A detailed description of the exposure set up can be found in the thesis of Delgado .
After 28 days of exposure, fish were sedated using 5 mg/L Aquacalm fish sedative (Metomidate HCl) and sacrificed by a sharp blow to the head. The entire intestinal tract was removed, placed in a sterile container, snap-frozen in liquid nitrogen and stored at − 80 °C until DNA extraction. For analysis of PAH uptake and biotransformation in fish, bile samples were drawn using a syringe, snap-frozen in liquid nitrogen and stored at − 80 °C until analysis.
Chemical analysis was carried out in order to assess PAH content of the Troll C crude oil and of samples taken from the header tank on day 14 and 28 (Intertek West Lab AS, ISO28540:2011). In the Troll C crude oil, the EPA (Environmental Protection Agency) PAHs were quantified, while in the header tank, 25 PAH compounds and groups of compounds were measured.
Condition Index (CI) = (weight (g)/[length (cm)]3) × 100.
Hepatosomatic Index (HSI) = [liver weight (g)/fish weight (g)] × 100.
PAH metabolites in bile were determined by two methods, the semi-quantitative fixed wavelength fluorescence (FF) and gas chromatography mass spectrometry (GC-MS).
Fixed wavelength fluorescence (FF) analysis of bile samples was conducted as described by Aas et al. . Bile samples were diluted in 50% methanol and analysed using a Lumina fluorescence spectrometer (Thermo Scientific). The concentration of PAH metabolites was expressed as mg pyrene fluorescence equivalents (PFE)/mL bile.
GC-MS analysis was carried out as previously described [24, 25]. In brief, the OH-PAHs were extracted with ethylacetate, dried with anhydrous sodium sulphate and concentrated prior to sylilation with N,O-Bis(trimethylsilyl)trifluoroacetamide (BSFTA). Trimethylsilyl (TMS) ethers of OH-PAHs were analysed using an HP5890 series II gas chromatograph (Shimadzu QP2010) equipped with a CP-Sil 8 CB-MS (Varian) column. Mass spectra were obtained at 70 eV in selected ion mode (SIM).
Gastrointestinal microbial community
The GI tract of three fishes from each exposure group were used for microbial community composition analysis. After thawing on ice, each GI tract sample was cut into small pieces using sterile scissors and tweezer and placed in 15 mL sterile tubes with 3 mL ATL lysis buffer (Qiagen). Samples were vortexed vigorously several times and kept at room temperature (approx. 22 °C), then centrifuged briefly at 10,000 rpm. The liquid phase was decanted into a fresh 15 mL tube and contents were lysed overnight at 55 °C. Following this overnight incubation, a 30 min RNA digestion was carried out at 37 °C by adding 20 μL RNAse A solution (100 mg/mL, Qiagen) to each sample. DNA was then purified using phenol:chlorophorm:isoamylalcohol (PCI, 25:24:1, Sigma) extraction and ethanol precipitation. Briefly, to each lysate (approx. 3 mL), 5 mL of PCI was added and mixed vigorously. Samples were kept on ice followed by centrifugation at 5000 rpm (4 °C) for 15 min. The supernatant was pipetted into a new tube before the extraction was repeated one more time with PCI (5 mL) and then once with 5 mL chlorophorm:isoamylalcohol (CI, 24:1, Sigma). The DNA was finally precipitated with ice-cold ethanol (2 volumes of 100%) in the presence of sodium acetate (0.3 mM, pH 5.3). Samples were gently mixed and then incubated at − 20 °C for a minimum of 4 h. Falcon tubes were centrifuged at 5000 rpm for 30 min to collect the DNA pellet. Ethanol was discarded, a washing step with 70% ethanol was carried out and pellets were re-suspended in Tris-EDTA buffer (10 mM Tris-HCl, 1 mM EDTA, pH 7.8). Finally, an additional clean-up was performed using Genomic DNA Clean & Concentrator (Zymo Research) according to the manufacturer’s instructions, to completely remove inhibitors and resuspend DNA in molecular grade water.
Polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) was carried out on all individual genomic DNA samples and on pooled genomic DNA samples which were from the same exposure group. The V3-V5 region of the 16S rDNA was amplified as described previously  (forward primer 341f: 5‘-CCTACGGGAGGCAGCAG-3’ with a GC-clamp 5′- CGCCCGGGGCGCGCCCCGGGCGGGGCGGGGGCACGGGGGG-3′ attached at the 5′-end and reverse primer 907r: 5′- CCGTCAATTCMTTTGAGTTT-3′). Each 50 μL PCR mix was composed of molecular grade water, 5 μL PCR buffer (5 Prime), 0.1 mM of each dNTP’s, 10 pM of each primer, 1 μL of genomic DNA and 0.25 μL of Taq polymerase (5 U/μL, 5 Prime). The PCR program was as follows: initial activation at 94 °C for 2 min, followed by 25 cycles of: (1) denaturation at 94 °C for 30 s, (2) annealing at 55 °C for 40 s and (3) elongation at 72 °C for 1 min. A final elongation step at 72 °C for 7 min was also included. The DGGE was carried out on a 6% acrylamide gel containing denaturing agents, urea and formamide, in a gradient of 20–80% prepared by using the IngenyPhorU gradient former. Equal amounts of PCR products were loaded into each well (confirmed by agarose gel electrophoresis). The run was performed using an IngenyPhorU (Ingeny International BV, Goes, The Netherlands) system (17 L volume, 1 x TAE as buffer, temperature at 60 °C) at 90 V for 18 h. Bands were visualized using post-staining (GelRed in TAE) for 60 min prior to imaging.
Amplicon sequencing and bioinformatics analysis
Amplicon libraries were prepared using BioScientific’s NEXTflex V4 library preparation kit following the manufacturer’s instructions with minor modifications. All DNA samples were diluted 10 times prior to first PCR step. Primer and barcode sequences are shown as additional file [see Additional file 2]. Instead of the recommended 10 cycles, 25 cycles were used during the first PCR reaction in order to obtain sufficient amount of products. PCR clean-up was then performed using MinElute Column Clean-up kit (Qiagen) and cleaned reaction mixes were diluted. The second PCR step was performed using 25 cycles as well. Approximately 40 μL of PCR products were loaded on a 2% agarose gel, bands corresponding to expected product size were excised and the second cleanup step was performed with MinElute Gel Extraction kit (Quiagen). DNA concentration was then measured with a NanoVue instrument (GE Healthcare). Products from three repeated library preparation rounds were pooled together prior to amplicon sequencing. Finally, the amplicons were sequenced by the Norwegian Sequencing Center (NSC) using Illumina MiSeq (Illumina, San Diego, CA, USA) V3 chemistry reagents (600 cycle, paired end). PhiX was blended in at 30% and a maximum of 1 bp mismatch was allowed during demultiplexing. All sequences have been deposited in the European Nucleotide Archive (ENA), acc. nr. PRJEB21667.
Read qualities were analysed using FastQC , primer and adapter sequences were removed using “cutadapt”  (standard parameters with primer sequences), before remaining phiX sequences were removed by mapping the reads to the phiX reference genome (Genbank acc. nr. J02482.1) using “BWA mem” with default settings . Reads aligning were removed using seqtk with the subseq command . Downstream processing and sequence analysis was performed using Mothur (v. 1.36.1) according to the Mothur Illumina MiSeq SOP . Briefly, paired-end reads were joined using “make.contigs” and filtered based upon a minimum average Phred quality score of 25, contig length (max 289 bp) and a zero-tolerance to ambiguous base pairs. Chimeras were removed using the UCHIME algorithm  implemented in Mothur, before the sequences (V4 region) were aligned and classified using the Silva SEED bacterial database (v. 119) as reference. The dataset was clustered into OTUs based on a 97% similarity sequence similarity using average neighbour clustering. The OTU table output of Mothur was imported into R studio  (based on R (v. 3.2.3) ) for further processing using the R package “phyloseq” (v. 1.14.0) . Results were visualized using the R package “ggplot2” (2.0.0) .
A common scaling procedure  normalized the OTU count in a given library with a factor corresponding to the ratio of the smallest library size in the dataset to the library size of the OTU in question. This replaces rarefying (i.e. random sub-sampling to the lowest number of reads) as it effectively results in the library scaling one would achieve by averaging an infinite number of repeated sub-samplings. We then excluded all OTUs present in only one sample from the dataset. Next, rarefaction curves were calculated on the observed number of reads to confirm a sufficient sequencing effort using an R script from the Phyloseq repository website . Finally, alpha diversity measures and ordination distances for non-metric multidimensional scaling (NMDS) plots were calculated in Phyloseq. Distances were calculated based on OTUs with more than 10 reads (> 99% of total number of reads).
Differences in within-sample (fish specimen) diversity (alpha diversity) were tested using both analysis of variance (ANOVA) and the Wilcoxon rank-sum test. Further, differences in bacterial community composition (beta diversity) between the exposure groups were assessed using Permutational Multivariate Analysis of Variance (PERMANOVA) using the “adonis” function in the R package “Vegan” . The test was run with 20,000 permutations, applying both OTU-based (Bray-Curtis/Jaccard) and phylogeny-based (weighted and unweighted UniFrac) measures. PERMANOVA was also applied to test for differences in PAH metabolites between the exposure groups. PERMANOVA assumes the multivariate dispersion in the compared groups to be homogeneous; this was verified using the “betadisper” function in “Vegan”. Similarity percentage (SIMPER) procedure implemented in “Vegan” was used to quantify the contribution of individual OTUs to the overall Bray-Curtis dissimilarity between the groups. “MetaStats”  was used to test for OTUs and orders with statistically significant differential abundances between the exposure groups. Only OTUs containing more than 0.1% of the total number of reads in the normalized and filtered dataset were used in this analysis. Finally, a correspondence analysis (CA) (“cca” in Vegan) was performed to explore the gut microbial community structure. Environmental parameters (order abundances) that significantly (p < 0.01) correlated with the ordination were fitted onto the CA plot using the “envfit” command in Vegan.
Chemical analysis of the Troll C crude oil used in this exposure confirmed that all 16 EPA PAHs were present, their sum was 1.6 g/L. The most abundant PAHs present were 2- and 3-ring PAHs and the highest concentrations were found for those of 2-methylnaphthalene and 3-methylnaphthalene. The PAH profile of the seawater in the header tank also showed presence of both low and high molecular weight PAHs. Chemical analysis results for Troll C crude and header tank are summarized in additional files [see Additional files 3 and 4, respectively].
Morphometric parameters (CI and HSI) and concentrations of PAH metabolites in fish bile, reported for the three individuals used for GI microbial community analysis, were intended to give an overall description of the effect of oil on the exposed fishes. Moreover, concentrations of PAH metabolites in fish bile were used as support parameter for the GI microbial community analysis.
CI and HSI values, reported in additional tables [see Additional file 5], showed no differences between fish exposed to crude oil and the control group. The lack of significant difference between exposure groups was confirmed by statistical analysis using data from all exposed fishes (data not shown).
The results of both FF and GC-MS analysis of PAH metabolites in bile confirmed the uptake and biotransformation of PAHs after exposure and results can be seen in additional files [see Additional files 6 and 7, respectively]. Distribution pattern of PAHs found in Troll C crude oil was reflected through the distribution pattern of PAH metabolites in the bile analysed by FF. Highest concentrations were measured for 2–3-ring PAH metabolites, while 4-ring and 5-ring PAH metabolites were present in lower concentration. GC-MS results further confirmed these observations, showing C2-OH-, and C3-OH-naphthalene to be found in the highest concentrations in bile from crude oil exposed fish. The level of PAH metabolites revealed a grouping of samples according to exposure level when using multi-variate analyses (Fig. 1). A non-metric multidimensional scaling (NMDS) plot showed a distinct clustering of control and low oil concentration exposed fish, clearly separated from the cluster of medium and high oil concentration exposed individuals.
Variation among individual fish was highest for the high exposure group as revealed by the scattering. The plot indicates a dose-response effect. Moreover, PERMANOVA analysis confirmed a significant difference between the control and low oil concentration exposed fish when compared to the medium and high oil concentration samples (p = 0.002, R2 = 0.59, Pseudo-F = 14.2, DF = 11) [see Additional file 8]. This, together with the NMDS plot, indicated a dichotomy between the lower (control + low) and upper (medium + high) exposure levels in the dataset.
Gastrointestinal microbial community
Three individual fish were randomly selected from each exposure group for DNA extraction and 16S rRNA gene-based community composition analysis. First, bacterial community compositions were compared by PCR-DGGE. This first screening showed that there were large variations among individuals as shown by discrete band patterns (Fig. 2a). However, some bands were clearly present in all samples. There was also clear distinction between the band patterns of the pooled samples of lower and upper oil exposure groups (Fig. 2b). One band (#1 on Fig. 2b) in particular appeared to be more dominant in the upper exposure group while two other bands (#2 and #3 on Fig. 2b) were less dominant in comparison to the lower exposure group.
The paired-end sequencing of the 12 amplicon libraries resulted in a mean raw read number of 1,313,000 per sample [see Additional file 9]. After removal of phiX and other non-bacterial sequences we obtained a total of about 6 million assembled reads and a median of approx. 540,000 reads per sample [see Additional file 10]. After normalizing OTU abundances by common scaling based on the smallest library size (352,363 reads), and removal of OTUs present in only one sample, 994 OTUs representing > 99% of the reads in the dataset were identified. Rarefaction curve analyses on the normalized data, depicting the relationship between read number and number of detected OTUs, confirmed that the number of OTUs detected per sample was not caused by uneven sequencing depth (Fig. 3a). The difference in microbial communities between the lower and upper exposure groups was also suggested by their separation (blue vs. red tones) in an NMDS plot based on Bray-Curtis dissimilarity between the samples (Fig. 3b), as well as from the DGGE patterns observed for pooled samples (Fig. 2b).
The 12 GI microbial community samples contained between 133 and 720 OTUs and varied in diversity estimated by Shannon (H) and Inverse Simpson (1/D) indices (Table 1).
While the number of observed OTUs appeared to decrease with increasing levels of oil exposure, no significant effect on the number of OTUs was determined by ANOVA (p = 0.37) as variation within treatment was high. Likewise, no significant differences in Shannon and Inverse Simpson diversity were detected between the four treatment levels in an ANOVA analysis. When comparing the treatment levels as partitioned into a lower (control + low) and an upper (medium + high) exposure group, however, a significant difference in Shannon diversity was identified using both ANOVA (p = 0.03, F = 6.04, DF = 1) and a Wilcoxon rank-sum test between the two groups (p = 0.04, W = 31). Moreover, PERMANOVA analysis of differences in community composition between the lower and upper exposure group was performed using four different beta diversity measures [see Additional file 8]. Significance values for Bray-Curtis dissimilarity and the Jaccard index were both less than 0.05, although not below the Bonferroni corrected p-value of 0.0125 (0.05/4). Nevertheless, this analysis also indicated a shift (decrease) in the overall intestinal microbial diversity between cod exposed to lower and upper levels of crude oil.
The 994 OTUs classified into 11 phyla with highly variable relative abundances among the cod individuals and treatment levels (Fig. 4a) [see Additional file 11]. A single abundant OTU was classified by a different method than the rest of the OTUs as described in an additional file [see Additional file 12]. The phylum Proteobacteria represented 45% (263 OTUs) of the sequences in the 12 cod samples, followed by the phyla Deferribacteres (20%, 1 OTU), Bacteroidetes (11%, 186 OTUs) and Fusobacteria (8%, 45 OTUs). Proteobacteria dominated 4 out of 6 GI communities in the lower exposure group and 2 out of 6 in the upper exposure group. The remaining two lower exposure samples, 44 and 45, were dominated by Bacteroidetes or Tenericutes, respectively. Phylum Deferribacterales had the largest relative abundance in the datasets obtained from the remaining 4 upper exposure samples (49, 50, 57 and 58). Firmicutes were present in relatively high abundance among lower exposure samples and in upper sample 49. Fusobacteria showed a similar trend with having high relative abundance in only one upper exposure sample, 59.
OTU level community compositions are shown in Fig. 4b. It is worth highlighting that a single OTU (OTU01), classified as a Photobacterium species, represented 30% of the total sequences in the 12 cod samples (Fig. 4b). This OTU represented 71% of all sequences classified as Vibrionales. OTUs assigned to the orders Deferribacterales (OTU02), Vibrionales (Aliivibrio species) (OTU03) and Bacteriodales (family Porphyromonadaceae) (OTU04) were also abundant.
The microbiome of the 12 sampled cod was further classified into 18 classes and 26 orders [see Additional file 13]. Further comparison of GI microbial community composition was performed on OTU and order level.
Eight OTUs had significantly different relative abundances between the lower and upper exposure groups according to MetaStats analysis (Table 2). This was in agreement with the SIMPER analysis, identifying these as some of the OTUs contributing the most to the observed (Bray-Curtis) dissimilarity between the two exposure groups [see Additional file 14]. Interestingly, the Deferribacterales OTU02 was the only OTU with a higher relative abundance in the upper exposure group. In contrast, OTU04, a Bacteroidales belonging to the family Porphyromonadaceae, was considerably more abundant in individuals from the lower exposure group.
Normalized read counts (log2 transformed) of bacterial orders showed that in total five orders were differentially abundant in the two exposure groups (Fig. 5). Besides Deferribacterales and Bacteroides already mentioned above, Fusobacteriales, Clostridiales and Alteromonadales had a significantly lower abundance in the upper exposure groups.
A correspondence analysis (CA) based on abundance of the different orders in the 12 cod individuals revealed a tightly clustered group of medium and high exposure samples (49, 50, 57 and 58) (Fig. 6). In concurrence with the findings from MetaStats analysis, fitting of environmental parameters (bacterial orders) demonstrated that these samples were significantly (p = 0.002) and positively correlated with Deferribacterales [see Additional file 15].
The other two samples from the upper exposure group (48 and 59) clustered together with control samples and the majority of low exposure samples. This was mostly due to the high abundance of Vibrionales. This order was indeed significantly correlated (p = 0.004) with these and most of the samples from the lower exposure group. There was also a significant correlation between abundance of Fusobacteriales and the lower exposure group samples (p = 0.005). Finally, samples 44 and 45 were positioned more peripherally on the CA plot driven by their abundances of Bacteroidales (phylum Bacteroidetes, sample 44) and Mycoplasmatales (phylum Tenericutes, sample 45), respectively.
Crude oil from the Troll platform, one of the largest oil fields on the Norwegian continental shelf, was used in order to investigate changes in GI microbial community composition of Atlantic cod upon a 28-day exposure, mimicking an environmental oil spill scenario. Chemical analysis showed that Troll C crude oil, similarly to other crude oils used in exposure studies of Atlantic cod, e.g., North Sea crude, contained significant amounts of PAHs.
Overall, this study showed that a 28-day exposure to Troll C crude oil at 0.01–0.1 mg/L resulted in no fish death and no significant changes in morphological parameters. This confirmed previous findings where no changes in CI and HSI were observed after 3-weeks exposure to produced water at 0.125% (total PAH 0.102 μg/L) in Atlantic cod . Generally, a high CI and a high HSI indicate good health in fish and HSI in particular is considered the most sensitive growth indicator of fish. Nevertheless, CI can take up to several months prior to showing significant change due to environmental stressors . Sub-lethal doses of crude oil were previously demonstrated to cause declining HSI, however, such effects are not always observed .
Following exposure to petrogenic PAHs, a wide range of conjugated oxidation products are accumulated in the bile of Atlantic cod [42, 43]. In this study, an increased level of PAH metabolites, showing a dose-response pattern, was found in bile samples of oil exposed fish. Bile metabolite measurement is considered as a sensitive tool for assessing exposure to oil and also produced water containing PAHs . Several laboratory and field studies documented significantly elevated bile metabolite levels even at very low exposure concentrations similarly to the experiment presented here. PERMANOVA and NMDS analysis of GC-MS measured PAH metabolite concentrations in fish bile showed a dichotomy between the two lower (0.0 and 0.01 mg/L oil) and the two upper (0.05 and 0.1 mg/L oil) exposure levels. Based on these results, statistical analyses performed on the GI microbial community was done on pooled data, i.e., on data combined into lower and upper exposure groups.
OTU richness with numbers ranging from 133 to 720 OTUs identified in the 12 cod microbiomes were similar to that observed in analysis of 11 specimens of wild-caught Atlantic cod, where the number of OTUs per sample ranged from 40 to 228 . A significantly higher number of OTUs (1850) were detected in 20 specimens of juvenile Atlantic salmon (salmon parr) at a similar sequencing depth . Despite a relatively long time (28 days) spent under laboratory conditions, Shannon and Inverse Simpson indices varied greatly among individuals. Given the assumption that identical environmental conditions and feeding routines are expected to homogenize GI microbial communities this result was perhaps unexpected. Nevertheless, the ranges were somewhat narrower (Shannon: 0.54–2.30 and Inverse Simpson: 1.31–5.14) when compared to 11 specimens of wild-caught Atlantic cod kept under laboratory conditions for 12 days (Shannon: 0.30–3.07 and Inverse Simpson: 1.09–11.18). Large inter-individual variation in the intestinal microbial composition has been observed in other studies as well [44, 45] including GI microbial communities of a group of identically reared sibling zebrafishes . Based on results emerging from next-generation sequencing studies of fish GI tract microbiota, the idea that each individual harbors a unique microbial ecosystem is being constantly reinforced . Such high inter-individual variability has already raised the question of experimental reproducibility of exposure studies, which may be affected by the influence of highly variable fish associated microflora . However, implications of this variability in terms of functional variability of the GI tract microbiome and host-microbiome interactions is yet unknown.
In the study presented here, the microbial contents of the entire GI tract were used from fishes that were starved 7 days prior to sampling. It is well known that different alimentary components harbor distinct microbial communities  and that starvation influences the composition of GI tract microbiota . In future studies, it is preferable to collect intestinal contents based on a prior visual assessment of the alimentary anatomy.
On the phylum level, pentachlorophenol (PCP) exposure of goldfish for 28 days resulted mainly in increase in Bacteroidetes abundance and reduction of Bacteroidetes/Firmicutes ratio . We did not observe such effect of oil on Atlantic cod in this study. Here, both Bacteroidetes and Firmicutes were present in lower relative abundances among upper exposure level individuals. Another recent study, which examined the effects of radio-labeled graphene nanoparticles of different size in adult zebrafish, also found that microbial communities were affected, and the size of the graphene nanoparticles seemed to be an important factor in determining the composition of GI tract microbiota . While bacterial community composition of control fish (n = 4) and large nanoparticle exposed fish (n = 4) were similar, significant differences were observed in the gut microbiota composition of the small nanoparticle exposed fish (n = 4) based on NMDS analysis.
The intestinal bacterial community in our three control samples resembles that found in 11 cod individuals from the inner Oslo fjord , as it is dominated by Vibrionales and Bacteroidales, followed by Clostridiales and Desulfovibrionales. The large proportion of Vibrionales (median: 55%) in our samples is close to what is observed by Star et al. (median: 50%)  and is also in agreement with findings in other marine carnivores . Further, members of the most abundant orders have previously been reported in Atlantic cod [49, 52,53,54,55]. Although the comparison with the Oslo fjord samples is based on a small sample sizes, the similarity between these communities indicates that our laboratory setting did not alter the species composition drastically compared to those found in wild populations.
Deferribacterales is a member of the recently discovered phylum Deferribacteres, and bacteria belonging to this group are considered acetoclastic iron-reducers. Iron is a nutrient that often limits bacterial growth, and animals typically try to reduce biologically available iron in order to avoid bacterial infections. Interestingly, it was found earlier, that exposure of juvenile Atlantic cod to HDF200 base oil (containing PAHs) for 30 days, resulted in downregulation of serotransferrin, a glycoprotein involved in iron-binding . Lower levels of serotransferrin could result in higher levels of biologically available iron affecting microbial community composition in the cod gut.
The phylum Deferribacteres is present -albeit at low levels of relative abundance- in both the Oslo fjord cod  and our control and low-level individuals. Nonetheless, here we find that it is the sole order that significantly increased in relative abundance with increasing oil exposure levels. Deferribacterales has a wide environmental distribution and members have been found in hydrothermal vent communities , the gut microbiota of the hydrothermal shrimp Rimicaris exoculate  and in rodents . Yet interestingly this order has also been found in high abundance in subsea petroleum reservoirs . Moreover, it also occurs in relatively high abundance (19.2%) in biodegraded crude oil samples while it is absent in non-biodegraded samples from the same site . This implies that Deferribacteres species thrive under oil exposed conditions. A recent metagenomic study of succession in a petroleum reservoir detected gene clusters encoding the pathway for anaerobic degradation of monoaromatic compounds in bins assigned to the genus Flexistipes (Deferribacterales) . This could suggest that other genera in the Deferribacterales might be capable of breaking down aromatic hydrocarbons as well. However, since none of the described species within the Deferribacteres are implicated in the active breakdown of complex recalcitrant hydrocarbons , it is of interest to see the results of such experiments. Until then, we have to be cautious to interpret that our findings are due to catabolism of aromatic compounds by Deferribacterales species found in the Atlantic cod gut microbiome. Nonetheless, the Deferribacterales metabolize small molecules that are also abundant in the gut, such as monosaccharides, amino acids, short chain fatty acids in the presence of a suitable electron donor such as iron and nitrate. Therefore, when PAH exposure triggers higher iron concentrations in the cod gut, these bacteria could metabolize more of the small molecules and outcompete others. Our results provide another line of evidence suggesting that specific members of this order thrive in an environment exposed to oil, regardless of the presence of extant, complex microbial communities.
It is recognized that the GI microbiome plays an essential role in gut health, immune function and nutrient uptake in fish. Although several studies have been investigating the intestinal microflora of various fish species, our current knowledge regarding the roles of gut microbiota in ecotoxicology remains limited. The exposure set-up in our experiment represents a realistic situation of chronic oil pollution. To the best of our knowledge, this is the first study demonstrating how exposure to sub-lethal concentrations of crude oil alters the intestinal microbiome of Atlantic cod. Despite the low number of samples, we identified 8 OTUs which were present in different relative abundances among lower and upper exposure groups. OTUs classified as Porphyromonadaceae, Rikenella, Ruminococcaceae, Alistipes and Clostridiales decreased in relative abundance. The only OTU which showed increased relative abundance at higher oil exposure was classified as Deferribacterales. Hence, Deferribacterales appear to be a potentially interesting microbial indicator of oil exposure in Atlantic cod. Further research with a larger sample number and functional analysis is necessary, in order to understand the implications of such shift and the potential changes of metabolic functions of intestinal microorganisms in response to environmental pollutants.
Analysis of variance
Continuous flow system
Denaturing gradient gel electrophoresis
European Nucleotide Archive
Environmental Protection Agency
Fixed wavelength fluorescence
Gas chromatography mass spectrometry
Non-metric multidimensional scaling
Operational taxonomic unit
Polycyclic aromatic hydrocarbons
Polymerase chain reaction
Permutational multivariate analysis of variance
Pyrene fluorescence equivalents
Ribosomal Ribonucleic acid
Single ion mode
International Arctic Petroleum Cooperation. Barents Sea Scenarios. In: Bourmistrov A, Mellemvik F, Bambulyak A, Gudmestad O, Overland I, Zolotukhin A, editors. Routledge Studies in Environmental Policies. Abingdon: Routledge; 2015.
Sundt RC, Ruus A, Jonsson H, Skarphéðinsdóttir H, Meier S, Grung M, Beyer J, Pampanin DM. Biomarker responses in Atlantic cod (Gadus morhua) exposed to produced water from a North Sea oil field: laboratory and field assessments. Mar Pol Bul. 2012;64(1):144–52.
Pampanin DM, Sydnes MO. Polycyclic aromatic hydrocarbons a constituent of petroleum: presence and influence in the aquatic environment. In: Vladimir K, Kolesnikov A, editors. Hydrocarbon. Rijeka: InTech; 2013. https://doi.org/10.5772/48176.
Le Goff J. Carcinogenicity of Petrogenic PAHs. In: Pampanin DM, Sydnes MO, editors. Petrogenic polycyclic aromatic hydrocarbons in the aquatic environment: analysis, Synthesis, Toxicity and Environmental Impact. Bentham Science Publishers, Sharjah, UAE; 2017. p. 65–110.
Fernandes D, Marqueno A, Porte C, Sole M. PAH metabolites in fish and invertebrates: analysis and endocrine disruptive potential. In: Pampanin DM, Sydnes MO, editors. Petrogenic polycyclic aromatic hydrocarbons in the aquatic environment: analysis, Synthesis, Toxicity and Environmental Impact: Bentham Science Publishers; 2017. p. 111–34.
Aas E, Baussant T, Balk B, Liewenborg B, Andersen OK. PAH metabolites in bile, cytochrome P4501A and DNA adducts as environmental risk parameters for chronic oil exposure: a laboratory experiment with Atlantic cod. Aquat Toxicol. 2000;51(2):241–58.
Pampanin DM, Brooks SJ, Grøsvik BE, Le Goff J, Meier S, Sydnes MO. DNA adducts in marine fish as biological marker of genotoxicity in environmental monitoring: the way forward. Mar Environ Res. 2017;124:49–62.
Stegeman JJ, Lech JJ. Cytochrome P-450 monooxygenase systems in aquatic species: carcinogen metabolism and biomarkers for carcinogen and pollutant exposure. Environ Health Perspect. 1991;90:101–9.
Ohnishi S, Kawanishi S. Double base lesions of DNA by a metabolite of carcinogenic benzo[a]pyrene. Biochem Biophys Res Commun. 2002;290(2):778–82.
Kleinow KM, Nichols JW, Hayton WL, McKim JM, Barron MG. Toxicokinetics in Fishes. In: Di Giulio RT, Hinton DE, editors. The Toxicology of Fishes. Boca Raton: CRC Press, Taylor & Francis Group; 2008. p. 56–134.
Bakke J, Struble C, Gustafsson JA, Gustafsson B. Catabolism of premercapturic acid pathway metabolites of naphthalene to naphthols and methylthio-containing metabolites in rats. Proc Natl Acad Sci U S A. 1985;82(3):668–71.
Claus SP, Guillou H, Ellero-Simatos S. The gut microbiota: a major player in the toxicity of environmental pollutants? NPJ Biofilms and Microbiomes. 2016;2:16003.
Zhang L, Nichols RG, Correll J, Murray IA, Tanaka N, Smith PB, Hubbard TD, Sebastian A, Albert I, Hatzakis E, Gonzalez FJ, Perdew GH, Patterson AD. Persistent organic pollutants modify gut microbiota–host metabolic homeostasis in mice through aryl hydrocarbon receptor activation. Environ Health Perspect. 2015;123(7):679–88.
Maurice CF, Haiser HJ, Turnbaugh PJ. Xenobiotics shape the physiology and gene expression of the active human gut microbiome. Cell. 2013;152(1–2):39–50.
Ringø E, Strøm E, Tabachek J-A. Intestinal microflora of salmonids: a review. Aquac Res. 2008;26(10):773–89.
Star B, THA H, Jentoft S, Jakobsen KS. Next generation sequencing shows high variation of the intestinal microbial species composition in Atlantic cod caught at a single location. BMC Microbiol. 2013;13:248.
Smith CC, Snowberg LK, Gregory Caporaso J, Knight R, Bolnick D. Dietary input of microbes and host genetic variation shape among-population differences in stickleback gut microbiota. The ISME J. 2015;9(11):2515–26.
Kan H, Zhao F, Zhang XX, Ren H, Gao S. Correlations of gut microbial community shift with hepatic damage and growth inhibition of Carassius auratus induced by pentachlorophenol exposure. Environ Sci Technol. 2015;49(19):11894–902.
Sanni S, Øysæd KB, Høivangli V, Gaudebert B. A continuous flow system (CFS) for chronic exposure of aquatic organisms. Mar Environ Res. 1998;46(1–5):97–101.
Delgado IAK. Biomarker responses in Atlantic cod (Gadus morhua) exposed to PAHs: data treatment, data interpretation and communication of results. M.Sc. Thesis. 2016. Accessed from: http://hdl.handle.net/11250/2410623. Accessed 26 Mar 2018.
Goede RW, Barton BA. Organismic indices and an autopsy-based assessment as indicators of health and condition of fish. In: Adams SM, editor. Biological indicators of stress in fish. Bethesda: American Fisheries Society; 1990. p. 93–108.
Nash RDM, Valencia AH, Geffen AJ. The origin of Fulton’s condition factor—setting the record straight. Fisheries. 2006;31:236–8.
Aas E, Beyer J, Goksøyr A. Fixed wavelength fluorescence (FF) of bile as a monitoring tool for polyaromatic hydrocarbon exposure in fish: an evaluation of compound specificity, inner filter effect and signal interpretation. Biomarkers. 2000;5(1):9–23.
Jonsson G, Beyer J, Wells D, Ariese F. The application of HPLC-F and GC-MS to the analysis of selected hydroxy polycyclic hydrocarbons in two certified fish bile reference materials. J Environ Monit. 2003;5(3):513–20.
Jonsson G, Taban IC, Jørgensen KB, Sundt RC. Quantitative determination of de-conjugated chrysene metabolites in fish bile by HPLC-fluorescence and GC-MS. Chemosphere. 2004;54(8):1085–7.
Brakstad OG, Bonaunet K. Biodegradation of petroleum hydrocarbons in seawater at low temperatures and bacterial communities associated with degradation. Biodegradation. 2006;17(1):71–82.
FastQC: A quality control tool for high throughput sequence data. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ Accessed 26 Mar 2018.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17(1):10–2.
Li H, Durbin R. Fast and accurate short read alignment with burrows–wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Seqkt: Toolkit for processing sequences in FASTA/Q formats. https://github.com/lh3/seqtk. Accessed 26 Mar 2018.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Env Microbiol. 2009;75(23):7537–41.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Racine JS. R studio: a platform independent IDE for R and SWEAE. J Appl Econ. 2012;27:167–72.
The R Project for Statistical Computing. www.r-project.org. Accessed 26 Mar 2018.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217.
Wickham H. ggplot2: Elegant Graphics for Data Analysis. 1st ed. New York: Springer-Verlag; 2009.
McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10(4):e1003531.
Phyloseq. https://github.com/joey711/phyloseq/issues/143. Accessed 26 Mar 2018.
Oksanen J. Multivariate Analysis of Ecological Communities in R: vegan tutorial. 2015. http://cc.oulu.fi/~jarioksa/opetus/metodi/vegantutor.pdf. Accessed 26 Mar 2018.
White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009;5(4):e1000352.
Hoque M, Yusoff F, Law A, Syed MA. Effect of hydrogen sulphide on liver somatic index and Fulton's condition factor in Mystus nemurus. J Fish Biol. 1998;52(1):23–30.
Holth TF, Eidsvoll DP, Farmen E, Sanders MB, Martínez-Gómez C, Budzinski H, Burgeot T, Guilhermino L, Hylland K. Effects of water accommodated fractions of crude oils and diesel on a suite of biomarkers in Atlantic cod (Gadus morhua). Aquat Toxicol. 2014;154:240–52.
Sanni S, Björkblom C, Jonsson H, Godal BF, Liewenborg B, Lyng E, Pampanin DM. I: biomarker quantification in fish exposed to crude oil as input to species sensitivity distributions and threshold values for environmental monitoring. Mar Environ Res. 2016;125:10–24.
Dehler CE, Secombes CJ, Martin SA. Environmental and physiological factors shape the gut microbiota of Atlantic salmon parr (Salmo salar L.). Aquaculture. 2017;467:149–57.
Lyons PP, Turnbull JF, Dawson KA, Crumlish M. Exploring the microbial diversity of the distal intestinal lumen and mucosa of farmed rainbow trout Oncorhynchus mykiss (Walbaum) using next generation sequencing (NGS). Aquac Res. 2015;48(1):77–91.
Melancon E, Gomez De La Torre Canny S, Sichel S, Kelly M, Wiles TJ, Rawls JF, Eisen JS, Guillemin K. Best practices for germ-free derivation and gnotobiotic zebrafish husbandry. Chapter 3. In: Detrich III HW, Westerfield M, Zon LI, editors. The Zebrafish: Disease Models and Chemical Screens, Methods in Cell Biology, vol. 138. 4th ed. Cambridge: Academic Press; 2017.
Clements KD, Angert ER, Montgomery WL, Choat JH. Intestinal microbiota in fishes: what's known and what's not. Mol Ecol. 2014;23(8):1891–8.
Vatsos IN. Standardizing the microbiota of fish used in research. Lab Anim. 2016;51(4):353–64.
Dhanasiri AK, Brunvold L, Brinchmann MF, Korsnes K, Bergh Ø, Kiron V. Changes in the intestinal microbiota of wild Atlantic cod Gadus morhua L. Upon Captive Rearing. Microb Ecol. 2011;61(1):20–30.
Lu K, Dong S, Petersen EJ, Niu J, Chang X, Wang P, Lin S, Gao S, Mao L. Biological uptake, distribution, and depuration of radio-labeled graphene in adult zebrafish: effects of graphene size and natural organic matter. ACS Nano. 2017;11(3):2872–85.
Sullam KE, Essinger SD, Lozupone CA, O’Connor MP, Rosen GL, Knight R, Kilham SS, Russell JA. Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. Mol Ecol. 2014;21(13):3363–78.
Ringø E, Sperstad S, Myklebust R, Refstie S, Krogdahl Å. Characterisation of the microbiota associated with intestine of Atlantic cod (Gadus morhua L.): the effect of fish meal, standard soybean meal and a bioprocessed soybean meal. Aquaculture. 2006;261(3):829–41.
Brunvold L, Sandaa R-A, Mikkelsen H, Welde E, Bleie H, Bergh Ø. Characterisation of bacterial communities associated with early stages of intensively reared cod (Gadus morhua) using denaturing gradient gel electrophoresis (DGGE). Aquaculture. 2007;272(1–4):319–27.
Fjellheim AJ, Playfoot KJ, Skjermo J, Vadstein O. Vibrionaceae dominates the microflora antagonistic towards Listonella anguillarum in the intestine of cultured Atlantic cod (Gadus morhua L.) larvae. Aquaculture. 2007;269(1–4):98–106.
Reid HI, Treasurer JW, Adam B, Birkbeck TH. Analysis of bacterial populations in the gut of developing cod larvae and identification of Vibrio logei, Vibrio anguillarum and Vibrio splendidus as pathogens of cod larvae. Aquaculture. 2009;288(1–2):36–43.
Pampanin DM, Larssen E, Oysæd KB, Sundt RC, Sydnes MO. Study of the bile proteome of Atlantic cod (Gadus morhua): multi-biological markers of exposure to polycyclic aromatic hydrocarbons. Mar Environ Res. 2014;101:161–8.
Miroshnichenko ML, Slobodkin AI, Kostrikina NA, L'Haridon S, Nercessian O, Spring S, Stackebrandt E, Bonch-Osmolovskaya EA, Jeanthon C. Deferribacter abyssi sp. nov., an anaerobic thermophile from deep-sea hydrothermal vents of the mid-Atlantic ridge. Int J Syst Evol Microbiol. 2003;53:1637–41.
Cowart DA, Durand L, Cambon-Bonavita MA, Arnaud-Haond S. Investigation of bacterial communities within the digestive organs of the hydrothermal vent shrimp Rimicaris exoculata provide insights into holobiont geographic clustering. PLoS One. 2017;12(3):e0172543.
Berry D, Kuzyk O, Rauch I, Heider S, Schwab C, Hainzl E, Decker T, Müller M, Strobl B, Schleper C, Urich T, Wagner M, Kenner L, Loy A. Intestinal microbiota signatures associated with inflammation history in mice experiencing recurring colitis. Front Microbiol. 2015;6:1408.
Gittel A, Kofoed MV, Sørensen KB, Ingvorsen K, Schramm A. Succession of Deferribacteres and Epsilonproteobacteria through a nitrate-treated high-temperature oil production facility. Syst Appl Microbiol. 2012;35(3):165–74.
Silva TR, Verde LCL, Santos Neto EV, Oliveira VM. Diversity analyses of microbial communities in petroleum samples from Brazilian oil fields. Int Biodeterior Biodegradation. 2013;81:57–70.
Vigneron A, Alsop EB, Lomans BP, Kyrpides NC, Head IM, Tsesmetzis N. Succession in the petroleum reservoir microbiome through an oil field production lifecycle. The ISME J. 2017;1(9):2141–54.
Alauzet C, Jumas-Bilak E. The phylum Deferribacteres and the genus Caldithrix. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes. 4th ed. Berlin Heidelberg: Springer-Verlag; 2014. p. 595–611.
The authors would like to acknowledge Karianne Skogland Enerstvedt, Mari Mæland Nilsen, Azin Bayat, and Mads Busvold for their help with the exposure set-up and sampling. The sequencing service was provided by the Norwegian Sequencing Centre (www.sequencing.uio.no), a national technology platform hosted by the University of Oslo and supported by the “Functional Genomics” and “Infrastructure” programs of the Research Council of Norway and the South-eastern Regional Health Authorities.
The Research Council of Norway, PETROMAKS 2 (project # 229153), is gratefully acknowledged for financial support. Additional funding was obtained from the research program: Bioactive, University of Stavanger, Stavanger, Norway. The funding bodies had no role in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The dataset (raw sequencing data) analysed during the current study are available in the European Nucleotide Archive (ENA) repository, http://www.ebi.ac.uk/ena/data/view/PRJEB21667.
All supporting data generated during this study are included in this published article and its additional files.
Ethics approval and consent to participate
The experiment was approved by the Norwegian Food Safety Authority, FOTS (case number: 7262) and was conducted in accordance with the European Convention for the protection of vertebrate animals (http://conventions.coe.int/treaty/en/treaties/html/123.htm) used for experimental and other scientific purposes.
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.
Figure S1. Continuous flow system supplying dispersed crude oil into the exposure tanks at different levels. (PDF 167 kb)
Table S1. Primer and barcode sequences. (XLSX 9 kb)
Table S2. Composition of Troll C crude oil. (XLSX 12 kb)
Table S3. Polycyclic aromatic hydrocarbon (PAH) content of oil dispersion in header tank. (XLSX 12 kb)
Table S4. Condition indices (CI) and hepatosomatic indices (HSI) of oil exposed Atlantic cod specimens. (XLSX 9 kb)
Table S5. Level of polycyclic aromatic hydrocarbon (PAH) metabolites in fish bile samples measured by fixed fluorescence (FF). (XLSX 8 kb)
Table S6. Level of polycyclic aromatic hydrocarbon (PAH) metabolites in fish bile samples measured by gas chromatography mass spectrometry (GC-MS). (XLSX 13 kb)
Table S7. Result of PERMANOVA analysis. (XLSX 10 kb)
Table S8. Library sizes. (XLSX 9 kb)
Table S9. Assembled reads. (XLSX 9 kb)
Table S10. Taxa levels and abundances. (XLSX 25 kb)
Supplementary text describing the classification method used for an unclassified highly abundant OTU and phylogenetic tree. (DOCX 380 kb)
Table S11. OTU table for all samples. (XLSX 222 kb)
Table S12. Result of SIMPER analysis. (XLSX 16 kb)
Table S13. Result of envfit analysis. (XLSX 9 kb)