- Open Access
Holobionts and ecological speciation: the intestinal microbiota of lake whitefish species pairs
Microbiomevolume 6, Article number: 47 (2018)
It is well established that symbionts have considerable impact on their host, yet the investigation of the possible role of the holobiont in the host’s speciation process is still in its infancy. In this study, we compared the intestinal microbiota among five sympatric pairs of dwarf (limnetic) and normal (benthic) lake whitefish Coregonus clupeaformis representing a continuum in the early stage of ecological speciation. We sequenced the 16s rRNA gene V3-V4 regions of the intestinal microbiota present in a total of 108 wild sympatric dwarf and normal whitefish as well as the water bacterial community from five lakes to (i) test for differences between the whitefish intestinal microbiota and the water bacterial community and (ii) test for parallelism in the intestinal microbiota of dwarf and normal whitefish.
The water bacterial community was distinct from the intestinal microbiota, indicating that intestinal microbiota did not reflect the environment, but rather the intrinsic properties of the host microbiota. Our results revealed a strong influence of the host (dwarf or normal) on the intestinal microbiota with pronounced conservation of the core intestinal microbiota (mean ~ 44% of shared genera). However, no clear evidence for parallelism was observed, whereby non-parallel differences between dwarf and normal whitefish were observed in three of the lakes while similar taxonomic composition was observed for the two other species pairs.
This absence of parallelism across dwarf vs. normal whitefish microbiota highlighted the complexity of the holobiont and suggests that the direction of selection could be different between the host and its microbiota.
Earth is dominated by approximately 1030 microbial cells , which is two- or three-fold more than the number of plant and animal cells combined . Therefore, it is important to consider that animal and plant evolution has and continues to occur in the presence of microbiota, which have either parasitic, mutualistic, or commensal interaction with a host . The ubiquity and importance of the microbiota is supported by its influence on host development, immunity, metabolism, behavior, and numerous other processes including speciation [4,5,6,7,8,9,10,11,12]. The host (animal or plant) and their microbiota are referred to as a “holobiont” [10, 13,14,15], which represents a unique biological entity evolving through selection, drift, mutation, and migration .
The concept of holobiont offers a new angle for the study of adaptive divergence ultimately leading to speciation. For instance, the role of microbiota on pre-zygotic isolation has recently been documented . Moreover, the host’s visual, auditory, and chemosensory signals implicated in mate choice could be influenced by its microbiota [18,19,20,21,22]. Also, host populations sharing similar environment or diet have been shown to share similar microbiomes, known as a “socially shared microbiome” . The socially shared microbiome could recognize specific signals of the host population and thus influence its evolution in ways that are microbe-specific and microbe-assisted, which may lead to post-zygotic isolation .
The intestinal microbiota could be particularly prone to playing a key role in the process of population divergence and speciation given its broad array of functional impacts on its host . The involvement of the intestinal microbiota in organismal functions comprises nutrition [24, 25], toxicity resistance , energy metabolism [9, 27, 28], morphology , and behavior [5, 8, 30, 31]. On the other hand, the intestinal microbiota can also promote host phenotypic plasticity, which may contribute to adaptation. For example, new intestinal microbiota genes can be acquired from the environment through acquisition of new bacteria [32, 33]. The intestinal microbiota can also adapt in response to variation in the host’s physiological and environmental conditions . Moreover, the short generation time of the intestinal microbiota and the horizontal transfer of genes can favor rapid microbiota evolution [35, 36].
While there are now a plethora of studies that have documented the positive influence of holobionts on hosts, including humans, relatively few studies have focused on fish microbiota in the wild even though they represent around 50% of the total vertebrate diversity [37, 38]. To date, about 20 studies have investigated fish intestinal microbiota in the wild (e.g., [39,40,41,42,43]). Of these, very few concerned speciation and to our knowledge, none analyzed specifically the adherent bacteria present in the fish epithelial mucosa [44,45,46,47,48,49]. Adherent bacteria are of particular interest because they may interact more closely with their host than bacteria present in the alimentary bolus .
Lake whitefish (Coregonus clupeaformis) comprises sympatric species pairs referred to as dwarf and normal whitefish that are found in five lakes of the St. John River drainage in the province of Québec, Canada, and in Maine, USA. A relatively recent period of post-glacial adaptive radiation occurred approximately 12,000 years before present (YBP), leading to parallel phenotypic and ecological divergence in different lakes of the dwarf whitefish derived from the ancestral normal whitefish . Dwarf and normal whitefish are partially reproductively isolated in each lake , differ in genetically based morphological, physiological, behavioral, ecological, and life history traits [52,53,54,55,56] and occupy the limnetic and benthic habitat, respectively. Dwarf and normal whitefish also differ in trophic niche, where dwarf whitefish (and limnetic whitefish in general) feed almost exclusively on zooplankton [57, 58] and normal whitefish are more generalist and feed on more diverse prey items including zoobenthos, molluscs, and fish prey [50, 58].
In this study, we investigate the within- and between-lake variation in the intestinal microbiota among these five sympatric pairs of dwarf and normal whitefish, representing a continuum in the early stage of ecological speciation. We sequenced the 16S rRNA gene of adherent bacteria present in the intestinal tissue and in order to test for differences between intestinal microbiota of dwarf and normal whitefish pairs. We chose adherent microbiota present on intestinal tissues because this microbiota may be more involved in host-microbiota interactions. In parallel, we also sequenced the 16S rRNA gene of water bacterial communities from the five lakes in order to test the association between the water bacterial community and the whitefish intestinal microbiota. Ultimately, our main goal was to test for the occurrence of parallelism in the microbiota of sympatric dwarf and normal whitefish across different environments, where evidence for parallelism would provide strong indirect evidence for the role of natural selection in shaping host microbiota.
Lake whitefish (44 dwarf and 64 normal fish) were sampled with gill nets from Cliff Lake, Indian Pond, and Webster Lake in Maine, USA, in June 2013, and from East and Témiscouata lakes in Québec, Canada, during summer 2013, from May to July (Table 1). Fish were dissected in the field in sterile conditions. The ventral belly surface was rinsed with 70% ethanol, and non-disposable tools were rinsed with ethanol and flamed over a blowtorch between samples. The intestine was cut at the hindgut level (posterior part of the intestine), and the digesta were aseptically removed. Then, the intestine was cut at the foregut level (anterior part of the intestine), removed from the peritoneal cavity, and clamped on both extremities in order to isolate the adherent bacteria in the laboratory. The clamped intestines were individually stored in sterile cryotubes and flash-frozen in liquid nitrogen. Water samples (2 L) were collected in each lake at four depths (at the top of the water column, at 5, 10, and 15 m corresponding to 1 m above the lake bottom) with a Niskin© (General Oceanics). Water samples were filtered first with a 3.0-μm mesh, followed by a 0.22-μm nitrocellulose membrane using a peristaltic pump (Cole-Parmer: Masterflex L/S Modular Drive). The 0.22-μm membranes were placed into cryotubes and flash-frozen with liquid nitrogen. All samples were transported to the laboratory and kept at − 80 °C until further processing.
DNA extraction, amplification, and sequencing of intestinal bacteria
Adherent bacterial DNA from the intestinal segment was isolated by rinsing the interior of the intestines three times with 3 ml of sterile 0.9% saline  and extracted using a modification of the QIAmp© Fast DNA stool mini kit (QIAGEN). In order to ensure efficient lysis of Gram-positive bacteria, temperature and digestion time were increased during the incubation steps. Moreover, to maximize DNA extraction, the volume of supernatant and all of the products used with the supernatant (Proteinase K, Buffer AL, and ethanol 100%) were doubled. Thus, 1200 μl were transferred into the column (in two subsequent steps) and bacterial DNA was eluted from the column with 100 μl of ultrapure water (DEPC-treated Water Ambion®). Bacterial DNA from the water samples was also extracted using a modified QIAmp© Fast DNA stool mini kit (QIAGEN) protocol. The 0.22-μm membranes were transferred with a 1-ml InhibitEX buffer to bead beating tubes (Mobio), incubated overnight at 50 °C, and then vortexed for 1 h. The same modified protocol used for the adherent bacterial DNA was used. In order to test the sterility during the extraction manipulation, seven blank extractions were done with buffer only. Moreover, the same extraction kit was used between fish microbiota and water bacterial community in order to avoid bias during extraction. Extracted DNA was quantified with a Nanodrop (Thermo Scientific) and stored at − 20 °C until use.
The partial DNA fragments of bacterial 16S rRNA genes were amplified by touchdown PCR for adherent bacterial DNA. Touchdown PCR is the optimal method to avoid eukaryotic contamination, potentially due to cross amplification with host DNA [60, 61]. A region ~ 250 bp in the 16S rRNA gene, covering the V3–V4 regions, was selected to construct the community library using specific primers with Illumina barcoded adapters Bakt_341F-long and Bakt_805R-long  in a dual indexed PCR approach. The touchdown PCR of adherent bacterial DNA used 25 μl of NEBNext Q5 Hot Start Hifi PCR Master Mix, 1 μl (0.2 μM) of each specific primer, 15 μl of sterile nuclease-free water, and 8 μl of DNA (around 170 ng/μL). The PCR program consisted of an initial denaturation step at 98 °C for 30 s, followed by 20 cycles at 98 °C for 10 s, 67–62 °C (touchdown PCR annealing step) for 30 s, and 72 °C for 45 s. After the initial touchdown PCR cycles, an additional 15 cycles were run at 98 °C for 10 s (denaturation), 62 °C for 30 s (annealing) and 72 °C for 45 s (extension), and a final extension of 72 °C for 5 min.
The PCR amplification for water bacterial DNA comprised a 50-μl PCR amplification mix containing 25 μl of NEBNext Q5 Hot Start Hifi PCR Master Mix, 1 μl (0.2 μM) of each specific primer, 21 μl of sterile nuclease-free water, and 2 μl of water bacterial DNA (around 5 ng/μL). The PCR program consisted of an initial denaturation step at 98 °C for 30 s, followed by 30 cycles, with 1 cycle at 98 °C for 10 s (denaturation), 56 °C for 30 s (annealing) and 72 °C for 45 s (extension), and a final extension of 72 °C for 5 min. Negative and positive controls were included for all PCRs. All the PCR results, including the negative controls, were purified using the AMPure bead calibration method. The purified samples were quantified using a fluorometric kit (QuantIT PicoGreen; Invitrogen), pooled in equimolar amounts, and sequenced paired-end using Illumina MiSeq Bakt_341F-long and Bakt_805R-long at the Plateforme d’Analyses Génomiques (IBIS, Université Laval, Québec, Canada). To prevent focusing, template building, and phasing problems due to the sequencing of low-diversity libraries such as 16S rRNA amplicons, 50% PhiX genome was spiked in the pooled library.
Raw forward and reverse reads were quality trimmed, assembled into contigs for each sample, and classified using Mothur v.1.36.0 [63, 64]. Contigs were quality trimmed with the following criteria: (i) when aligning paired ends, a maximum of two mismatches were allowed; (ii) ambiguous bases were excluded; (iii) homopolymers of more than 8 bp were removed; (iv) sequences with lengths less than 400 bp and greater than 450 bp were removed; (v) sequences from chloroplasts, mitochondria, and non-bacterial were removed; and (vi) chimeric sequences were removed using the UCHIME algorithm . Moreover, the database SILVA was used for the alignment and the database RDP (v9) was used to classify the sequences with a 0.03 cutoff level. The Good’s coverage index, Shannon index, inverse Simpson diversity, and weighted UniFrac tests were estimated with Mothur. The Good’s coverage index estimates the quality of the sequencing depth whereas alpha diversity (diversity within the samples) was estimated with the inverse Simpson index and the Shannon index. Beta diversity (diversity between samples) was calculated using a weighted UniFrac test , which was performed using thetayc distance.
A matrix containing the number of bacterial sequences was constructed for each genus in each fish sample from the two Mothur taxonomy files (stability.an.shared and stability.an.cons.taxonomy). Therefore, OTUs (operational taxonomic units) with the same taxonomy were merged. This genus-merged matrix was used to perform the taxonomic composition analysis at the phylum and genus level, the principal coordinate analyses (PCoA), the permutational analysis of variance (PERMANOVA), the Metastats analysis, and the network analysis. Moreover, to determine if there was a significant difference at the alpha diversity level between species within and among lakes, we used a generalized linear model (GLM) with a Gaussian family followed by an ANOVA. In order to build the PCoAs, a Jaccard distance matrix was made from the genus-merged matrix after Hellinger transformation using the vegan package  in R (R Core Team 2016). The PERMANOVA analysis (number of permutations = 10,000) was also performed with the vegan package in R to test the species effects, the lake effects, and their interaction. The METASTATS software with standard parameters was also used (p ≤ 0.05 and number of permutations = 1000) to detect differential abundance of bacteria at the genus level between dwarf and normal whitefish . Network analyses, based on a Spearman’s correlation matrix, were performed to document the interaction between dwarf and normal whitefish microbiota. The Spearman’s correlation matrix was calculated using R on the Hellinger transformed matrix. Moreover, P values and Bonferroni corrections were calculated for Spearman’s correlations for each sample. Then, the different networks were visualized using Cytoscape version 3.2.1, a software for visualizing networks . Finally, PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States, version 1.0.0) was used to predict putative functions for the whitefish microbiota based on the 16S rRNA sequence dataset . To this end, our OTU data was assigned against the Greengenes database (released August 2013) and we used the Mothur command “make.biom” to obtain a data file compatible with PICRUSt.
A total of 1,603,342 sequences were obtained after trimming for the entire dataset composed of 108 whitefish intestinal microbiota (44 dwarf and 64 normal whitefish) and 36 bacterial water samples (Additional file 1: Table S1). Among these sequences, 24,308 different operational taxonomic units (OTUs) were identified with a 97% identity threshold, representing 544 genera. The average Good’s coverage estimation, used to estimate the quality of the sequencing depth, was 99% ± 2% of coverage index.
Very few sequences were obtained from the five PCR negative controls (Additional file 2: Table S2). Although there were no bands after PCR amplification, 95 sequences in total were obtained from the five PCR negative controls, representing 0.006% of the total dataset. Sixty-one different species were identified with a range of 1–11 reads per bacterial species. Some of these sequences represented bacteria that are typically associated with fish, seawater, or freshwater environments, but also with fish pathogens (Additional file 2: Table S2). None were associated to humans or to the laboratory environment. This suggests that contamination was very low, but not completely absent, as typically observed in similar studies [71,72,73].
Whitefish intestinal microbiota vs. water bacterial communities
Highly different communities at the genus level were observed with weighted UniFrac and PERMANOVA tests between the water bacterial community and whitefish microbiota within each lake and among the lakes (Table 2). Moreover, water bacterial communities as well as dwarf and normal whitefish microbiota had distinct dominant phyla composition (Fig. 1a). The water bacterial community was composed of Proteobacteria (38.7%), Actinobacteria (33.5%), Bacteroidetes (10.6%), Verrucomicrobia (4.4%), OD1 (2.0), and Firmicutes (1.9%). The five most abundant phyla of dwarf intestinal microbiota were Proteobacteria (40.6%), Firmicutes (17.8%), Actinobacteria (6.1%), OD1 (5.5%), and Bacteroidetes (3.4%), whereas the five most abundant phyla of normal microbiota were Proteobacteria (39.0%), Firmicutes (20.1%), Fusobacteria (4.1%), Actinobacteria (4.1%), and Tenericutes (2.5%). Thus, the phylum Proteobacteria dominated all sample types, but other phyla differed between the fish microbiota and water bacterial communities. Moreover, even if Proteobacteria, Firmicutes, and Actinobacteria were present in similar abundances between dwarf and normal microbiota, the phyla OD1 and Bacteroidetes were more present in dwarf whitefish and the phyla Fusobacteria and Firmicutes were more present in the normal whitefish.
Dwarf vs. normal whitefish microbiota: parallelism or not parallelism?
There was a significant difference between the dwarf and the normal whitefish microbiota at the genus level across all lake populations combined (Table 2). When treating each lake separately, the PERMANOVA tests revealed significant differences between dwarf and normal whitefish in Cliff, East, and Témiscouata lakes whereas no significant differences were found in Indian and Webster lakes (Table 2). Moreover, there is a gradient of genetic population distance between dwarf and normal whitefish from different lakes (Table 1) [56, 74]. Namely, sympatric whitefish from Cliff Lake are the most genetically differentiated (FST = 0.28) whereas those from Témiscouata Lake are the least differentiated (FST = 0.01). Thus, if there was some association between the extent of genetic divergence and the difference in microbiota, dwarf and normal whitefish from Cliff should have the most differentiated intestinal microbiota and Témiscouata should have the least differentiated ones. This was not the case as species specific microbiota was observed in the latter lake, whereas no significant difference was found in both Indian and Webster lakes where genetic differentiation between dwarf and normal whitefish is more pronounced (FST Indian = 0.06 and FST Webster = 0.11).
The weighted UniFrac, which took into account the bacterial abundance rather than simply the presence or absence of taxa in the samples, were significant in all lake populations (Table 2). Therefore, the taxonomic composition of the microbiota was not always different between whitefish species depending on the lake but the abundance of microbiota always differed between whitefish species within each lake. No global differentiation was observed between whitefish species or lakes when all samples were included in the PCoA (Fig. 2a). However, the analysis revealed partially overlapping clusters corresponding to dwarf and normal whitefish in Cliff, East, Témiscouata, and Webster lakes (Fig. 2b, f). Dwarf and normal whitefish clusters were close to each other but nevertheless distinct. For example, in Cliff Lake, the dwarf cluster was more separated by axis one, whereas the normal cluster was more differentiated by axis two. In East, Témiscouata, and Webster Lakes, the opposite pattern was observed: dwarf and normal clusters were better separated by axis two and axis one, respectively. However, only three dwarf whitefish from Webster Lake could be collected resulting in low power of discrimination in that lake. Finally, dwarf and normal whitefish clusters almost completely overlapped in Indian Pond.
Based on the network analysis, the five networks corresponding to each lake gave results that were similar to those obtained with the PCoA analysis, further supporting the observation that the dwarf-normal difference in microbiota varies according to the lake (Fig. 3). Although the network analysis containing all the fish samples revealed no clear pattern, lake-specific networks tended to cluster dwarf and normal samples separately in Cliff and Témiscouata Lakes. Even if the pattern is less clear for East Lake, the dwarf whitefish microbiota from this lake tended to cluster together (but not the normal whitefish microbiota). Also, no clear difference was observed in Indian Pond and as in previous analyses, interpreting patterns observed in Webster Lake was hampered by the small sample size of dwarfs, although microbiota of normal whitefish clustered together.
Functional annotation of whitefish microbiota
Putative microbiota functions were predicted using PICRUSt by assignment of the predicted metagenome (Fig. 4). The gene category, which represented a set of genes influencing the same functional profile, varied widely according to the whitefish species or lake. Only one gene category, cell communication, was stable and had very low gene abundance. Some gene categories, including membrane transport, transcription, or energy metabolism, had high gene abundance in all dwarf and normal whitefish. However, the predicted microbiota functions revealed no significant functional differences between dwarf and normal whitefish microbiota within a given lake except for Témiscouata Lake (Table 2). Globally, there was no significant functional difference between dwarf and normal whitefish microbiota across all lakes combined. Instead, gene abundance differed among lakes and the interaction term between lake populations and species was significant, indicating a strong lake population effect but no significant functional differences between species (Table 2).
Complementary analysis on whitefish microbiota: diversity, core intestinal microbiota, and Metastats
There was no difference between the dwarf and the normal whitefish in terms of bacterial diversity. Thus, the inverse Simpson index was not significant either between species within lakes or between lakes (Table 3). Similar results were also obtained using the Shannon index.
The core intestinal microbiota was defined as the microbial component shared by 80% of the samples. Three genera were shared among all the lake whitefish populations: OD1, Methylobacterium, and Clostridium. Additionally, all dwarf whitefish populations shared Flavobacterium, TM7, and Pseudomonas, whereas all normal whitefish populations shared Aeromonas. Within a given lake, more genera were shared between dwarf and normal whitefish, their number varying between four and 11 depending on the lake (Fig. 1b). Moreover, dwarf whitefish individuals shared more genera than normal whitefish did in Cliff, Indian, Témiscouata, and Webster Lakes. In East Lake, the same number of shared genera was observed between both species. Although the number of shared genera among populations of each species or among lakes was modest, they represented on average 49.5% of all dwarf whitefish shared sequences and 39% of all normal whitefish shared sequences (Table 1).
The Metastats analysis did not allow identifying any genera that were only present in one species. However, several genera were found in only one species within a given lake. These genera were blasted to identify the bacterial taxa being represented (Additional file 3: Table S3). Most of them were bacteria from the environment found in soil, plant, or freshwater. Interestingly, several bacteria previously found in seawater and human clinical specimens (but not found here in the negative control) were also found in intestinal whitefish microbiota, such as Arsenicicoccus piscis, Lactococcus lactis, or Plesiomonas shigelloides [75,76,77]. We also found bacteria known to be pathogenic in fish and humans, such as Flavobacterium spartansii and Clostridium baratii as well as Bifidobacterium thermophilum, which is a probiotic bacterium [78,79,80].
We investigated the intestinal microbiota of sympatric dwarf and normal whitefish pairs in order to (i) test for differences in whitefish intestinal microbiota and water bacterial community from the same lake, (ii) test for differences in intestinal microbiota between dwarf and normal whitefish from the same lake, and (iii) test for the occurrence of parallelism in those patterns. Below, we discuss the main results obtained for each of these objectives, as well as their relevance in the context of ecological speciation.
In order to improve the laboratory protocol and avoid bacterial contamination, meticulous care was taken by working in sterile conditions, performing blank extractions, using positive and negative PCR controls, and sequencing negative PCR controls. These controls revealed very few sequences in negative PCR controls (representing 0.006% of our dataset; Additional file 2: Table S2). These low-contamination sequences were typically associated with fish or fish environments and were represented, in a large majority, by one unique sequence. This contamination is therefore too low to influence the fish mucosa dataset and as such is unlikely to explain the lack of consistent parallelism observed in our dataset. Of the few previous studies that sequenced PCR negative controls, many found contamination without bands following PCR amplification [71,72,73]. Therefore, the PCR negative controls seemed not to be an adequate quality step and in order to know and reduce the risk of contamination, sequencing of PCR negative controls in the case of 16s rRNA gene pyrosequencing should be applied systematically, as we have done here.
Whitefish microbiota vs. water bacterial community within a given lake
The whitefish intestinal microbiota was not reflective of the whitefish environment within each lake tested. Therefore, host physiology, immunity, and genetic background may play a role in determining the internal intestinal microbiota [34, 45, 47, 81]. The taxonomy between the fish intestinal microbiota and the bacterial water community was highly distinct among lakes. The water and the fish bacterial community shared 23, 21, 29, 27, and 23% of genera for Cliff, East, Indian, Témiscouata, and Webster lake populations, respectively. These values are substantially greater than the 5% shared OTUs reported recently between Trinidadian guppies (Poecilia reticulata) and their environment . However, this could be due to the fact that these authors compared the fish microbiota with the bacterial community from both water and sediments. There are two major ways to colonize the fish intestine: via maternal microbial transmission [72, 82] or via the environment, which is the primary mechanism of microbiota acquisition for fish . However, Smith et al. showed that the intestinal microbiota of three-spined stickleback (Gasterosteus aculeatus) tends to be more similar to food-associated bacteria rather than water-associated bacteria . Although we did not sample the whitefish prey, our data demonstrate that around 25% of bacterial genera were shared between water and whitefish microbiota. Moreover, some of the main genera from whitefish microbiota were found at very low frequency in the environment. Therefore, even if the shared bacteria could come from the whitefish diet, it is quite likely that an important proportion of the intestinal microbiota could be attributed to the colonization of bacteria from the water.
Whitefish intestinal vs. kidney microbiota and host effect
In this study, only the bacteria that formed a stable and specific association with the whitefish were analyzed. In fact, only the intestinal adherent microbiota of whitefish was selected, allowing for an indirect investigation of the host effect. In freshwater fishes, the dominant Proteobacteria is reported to be the most abundant phylum . Also, the occurrence Firmicutes, Bacteroidetes, Actinobacteria, Acidobacteria, Chlamydiae, Fusobacteria, Planctomycetes, Spirochaetes, TM7, Verrucomicrobia, and Tenericutes has been reported in many freshwater fishes [38, 41, 42, 84, 85]. However, the phyla OD1, which was present at a relatively low frequency in both dwarf and normal whitefish, has usually been reported in freshwater samples but not freshwater fish, further supporting the acquisition of part of whitefish microbiota from the environment [86, 87].
Globally, we observed a total of 421 different genera in the intestinal mucosa from 108 fish. This is comparable to the level of diversity reported in other recent studies that analyzed 30 intestinal contents of five wild African cichlid fish species (tribe Perissodini) and 72 feces of the wild Amazonian fish tambaqui (Colossoma macropomum) that reported 121 and 525 genera, respectively [47, 88]. Therefore, the number of genera adherent to whitefish intestinal mucosa was similar to the number of genera found in feces or intestinal content in other wild freshwater fish. In a previous study of the kidney bacterial community in lake whitefish , the observed genera diversity (579 genera from 133 apparently healthy fish) was higher than that observed here for the intestinal mucosa. However, many more OTUs (24,308 OTUs) were found in the intestinal mucosa than in the kidneys (2168 OTUs). In both studies, mature fish were sampled in the same environment and they were sampled at the same period of time but in different years. The difference in genera diversity may result from both host genetic and immunity effects. Although the intestinal tract of animals contains the largest number of bacteria, which explains the difference between the intestinal mucosa and the kidney microbiomes at the OTU level, bacterial selection by the host may stabilize the number of intestinal genera [14, 16, 17, 81]. Such host-driven selection was highlighted in a zebrafish (Danio rerio) intestinal microbiota study, where the number of OTUs decreased during zebrafish development until reaching an equilibrium at fish maturity .
Interestingly, our data revealed no difference in diversity between intestinal microbiota of dwarf and normal whitefish found in sympatry within a given lake. This is in contrast with our previous study on kidney tissues where normal whitefish harbored a higher diversity than dwarf whitefish in all five lakes studied . We had proposed that this difference may come from the distinct trophic niche of the two whitefish species. Dwarf whitefish feed almost exclusively on zooplankton [57, 58], whereas normal whitefish are generalists and feed on zoobenthos, molluscs, and fish prey [50, 58]. Moreover, Bolnick et al. observed a less diverse intestinal microbiota when the food was more diversified in both three-spined stickleback and Eurasian perch (Gasterosteus aculeatus and Perca fluviatilis), suggesting that the host had an effect on bacterial diversity . Thus, the strikingly different diets between dwarf and normal whitefish had no apparent effect on the diversity of the adherent intestinal microbiota. As mentioned above, host genetic effects could select commensal bacteria in its intestine, which could perhaps explain the similar diversity level observed between dwarf and normal whitefish. Indeed, while the intestinal microbiota lives in a tight symbiotic relationship with the host, this is less so the case for kidney where the kidney microbiota has more of a pathogenic relationship with the host [16, 49]. Therefore, the comparison between symbiotic and pathogenic relationship could highlight the important host effect on the stabilization of the intestinal microbiota but not in the kidney.
Sequencing the microbial world has revealed an overwhelming intestinal microbiota impact on the host and has allowed documenting the core intestinal microbial communities in mammalian and teleost fish [3, 39, 40, 42, 45, 91,92,93]. The core intestinal microbiota corresponds to the OTUs or the genera shared among close host relatives and could be horizontally transmitted and/or selected as a common set of bacteria [3, 47]. For example, Roeselers et al. documented the occurrence of core intestinal microbiota between the domesticated and wild Zebrafish (Danio rerio) . Here, our core microbiota data represented between 22 and 65% (mean ~ 44%) of genera shared between both species in each lake (Table 1). This percent of shared sequences is higher than that reported by Baldo et al., which found that the intestinal microbiota of cichlid species shared between 13 and 15% of sequences, but was equivalent to Sullam et al., which reported around 50% of shared sequences in the intestinal microbiota of Trinidadian guppy ecotypes [45, 47]. Therefore, the conservation of the core microbiota was strong within each whitefish species for each lake, further supporting the hypothesis of a strong host selective effect on its microbiota.
No clear evidence for parallelism in intestinal microbiota between dwarf and normal whitefish
Parallelism is the evolution of similar traits in independent populations . In the case of lake whitefish, the test for patterns of parallelism at many different levels may help identify the main factors that are at play in driving the process of ecological speciation in this system of repeated sympatric pairs. Here, given the many differences in their ecology and life history traits, we expected to observe some parallelism in differential intestinal microbiota between dwarf and normal whitefish species pairs. Indeed, parallelism between dwarf and normal whitefish has previously been documented for morphological, physiological, behavioral, and ecological traits [53, 55, 95,96,97,98,99,100,101]. Parallelism was also documented at the gene expression level, whereby dwarf whitefish consistently show significant overexpression of genes implicated with survival functions whereas normal whitefish show overexpression of genes associated with growth functions [56, 96]. Therefore, the apparent lack of parallelism in intestinal microbiota is somewhat surprising, especially given the known difference in trophic niches occupied by dwarf and normal whitefish. Indeed, fish diet is known to alter microbiota composition [83, 102,103,104,105]. Moreover, microbiotas have been reported to change in parallel with their host phylogeny [15, 17]. This phenomena coined “phylosymbiosis” has been reported in organisms as phylogenetically diverse as hydra, fish, and primates [40, 106, 107]. Here, we performed seven different types of analyses to test whether there were differences in the intestinal microbiota of the five whitefish species pairs that could have highlighted the occurrence of parallelism. However, while a clear difference between dwarf and normal whitefish microbiota composition was observed in three lakes, these differences were not parallel among lakes. Moreover, there was no difference between dwarf and normal whitefish from the other two lakes. Although the bacterial abundance (weighted UniFrac) differed between species in all five lakes, again, those differences were not parallel across lakes.
All in all, we found no clear evidence of parallelism in the intestinal microbiota across the five dwarf and normal whitefish species pairs. Instead, our results suggested that the main source of variation in whitefish microbiota was the lake of origin. As mentioned above, an important proportion of the intestinal microbiota could be attributed to the colonization by bacteria from the water. However, each lake studied had a distinct water bacterial community (PERMANOVA, water bacterial community of all the lakes = 0.0025). Although the whitefish host could modulate the intestinal microbiota, the lake bacterial variation could positively or negatively influence the intestinal microbiota of whitefish species. Cliff, Webster, and Indian lakes harbor the most genetically divergent species pairs, whereas East and Témiscouata species pairs are the least differentiated [51, 74]. These two groups of lakes are characterized by important environmental differences . More specifically, lakes with the most divergent populations are characterized by the greatest oxygen depletion and lower zooplankton densities, suggesting harsher environmental conditions favoring more pronounced competition for resources between the two species . On the contrary, lakes with the less divergent populations were characterized by more favorable environmental conditions . Among the three lakes with the most genetically divergent species pairs, dwarf and normal whitefish differed in their intestinal microbiota only in Cliff Lake. East and Témiscouata species pairs (the two least differentiated populations) were also characterized by distinct intestinal microbiota. These observations suggest that while the lake of origin explains the composition of whitefish intestinal microbiota better than the species, there is no clear association between lake abiotic and biotic characteristics and the fish microbiota, suggesting that other factors that still need to be elucidated are at play.
Whitefish microbiotas and their possible role in ecological speciation
Most of adherent bacteria living on the intestinal mucosa are not randomly acquired from the environment , but are rather retained by different factors in the host . These symbiotic bacteria may play an essential role in the ecology and evolution of their hosts. Indeed, certain symbionts may affect evolutionary trajectories by conferring fitness advantages [26, 109]. For example, the microbiota of the desert woodrats (Neotoma lepida) enables its host to feed on creosote toxic compounds, suggesting a fitness advantage by limiting resource competition . Symbionts can also influence speciation in several ways. First, there are two main processes that could influence pre-zygotic isolation: (i) microbe-specific, which involves bacterium-derived products such as metabolites and (ii) microbe-assisted, which involves bacterial modulation of the host-derived odorous products [14, 17]. In a recent study, Damodaram et al. showed that the attraction of male to female fruit flies is abolished when female flies are fed with antibiotics, implying a role of the fly’s microbiota in mate choice . Second, symbionts can influence post-zygotic reproductive isolation with, for example, cytoplasmic incompatibilities leading to hybrid inviability . These authors made crosses between two species of Nasonia wasp (Nasonia vitripennis and Nasonia giraulti) to create F2 hybrid larvae raised with their symbionts (conventional rearing) and without the symbionts (germ free). The F2 lethality was clearly more important with symbionts (conventional rearing) than without symbionts (germ free). Moreover, this lethality was not seen in pure larvae of both species reared with symbionts. Symbionts can also increase the host phenotype plasticity . For example, a facultative endo-symbiotic bacterium called pea aphid U-type symbiont (PAUS) allowed the pea aphid (Acyrthosiphon pisum) to acquire a new phenotype: the digestive capability of alfalfa (Medicago sativum) . This new phenotype supports a niche expansion that leads to geographic isolation between aphid populations and therefore indirectly confers a mechanism for pre-zygotic isolation. Given the absence of clear association between whitefish intestinal microbiota and whitefish species, it thus seems unlikely that any of these processes are at play in the speciation of the whitefish species pairs. This absence of parallelism across dwarf vs. normal whitefish microbiota highlights the complexity of the holobiont and suggests that the direction of selection could be different between the host and its microbiota.
In summary, we analyzed the intestinal microbiota in the context of population divergence and speciation in the natural environments. We selected the whitefish mucosa; only the bacteria which formed a stable and specific association with the whitefish were analyzed. To our knowledge, this is the very first study which sequenced the intestinal adherent microbiota in natural fish host populations. Our main goal was to test for the occurrence of parallelism in the microbiota of dwarf and normal whitefish that evolved in parallel across different environments. However, no clear evidence for parallelism was observed at the bacterial level. We found distinct microbiota between the dwarf and the normal species in three of the five lake populations suggesting more selective pressure from the environment. This absence of parallelism across dwarf vs. normal whitefish microbiota highlighted the complexity of the holobiont and suggests that the direction of selection could be different between the host and its microbiota. Furthermore, the comparison of the adherent microbiota with the water bacterial environment and whitefish kidney bacterial community  provided evidence for selection of the adherent bacteria composition made by the host as well as bacterial diversity stabilization. Finally, an experiment without environmental variation would be useful to limit the effect of this in order to determine whether differences between whitefish species remain as large as observed here.
- 16 s rRNA:
16S Ribosomal RNA
Analysis of variance
Comités de protection des animaux
Generalized linear model
Operational taxonomic unit
Principal Coordinates Analysis
Polymerase chain reaction
Permutational analysis of variance
Phylogenetic Investigation of Communities by Reconstruction of Unobserved States
Ribosomal Database Project
Years before present
Whitman WB, Coleman DC. Prokaryotes: the unseen majority. Proc Natl Acad Sci U S A. 1998;95:6578–83.
Huse SM, Dethlefsen L, Huber JA, Welch DM. Exploring microbial diversity and taxonomy using SSU rRNA hypervariable tag sequencing. PLoS Genet. 2008;4:e1000255.
Rawls JF, Mahowald MA, Ley RE, Gordon JI. Reciprocal gut microbiota transplants from zebrafish and mice to germ-free recipients reveal host habitat selection. Cell. 2006;127:423–33.
Brucker RM, Bordenstein SR. The capacious hologenome. Zoology. 2013;116:260–1.
Archie EA, Theis KR. Animal behaviour meets microbial ecology. Anim Behav. 2011;82:425–36.
McFall-Ngai M, Heath-Heckman EAC, Gillette AA, Peyer SM, Harvie EA. The secret languages of coevolved symbioses: insights from the Euprymna scolopes–Vibrio fischeri symbiosis. Semin Immunol. 2012;24:3–8.
Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, Pettersson S. Host-gut microbiota metabolic interactions. Science. 2012;336:1262–7.
Ezenwa VO, Gerardo NM, Inouye DW, Medina M, Xavier JB. Animal behavior and the microbiome. Science. 2012;338:198–9.
McCutcheon JP, von CD D. An interdependent metabolic patchwork in the nested symbiosis of mealybugs. Curr Biol. 2011;21:1366–72.
Gilbert SF, Sapp J, Tauber AI. A symbiotic view of life: we have never been individuals. Q Rev Biol. 2012;87:325.
Wang J, Kalyan S, Steck N, Turner LM, Harr B, Baines JF. Analysis of intestinal microbiota in hybrid house mice reveals evolutionary divergence in a vertebrate hologenome. Nat Commun. 2015;6:6440.
Lee YK, Mazmanian SK. Has the microbiota played a critical role in the evolution of the adaptive immune system? Science. 2010;330:1768–72.
Rosenberg E, ZilberRosenberg I. Symbiosis and development: the hologenome concept. Birth Defects Res C Embryo Today. 2011;93:56–66.
Brucker RM, Bordenstein SR. Speciation by symbiosis. Trends Ecol Evol. 2012;27:443–51.
Bordenstein SR, Theis KR. Host biology in light of the microbiome: ten principles of holobionts and hologenomes. PLoS Biol. 2015;13:e1002226.
Rosenberg E, Zilber I. The hologenome concept: human, animal and plant microbiota. Switzerland: Springer; 2013.
Shropshire JD, Bordenstein SR. Speciation by symbiosis: the microbiome and behavior. MBio. 2016;7:e01785–15.
Carlson AD, Copeland J, Raderman R, Bulloch A. Role of interflash intervals in a firefly courtship (Photinus macdermotti). Anim Behav. 1976;24:786–92.
Cator LJ, Arthur BJ, Harrington LC, Hoy RR. Harmonic convergence in the love songs of the dengue vector mosquito. Science. 2009;323:1077.
Matsuura K. Nestmate recognition mediated by intestinal bacteria in a termite, Reticulitermes speratus. Oikos. 2001;92:20–6.
De Cock R, Matthysen E. Sexual communication by pheromones in a firefly, Phosphaenus hemipterus (Coleoptera: Lampyridae). Anim Behav. 2005;70:807–18.
Damodaram K, Ayyasamy A, Kempraj V. Commensal bacteria aid mate-selection in the fruit fly, Bactrocera dorsalis. Microb Ecol. 2016;72:725–9.
Pfennig DW, Wund MA, Snell EC. Phenotypic plasticity’s impacts on diversification and speciation. Trends Ecol Evol. 2010;25:459–67.
Kohl KD, Stengel A, Dearing MD. Inoculation of tannin-degrading bacteria into novel hosts increases performance on tannin-rich diets. Environ Microbiol. 2015;18:1720–9.
David LA, Maurice CF, Carmody RN, Gootenberg DB, Turnbaugh PJ. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505:559–63.
Kohl KD, Weiss RB, Cox J, Dale C, Dearing MD. Gut microbes of mammalian herbivores facilitate intake of plant toxins. Ecol Lett. 2014;17:1238–46.
Taylor MW, Radax R, Steger D, Wagner M. Sponge-associated microorganisms: evolution, ecology, and biotechnological potential. Microbiol Mol Biol Rev. 2007;71:295–347.
Kamra DN. Rumen microbial ecosystem. Appl Environ Microbiol. 2005;2:303–8.
Broderick NA, Buchon N, Lemaitre B. Microbiota-induced changes in Drosophila melanogaster host gene expression and gut morphology. MBio. 2014;5:e01117–4.
Lewis Z, Lizé A. Insect behaviour and the microbiome. Curr Opin Insect Sci. 2015;9:86–90.
Cryan JF, Dinan TG. Mind-altering microorganisms: the impact of the gut microbiota on brain and behaviour. Nat Rev Neurosci. 2012;13:701–12.
Dantas G, Sommer M, Degnan PH, Goodman AL. Experimental approaches for defining functional roles of microbes in the human gut. Annu Rev Microbiol. 2013;67:459–75.
Hehemann JH, Correc G, Barbeyron T, Helbert W, Czjzek M, Michel G. Transfer of carbohydrate-active enzymes from marine bacteria to Japanese gut microbiota. Nature. 2010;464:908–12.
Alberdi A, Aizpurua O, Bohmann K, Zepeda-Mendoza ML, Gilbert MTP. Do vertebrate gut metagenomes confer rapid ecological adaptation? Trends Ecol Evol. 2016;31(9):689–99.
Bromham L. Why do species vary in their rate of molecular evolution? Biol Lett. 2009;5:401.
Gogarten JP, Townsend JP. Horizontal gene transfer, genome innovation and evolution. Nat Rev Microbiol. 2005;3:679–87.
Nelson JS. Fishes of the world. 4th ed. New York: John; 2006.
Sullam KE, Essinger SD, Lozupone CA, O’Connor MP, Rosen GL, Knight R, Kilham SS, Russell GA. Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. Mol Ecol. 2012;21(3):3363–78.
Llewellyn MS, McGinnity P, Dionne M, Letourneau J, Thonier F, Carvalho GR, Creer S, Derome N. The biogeography of the Atlantic salmon (Salmo salar) gut microbiome. ISME J. 2015;10:1280–4.
Miyake S, Ngugi DK, Stingl U. Diet strongly influences the gut microbiota of surgeonfishes. Mol Ecol. 2015;24:656–72.
Ye L, Amberg J, Chapman D, Gaikowski M, Liu WT. Fish gut microbiota analysis differentiates physiology and behavior of invasive Asian carp and indigenous American fish. ISME J. 2014;10:2076.
Roeselers G, Mittge EK, Stephens WZ, Parichy DM, Cavanaugh CM, Guillemin K, Rawls JF. Evidence for a core gut microbiota in the zebrafish. ISME J. 2011;5:1595–608.
Shiina A, Itoi S, Washio S, Sugita H. Molecular identification of intestinal microflora in Takifugu niphobles. Comp Biochem Physiol Part D Genomics Proteomics. 2006;1:128–32.
Baldo L, Pretus JL, Riera JL, Musilova Z, Nyom ARB, Salzburger W. Convergence of gut microbiotas in the adaptive radiations of African cichlid fishes. ISME J. 2017; https://doi.org/10.1038/ismej.2017.62.
Sullam KE, Rubin B, Dalton CM, Kilham SS, Flecker AS, Russell JA. Divergence across diet, time and populations rules out parallel evolution in the gut microbiomes of Trinidadian guppies. ISME J. 2015;9:1508–22.
Hata H, Tanabe AS, Yamamoto S, Toju H, Kohda M, Hori M. Diet disparity among sympatric herbivorous cichlids in the same ecomorphs in Lake Tanganyika: amplicon pyrosequences on algal farms and stomach contents. BMC Biol. 2014;12:90.
Baldo L, Riera JL, Tooming-Klunderud A, Albà MM, Salzburger W. Gut microbiota dynamics during dietary shift in Eastern African cichlid fishes. PLoS One. 2015;10:e0127462.
Smith C, Snowberg LK, Caporaso JG, Knight R, Bolnick DI. Dietary input of microbes and host genetic variation shape among-population differences in stickleback gut microbiota. ISME J. 2015;9:2515.
Sevellec M, Pavey SA, Boutin S, Filteau M, Derome N, Bernatchez L. Microbiome investigation in the ecological speciation context of lake whitefish (Coregonus clupeaformis) using next-generation sequencing. J Evol Biol. 2014;27:1029.
Bernatchez L. Ecological theory of adaptive radiation: an empirical assessment from coregonine fishes (Salmoniformes). In: Evolution illuminated: salmon and their relatives. New York: Oxford University Press; 2004. p. 175–207.
Gagnaire PA, Pavey SA, Normandeau E, Bernatchez L. The genetic architecture of reproductive isolation during speciation-with-gene-flow in lake whitefish species pairs assessed by rad sequencing. Evolution. 2013;67:2483.
Pavey SA, Sevellec M, Adam W, Normandeau E, Lamaze FC, Gagnaire PA, Filteau M, Herbert FO, Maaroufi H, Bernatchez L. Nonparallelism in MHCIIβ diversity accompanies nonparallelism in pathogen infection of lake whitefish (Coregonus clupeaformis) species pairs as revealed by next-generation sequencing. Mol Ecol. 2013;22:3833–49.
Laporte M, Dalziel AC, Martin N, Bernatchez L. Adaptation and acclimation of traits associated with swimming capacity in lake whitefish (coregonus clupeaformis) ecotypes. BMC Evol Biol. 2016;16:160.
Jeukens J, Renaut S, StCyr J, Nolte AW, Bernatchez L. The transcriptomics of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis spp., Salmonidae) divergence as revealed by next-generation sequencing. Mol Ecol. 2010;19:5389–403.
Dalziel AC, Martin N, Laporte M, Guderley H, Bernatchez L. Adaptation and acclimation of aerobic exercise physiology in lake whitefish ecotypes (Coregonus clupeaformis). Evolution. 2015;69:2167–86.
Bernatchez L, Renaut S, Whiteley AR, Derome N, Jeukens J, Landry L, Lu G, Nolte AW, Østbye K, Rogers SM, St-Cyr J. On the origin of species: insights from the ecological genomics of lake whitefish. Philos Trans R Soc Lond Ser B Biol Sci. 2010;365:1783–800.
Bodaly RA. Morphological and ecological divergence within the lake whitefish (Coregonus clupeaformis) species complex in Yukon territory. J Fish Board Can. 1979;36:1214.
Bernatchez L, Chouinard A, Lu G. Integrating molecular genetics and ecology in studies of adaptive radiation: whitefish, Coregonus sp., as a case study. Biol J Linnean Soc. 1999;68:173.
Ringø E, Sperstad S, Myklebust R, Refstie S. Characterisation of the microbiota associated with intestine of Atlantic cod (Gadus morhua L.): the effect of fish meal, standard soybean meal. Aquaculture. 2006;255(1–4):279.
Don RH, Cox PT, Wainwright BJ, Baker K. “Touchdown” PCR to circumvent spurious priming during gene amplification. Nucleic Acids Res. 1991;19(14):4008.
Korbie DJ, Mattick JS. Touchdown PCR for increased specificity and sensitivity in PCR amplification. Nat Protoc. 2008;3(9):1452.
Klindworth A, Pruesse E, Schweer T. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2012;41:e1.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112.
Schloss PD, Westcott SL, Ryabin T. Introducing Mothur. Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228.
Oksanen J, Kindt R, Legendre P, O'hara B, Stevens HH. Vegan: community ecology package. R package ver. 2.0–10 2013.
White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009;5:e1000352.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498.
Langille M, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Thurber RLV, Knight R, Beiko RG, Huttenhower C. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814.
Baldwin EA, Walther M, MacLean AM, Gohl DM, Beckman KB, Chen J, White B, Creedon DJ, Chia N. Persistent microbial dysbiosis in preterm premature rupture of membranes from onset until delivery. PeerJ. 2015;3:e1398.
Wilkins L, Rogivue A, Fumagalli L, Wedekind C. Declining diversity of egg-associated bacteria during development of naturally spawned whitefish embryos (Coregonus spp.). Aquat Sci. 2015;77:481–97.
Gu Y, Mao X, Zha L, Hou Y, Yun L. Development of a candidate method for forensic microbial genotyping using multiplex pyrosequencing combined with a universal biotinylated primer. Forensic Sci Int. 2015;246:e1.
Renaut S, Nolte AW, Rogers SM, Derome N, Bernatchez L. SNP signatures of selection on standing genetic variation and their association with adaptive phenotypes along gradients of ecological speciation in lake whitefish species pairs (Coregonus spp.). Mol Ecol. 2011;20:545–59.
Aldova E, Melter O, Chýle P, Slosarek M, Kodym P. Plesiomonas shigelloides in water and fish. Cent Eur J Public Health. 1999;7:172–5.
Itoi S, Abe T, Washio S, Ikuno E, Kanomata Y, Sugita H. Isolation of halotolerant Lactococcus lactis subsp. lactis from intestinal tract of coastal fish. Int J Food Microbiol. 2008;121:116–21.
Hamada M, Iino T, Iwami T, Tamura T, Harayama S, Suzuki K. Arsenicicoccus piscis sp. nov., a mesophilic actinobacterium isolated from the intestinal tract of a fish. Actinomycetologica. 2009;23:40.
Loch TP, Faisal M. Flavobacterium spartansii sp. nov., a pathogen of fishes, and emended descriptions of Flavobacterium aquidurense and Flavobacterium araucananum. Int J Syst Evol Microbiol. 2014;64:406–12.
McCroskey LM, Hatheway CL, Woodruff BA, Greenberg JA, Jurgenson P. Type F botulism due to neurotoxigenic Clostridium baratii from an unknown source in an adult. J Clin Microbiol. 1991;29:2618–20.
Domig KJ, Mayrhofer S, Zitz U, Mair C, Petersson A, Amtmann E, Mayer HK, Kneifel W. Antibiotic susceptibility testing of Bifidobacterium thermophilum and Bifidobacterium pseudolongum strains: broth microdilution vs. agar disc diffusion assay. Int J Food Microbiol. 2007;120:191.
Macke E, Tasiemski A, Massol F, Callens M, Decaestecker E. Life history and eco-evolutionary dynamics in light of the gut microbiota. Oikos. 2017;126:508–31.
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.
Nayak SK. Role of gastrointestinal microbiota in fish. Aquac Res. 2010;41:1553–73.
Li X, Yan Q, Xie S, Hu W, Yu Y, Hu Z. Gut microbiota contributes to the growth of fast-growing transgenic common carp (Cyprinus carpio L.). PLoS One. 2013;8:e64577.
Eichmiller JJ, Hamilton MJ, Staley C, Sadowsky MJ, Sorensen PW. Environment shapes the fecal microbiome of invasive carp species. Microbiome. 2016;4:44.
Smith KF, Schmidt V, Rosen GE, Amaral L. Microbial diversity and potential pathogens in ornamental fish aquarium water. PLoS One. 2012;7:e39971.
Qiu T, Liu L, Gao M, Zhang L, Tursun H, Wang X. Effects of solid-phase denitrification on the nitrate removal and bacterial community structure in recirculating aquaculture system. Biodegradation. 2016;27:165.
Sylvain FÉ, Cheaib B, Llewellyn M, Correia TG, Fagundes DB, Val AL, Derome N. pH drop impacts differentially skin and gut microbiota of the Amazonian fish tambaqui (Colossoma macropomum). Sci Rep. 2016;6:32032.
Stephens WZ, Burns AR, Stagaman K, Wong S, Rawls JF, Guillemin K, Bohannan BJ. The composition of the zebrafish intestinal microbial community varies across development. ISME J. 2016;10:644.
Bolnick DI, Snowberg LK, Hirsch PE, Lauber CL, Knight R, Caporaso JG, Svanbäck R. Individuals’ diet diversity influences gut microbial diversity in two freshwater fish (threespine stickleback and Eurasian perch). Ecol Lett. 2014;17:979.
Tap J, Mondot S, Levenez F, Pelletier E, Caron C, Furet JP, Ugarte E, Muñoz-Tamayo R, Paslier DL, Nalin R, Dore J, Leclerc M. Towards the human intestinal microbiota phylogenetic core. Environ Microbiol. 2009;11:2574.
Wong S, Waldrop T, Summerfelt S, Davidson J, Barrows F, Kenney PB, Welch T, Wiens GD, Snekvik K, Rawls JF, Good C. Aquacultured rainbow trout (Oncorhynchus mykiss) possess a large core intestinal microbiota that is resistant to variation in diet and rearing density. Appl Environ Microbiol. 2013;79:4974.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.
Schluter D, Nagel LM. Parallel speciation by natural selection. Am Nat. 1995;146:292–301.
Laporte M, Rogers SM, Dion AM, Normandeau E, Gagnaire PA, Dalziel AC, Chebib J, Bernatchez L. RAD-QTL mapping reveals both genome-level parallelism and different genetic architecture underlying the evolution of body shape in lake whitefish (Coregonus clupeaformis) species pairs. G3-Genes Genom Genet. 2015;5:1481.
StCyr J, Derome N, Bernatchez L. The transcriptomics of life-history trade-offs in whitefish species pairs (Coregonus sp.). Mol Ecol. 2008;17:1850–70.
Landry L, Bernatchez L. Role of epibenthic resource opportunities in the parallel evolution of lake whitefish species pairs (Coregonus sp.). J Evol Biol. 2010;23:2602–13.
Dion-Cote AM, Renaut S, Normandeau E, Bernatchez L. RNA-seq reveals transcriptomic shock involving transposable elements reactivation in hybrids of young lake whitefish species. Mol Biol Evol. 2014;31:1188.
Lu G, Bernatchez L. Experimental evidence for reduced hybrid viability between dwarf and normal ecotypes of lake whitefish (Coregonus clupeaformis Mitchill). Philos Trans R Soc Lond Ser B Biol Sci. 1998;265:1400.
Dalziel AC, Laporte M, Rougeux C, Guderley H, Bernatchez L. Convergence in organ size but not energy metabolism enzyme activities among wild lake whitefish (Coregonus clupeaformis) species pairs. Mol Ecol. 2016;26:225–44.
Jeukens J, Bittner D, Knudsen R, Bernatchez L. Candidate genes and adaptive radiation: insights from transcriptional adaptation to the limnetic niche among coregonine fishes (Coregonus spp., Salmonidae). Mol Biol Evol. 2009;26:155–66.
Haygood AM, Jha R. Strategies to modulate the intestinal microbiota of tilapia (Oreochromis sp.) in aquaculture: a review. Rev Aquacult. 2016;1:1–14.
Zarkasi KZ, Taylor RS, Abell G, Tamplin ML, Glencross BD, Bowman JP. Atlantic salmon (Salmo salar L.) gastrointestinal microbial community dynamics in relation to digesta properties and diet. Mol Ecol. 2016;71:589.
David CJ, Veena R, Kumaresan G. High cholesterol diet induces obesity in zebrafish. PLoS One. 2016;8:e66970.
Koo H, Hakim JA, Powell ML, Kumar R, Eipers PG, Morrow CD, Crowley M, Lefkowitz EJ, Watts SA, Bej K. Metagenomics approach to the study of the gut microbiome structure and function in zebrafish Danio rerio fed with gluten formulated diet. J Microbiol Methods. 2017;135:69–76.
Ochman H, Worobey M, Kuo CH, Ndjango J, Peeters M, Hahn BH, Hugenholtz P. Evolutionary relationships of wild hominids recapitulated by gut microbial communities. PLoS Biol. 2010;8:e1000546.
Franzenburg S, Walter J, Künzel S, Wang J, Baines JF, Bosch TCG, Fraune S. Distinct antimicrobial peptide expression determines host species-specific bacterial associations. Proc Natl Acad Sci U S A. 2013;110:e3730.
Landry L, Vincent WF, Bernatchez L. Parallel evolution of lake whitefish dwarf ecotypes in association with limnological features of their adaptive landscape. J Evol Biol. 2007;20:971.
Tsuchida T, Koga R, Fukatsu T. Host plant specialization governed by facultative symbiont. Science. 2004;303:1989.
We thank G. Côté and the members of the Bernatchez laboratory for fieldwork, C. Hernandez-Chàvez and Alex Bernatchez for their laboratory advice and support, B. Boyle for his help with the Illumina MiSeq sequencing, and B.J.G. Sutherland and A. Xuereb for comments on an earlier version of this manuscript.
Funding for this project was provided by the National Sciences and Engineering Council of Canada (NSERC), a Canadian Research Chair in genomics and conservation of aquatic resources to LB. It is also a contribution to the research program of Québec-Océan.
Availability of data and materials
Sequencing results are available in the Sequence Read Archive (SRA) database at NCBI under the BioProject ID PRJNA394764.
Ethics approval and consent to participate
All animals used in this study were treated in accordance with the approved University of Laval’s CPA (comités de protection des animaux).
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.
Table S1. Steps used to reduce sequencing and PCR errors. (XLSX 9 kb)
Table S2. Bacterial taxa found in the PCR negative control. (XLSX 19 kb)
Table S3. Bacterial species specific to a single whitefish species within lake obtained with Metastats and BLAST. (XLSX 17 kb)