Research article | Open | Published:
Next generation sequencing shows high variation of the intestinal microbial species composition in Atlantic cod caught at a single location
BMC Microbiologyvolume 13, Article number: 248 (2013)
The observation that specific members of the microbial intestinal community can be shared among vertebrate hosts has promoted the concept of a core microbiota whose composition is determined by host-specific selection. Most studies investigating this concept in individual hosts have focused on mammals, yet the diversity of fish lineages provides unique comparative opportunities from an evolutionary, immunological and environmental perspective. Here we describe microbial intestinal communities of eleven individual Atlantic cod (Gadus morhua) caught at a single location based on an extensively 454 sequenced 16S rRNA library of the V3 region.
We obtained a total of 280447 sequences and identify 573 Operational Taxonomic Units (OTUs) at 97% sequence similarity level, ranging from 40 to 228 OTUs per individual. We find that ten OTUs are shared, though the number of reads of these OTUs is highly variable. This variation is further illustrated by community diversity estimates that fluctuate several orders of magnitude among specimens. The shared OTUs belong to the orders of Vibrionales, which quantitatively dominate the Atlantic cod intestinal microbiota, followed by variable numbers of Bacteroidales, Erysipelotrichales, Clostridiales, Alteromonadales and Deferribacterales.
The microbial intestinal community composition varies significantly in individual Atlantic cod specimens caught at a single location. This high variation among specimens suggests that a complex combination of factors influence the species distribution of these intestinal communities.
The intestinal microbial community provides a variety of crucial functions for their vertebrate hosts e.g. , though the factors that influence the colonization of this habitat are less understood. Common patterns among microbial communities of different hosts have promoted the concept of a core set of species, which provides a minimal functionality in the healthy gut and which is determined by host-specific selection [2, 3]. For example, host transcriptional responses to microbial colonization appear to be conserved among a wide range of vertebrates, including fish . Moreover, within the intestinal community of humans, some species are more prevalent [3, 5, 6] and functional gene profiles are highly similar among individuals . Nevertheless, the utility of the core microbiota concept at a fine taxonomic level has recently been questioned due to limited evidence of universally abundant species in humans [8, 9].
Fish provide unique opportunities to investigate the factors that influence the composition of the vertebrate intestinal microbiota due to their high species diversity , dietary variation or habitat preferences , and divergent immune architecture. For instance, considering the differences in immune systems as an example, Atlantic cod lacks the antigen presenting major histocompatibility complex (MHC) II system, which was thought to be conserved among all jawed vertebrates . This lack of MHC II may affect the interactions of Atlantic cod with its microbial community . A extensive meta-analysis -based on uncultured and cultured sampling methods- indicates that the composition of the intestinal communities in teleosts is influenced by both abiotic and biotic factors . Nevertheless, this meta-analysis is predominantly based on pooled Sanger sequencing data, and studies investigating microbial communities in fish using high-throughput sequencing are relatively rare. Moreover, the studies that employed these methods so far have focused on fresh water species held in semi-controlled environments [14–16]. One exception investigating natural populations of zebrafish, identified a core intestinal microbiota based on shared Operational Taxonomic Units (OTUs), despite substantial differences in host provenance and domestication status . This study pooled 4, 6 and 20 individuals respectively, before sequencing . Therefore, to our knowledge, a characterization of the microbial community using high-through methodologies in wild-caught, individual fish is still lacking. Here we investigate the intestinal microbial communities of 11 wild-caught Atlantic cod collected at a single location and quantify a core microbiota based on shared membership in a 454 sequenced 16S rRNA V3 region amplicon dataset.
Results and discussion
We obtained 280447 sequences of approximately 200 basepair (bp) of the 16S rRNA V3 region and identified 573 OTUs at 97% sequence similarity. Rarefaction curve analysis, depicting the relationship between the number of detected OTUs and read number, shows that the final number of OTUs per sample (ranging from 40 to 228) is not caused by uneven sequencing depth (Figure 1a). This variation among samples also appears in estimates of community diversity (based on Shannon and Inverse Simpson indices), which vary an order of magnitude (Table 1). Analogous to the intestinal community composition in zebrafish  or human e.g. , the samples are typically dominated by a few abundant OTUs, while the majority of OTUs is present at rare frequency (e.g., 62% of the 573 OTUs occur once). At 97% sequence similarity, 10 OTUs are shared that are highly abundant based on the number of reads (Figure 1b). The number of reads assigned to these OTUs varies substantially among individuals, and no more than five OTUs are shared using a detection cut-off value of at least five reads (reflecting a 99% detection probability assuming a binominal distribution, Additional file 1: Table S1). Moreover, for sequence similarity values above 80% the number of shared OTUs is fairly constant, indicating that this number is not a result of restrictive cut-off values when clustering (Figure 1c). Overall, the shared OTUs represent a fraction of the overall sequence diversity for a wide range of cut-off values (Figure 1c).
At order level classification, the intestinal community of Atlantic cod is dominated by Vibrionales followed by variable numbers of Bacteroidales, Erysipelotrichales and Clostridiales (Figure 2a). The high proportion Vibrionales (50%) agrees with those proportions found in a meta-analysis based on GenBank sequences of other marine carnivores  and bacteria from this order have also previously been isolated from the Atlantic cod gut e.g. . Nevertheless, the intestinal community also contains a substantial proportion of Bacteriodales (17%). Such abundance has previously been proposed to be a characteristic of the microbial community of marine herbivores, and this finding suggests that the distinction between herbivorous and carnivorous fish may be more subtle . Members of the most abundant orders agree with those reported previously in Atlantic cod using both culture-dependent and culture-independent techniques [18–22]. In addition, using high throughput sequencing, several more orders are detected that are relatively rare. Shared OTUs belong to the orders Vibrionales, Bacteroidales, Erysipelotrichales, Clostridiales, Alteromonadales and Deferribacterales (Figure 2b). Overall, taxonomical diversity (based on number of OTUs per order) does not necessarily correlate to the number of reads per order.
To our knowledge, our dataset provides the first characterization using high throughput sequencing of individual intestinal microbial community structure in a natural population of marine fish. It is possible that our sampling retrieved fish from different populations. Nevertheless, tagging studies in the Norwegian Skagerrak coastal region have shown that adult Atlantic cod have confined home ranges  and genetic studies have revealed fine scale geographical population structure [24, 25]. Considering our single sample location, far from the fjord exit, and comparable size of individuals (Additional file 1: Table S2), we assume that our individuals were retrieved from a local population experiencing similar environmental conditions.
Among our samples, we find 10 shared OTUs, with profound variation in the number of reads per individual. This number of shared OTUs may be an overestimation as the fish were kept in a single tank after capture and thus experienced the same environmental conditions before sampling. We presume such a similar environment is more likely to homogenize microbial communities, rather than promote individual differences. Nevertheless, this shared number of OTUs appears relatively low compared to the number of shared OTUs (21 OTUs, at 97% sequence identity cut-off) among populations of zebrafish from radically different environmental conditions, coming either from natural populations in India or from artificial environments in two separate laboratories in the USA . For now, this difference in shared OTUs between our study and the study focusing on zebrafish is difficult to interpret due to methodological variation e.g. pooled versus individual samples, V1-V2 versus V3 16S rRNA region, . It will be interesting to investigate if these differences in shared OTUs membership are environmentally determined (e.g., a largely different food preference and habitat) or are species specific (e.g., the unusual Atlantic cod immune system which might affect its host-microbe interactions [12, 13]).
Community diversity estimates based on 454 amplicon data are influenced by methodological factors such as fragment length, PCR bias and choice of 16S rRNA gene region. Specifically, shorter amplicon lengths (e.g. < 400 bp) may result in relatively higher diversity estimates compared to longer fragments  and arguably provide a better assessment of community structure . In contrast, species richness estimate based on analyses of the 16 s rRNA V3 region appears to slightly underestimate diversity relative to the full-length gene . Such methodological issues make it difficult to compare community diversity across different studies , although metrics that use both richness and relative abundance (i.e. Shannon and Inverse Simpson indices) appear robust , in particular considering our extensive sequencing depth . Interestingly, these metrics fluctuate several orders of magnitude among our different specimens, and show large individual variation in community composition and diversity. The most diverse individuals appear to have a comparable community complexity relative to those found in humans [7, 32].
A variety of properties, such as shared OTU membership, shared phylogeny, persistence or connectivity can be used to define microbial cores . Here we investigated a core microbiota based on shared membership. Definitions for such a core have been proposed ranging from a lineage present in more than half the population  to an abundant lineage shared among all individuals . We argue that the utility of such concept depends on the specificity with which it describes a biological phenomenon and favor the idea that a lineage should be reliably identified among all individuals in order to belong to a core microbiota, hence with a detection probability of at least 99%. According to such definition, the presence of a core microbiota at 97% sequence similarity consists of five OTUs in this population of Atlantic cod. These OTUs belong to orders Vibrionales, Bacteroidales, Erysipelotrichales, Clostridiales and Alteromonadales.
It is possible that the observation of a shared OTU membership can be explained by other factors other than host-specific selection. For example, between teleost fish, the colonization and community structure of the microbial gut community appears better explained by environmental factors such as food choice or habitat (i.e. salinity) than by host phylogeny [11, 34]. Considering our single sample location, it is currently unclear if the observed core microbiota in Atlantic cod is explained by host-specific selection or driven by shared environmental factors. Interestingly, human microbial gut communities are functionally remarkably similar, despite extensive variation in taxonomic composition [7–9]. This functional redundancy may provide support for a ‘founder takes all’ process of colonization, in which a successful colonizer can prevent the subsequent colonization by other, functionally similar strains through high density blocking . Such a stochastic process could lead to the high variation in community composition that we observe among our different specimens.
Based on the extensive 454 sequencing of a 16S rRNA V3 region amplicon library, we find that the OTU based community diversity estimates of the intestinal microbial community in wild-caught Atlantic cod vary significantly among individuals collected at a single location. This individual level variation suggests that a complex combination of factors influences the microbial species distribution in these intestinal communities. Importantly, such variation has gone unobserved in previous studies of natural populations of teleosts whereby samples of pooled individuals were analyzed [11, 17], which may affect estimates of the number of shared OTUs among hosts.
Live Atlantic cod were collected at a single location (N59.871278, W10.587208) using a fish trap in the Oslo fjord, Norway (Additional file 1) and transported to an animal facility approved by the Norwegian Animal Research Authority (NARA, http://oslovet.norecopa.no/dokument.aspx?dokument=67, approval number 155/2008). The specimens were kept in a common tank (2000 l), at ambient water temperature and light conditions (i.e., 6°C and L:D 8:16, respectively) without feed for between seven and twelve days before sampling to help reduce variation in community composition due to the presence of food items . The fish were humanely sacrificed by a blow to the head (without any administration of other sedatives) before sampling. The experiments were approved by NARA’s authorized representative at the facility and were 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.
Flushed contents from the intestinal tract, including the rectal portion and excluding the stomach, were used for DNA isolation (Additional file 1) and PCR amplification after which individually barcoded amplicons were pooled and sequenced using 454 technology . Sequence data was binned into the individual samples using the barcoded tags and subsequently cleaned from artifacts. Diversity estimates were calculated using Mothur . Representative OTU sequences were compared to the SILVA SSU ref NR V108 database  (http://www.arb-silva.de) using BlastN  and classified using the LCA algorithm in Megan .
Ivanov II, Littman DR: Modulation of immune homeostasis by commensal bacteria. Curr Opin Microbiol. 2011, 14 (1): 106-114. 10.1016/j.mib.2010.12.003.
Maynard CL, Elson CO, Hatton RD, Weaver CT: Reciprocal interactions of the intestinal microbiota and immune system. Nature. 2012, 489 (7415): 231-241. 10.1038/nature11551.
Tap J, Mondot S, Levenez F, Pelletier E, Caron C, Furet J-P, Ugarte E, Muñoz-Tamayo R, Paslier DLE, Nalin R, et al: Towards the human intestinal microbiota phylogenetic core. Environ Microbiol. 2009, 11 (10): 2574-2584. 10.1111/j.1462-2920.2009.01982.x.
Rawls JF, Samuel BS, Gordon JI: Gnotobiotic zebrafish reveal evolutionarily conserved responses to the gut microbiota. Proc Natl Acad Sci USA. 2004, 101 (13): 4596-4601. 10.1073/pnas.0400706101.
Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, et al: A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010, 464 (7285): 59-65. 10.1038/nature08821.
Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto J-M, et al: Enterotypes of the human gut microbiome. Nature. 2011, 473 (7346): 174-180. 10.1038/nature09944.
The Human Microbiome Consortium: Structure, function and diversity of the healthy human microbiome. Nature. 2012, 486 (7402): 207-214. 10.1038/nature11234.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, et al: A core gut microbiome in obese and lean twins. Nature. 2009, 457 (7228): 480-484. 10.1038/nature07540.
Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R: Diversity, stability and resilience of the human gut microbiota. Nature. 2012, 489 (7415): 220-230. 10.1038/nature11550.
Nelson JS: Fishes of the world. 2006, Hoboken, New Yersey: John Wiley & Sons
Sullam KE, Essinger SD, Lozupone CA, O’Connor MP, Rosen GL, Knight ROB, Kilham SS, Russell JA: Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. Mol Ecol. 2012, 21 (13): 3363-3378. 10.1111/j.1365-294X.2012.05552.x.
Star B, Nederbragt AJ, Jentoft S, Grimholt U, Malmstrom M, Gregers TF, Rounge TB, Paulsen J, Solbakken MH, Sharma A, et al: The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011, 477 (7363): 207-210. 10.1038/nature10342.
Star B, Jentoft S: Why does the immune system of Atlantic cod lack MHC II?. Bioessays. 2012, 34 (8): 648-651. 10.1002/bies.201200005.
Geraylou Z, Souffreau C, Rurangwa E, D'Hondt S, Callewaert L, Courtin CM, Delcour JA, Buyse J, Ollevier F: Effects of arabinoxylan-oligosaccharides (AXOS) on juvenile Siberian sturgeon (Acipenser baerii) performance, immune responses and gastrointestinal microbial community. Fish Shellfish Immun. 2012, 33 (4): 718-724. 10.1016/j.fsi.2012.06.010.
Wu S, Wang G, Angert ER, Wang W, Li W, Zou H: Composition, diversity, and origin of the bacterial community in Grass carp intestine. Plos One. 2012, 7 (2): e30440-10.1371/journal.pone.0030440.
van Kessel M, Dutilh B, Neveling K, Kwint M, Veltman J, Flik G, Jetten M, Klaren P, Op den Camp H: Pyrosequencing of 16S rRNA gene amplicons to study the microbiota in the gastrointestinal tract of carp (Cyprinus carpio L.). AMB Express. 2011, 1 (1): 41-10.1186/2191-0855-1-41.
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 (10): 1595-1608. 10.1038/ismej.2011.38.
Ringø E, Sperstad S, Myklebust R, Refstie S, Krogdahl A: 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-841. 10.1016/j.aquaculture.2006.06.030.
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.
Brunvold L, Sandaa R-A, Mikkelsen H, Welde E, Bleie H, Bergh O: 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-327.
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.
Dhanasiri AKS, Brunvold L, Brinchmann MF, Korsnes K, Bergh O, Kiron V: Changes in the intestinal microbiota of wild Atlantic cod Gadus morhua L. upon captive rearing. Microb Ecol. 2011, 61 (1): 20-30. 10.1007/s00248-010-9673-y.
Espeland SH, Gundersen AF, Olsen EM, Knutsen H, Gjøsæter J, Stenseth NC: Home range and elevated egg densities within an inshore spawning ground of coastal cod. ICES J Mar Sci. 2007, 64 (5): 920-928. 10.1093/icesjms/fsm028.
Knutsen H, Jorde PE, Andre C, Stenseth NC: Fine-scaled geographical population structuring in a highly mobile marine species: the Atlantic cod. Mol Ecol. 2003, 12 (2): 385-394. 10.1046/j.1365-294X.2003.01750.x.
Olsen EM, Knutsen H, Gjosaeter J, Jorde PE, Knutsen JA, Stenseth NC: Small-scale biocomplexity in coastal Atlantic cod supporting a Darwinian perspective on fisheries management. Evol Appl. 2008, 1 (3): 524-533. 10.1111/j.1752-4571.2008.00024.x.
Engelbrektson A, Kunin V, Wrighton KC, Zvenigorodsky N, Chen F, Ochman H, Hugenholtz P: Experimental factors affecting PCR-based estimates of microbial species richness and evenness. ISME J. 2010, 4 (5): 642-647. 10.1038/ismej.2009.153.
Huber JA, Morrison HG, Huse SM, Neal PR, Sogin ML, Mark Welch DB: Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure. Environ Microbiol. 2009, 11 (5): 1292-1302. 10.1111/j.1462-2920.2008.01857.x.
Youssef N, Sheik CS, Krumholz LR, Najar FZ, Roe BA, Elshahed MS: Comparison of species richness estimates obtained using nearly complete fragments and simulated pyrosequencing-generated fragments in 16S rRNA gene-based environmental surveys. Appl Environ Microbiol. 2009, 75 (16): 5227-5236. 10.1128/AEM.00592-09.
Schloss PD: The effects of alignment quality, distance calculation method, sequence filtering, and region on the analysis of 16S rRNA gene-based studies. PLoS computational biology. 2010, 6 (7): e1000844-10.1371/journal.pcbi.1000844.
Pinto AJ, Raskin L: PCR biases distort bacterial and archaeal community structure in pyrosequencing datasets. Plos One. 2012, 7 (8): e43093-10.1371/journal.pone.0043093.
Lundin D, Severin I, Logue JB, Östman Ö, Andersson AF, Lindström ES: Which sequencing depth is sufficient to describe patterns in bacterial α- and β-diversity?. Environ Microbiol Rep. 2012, 4 (3): 367-372. 10.1111/j.1758-2229.2012.00345.x.
Dethlefsen L, Huse S, Sogin ML, Relman DA: The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. Plos Biology. 2008, 6 (11): e280-10.1371/journal.pbio.0060280.
Shade A, Handelsman J: Beyond the Venn diagram: the hunt for a core microbiome. Environ Microbiol. 2012, 14 (1): 4-12. 10.1111/j.1462-2920.2011.02585.x.
Nayak SK: Role of gastrointestinal microbiota in fish. Aquac Res. 2010, 41 (11): 1553-1573. 10.1111/j.1365-2109.2010.02546.x.
Waters JM, Fraser CI, Hewitt GM: Founder takes all: density-dependent processes structure biodiversity. TREE. 2013, 28 (2): 78-85.
Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen ZT, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al: Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75 (23): 7537-7541. 10.1128/AEM.01541-09.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35 (21): 7188-7196. 10.1093/nar/gkm864.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Huson D, Richter D, Mitra S, Auch A, Schuster S: Methods for comparative metagenomics. BMC Bioinformatics. 2009, 10 (Suppl 1): S12-10.1186/1471-2105-10-S1-S12.
We thank Jonathan E. Colman for the collection of specimens and the Norwegian Sequencing Centre (http://www.sequencing.uio.no) for 454 amplicon sequencing. We are also grateful for the efforts of two anonymous reviewers, whose insightful comments and suggestions helped improve this manuscript. This work was initiated during the Research Council of Norway (RCN) project “Genome sequencing of cod by exclusive use of ultra-high-throughput sequencing technology (#187940)” and supported by RCN project “Functional and comparative immunology of a teleost's world without MHC II (#222378)”. The raw amplicon sequences are deposited at the SRA archive with run number SRR579544.
The author declare that they have no competing interests.
BS and SJ conceived and designed the experiments. TH, BS and SJ performed the experiments. TH extracted DNA and created the amplicon libraries. BS and TH analyzed the data. BS, TH, SJ, and KJ wrote the manuscript. All authors read and approved the final manuscript.
Bastiaan Star, Thomas HA Haverkamp contributed equally to this work.