- Research article
- Open Access
Next generation sequencing shows high variation of the intestinal microbial species composition in Atlantic cod caught at a single location
- Bastiaan Star†1Email author,
- Thomas HA Haverkamp†1Email author,
- Sissel Jentoft1 and
- Kjetill S Jakobsen1
© Star et al.; licensee BioMed Central Ltd. 2013
Received: 11 March 2013
Accepted: 7 November 2013
Published: 9 November 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
Alpha diversity estimates of the Atlantic cod intestinal microbial community
Inverse Simpsons index
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 .
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.
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- The Human Microbiome Consortium: Structure, function and diversity of the healthy human microbiome. Nature. 2012, 486 (7402): 207-214. 10.1038/nature11234.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Nelson JS: Fishes of the world. 2006, Hoboken, New Yersey: John Wiley & SonsGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Nayak SK: Role of gastrointestinal microbiota in fish. Aquac Res. 2010, 41 (11): 1553-1573. 10.1111/j.1365-2109.2010.02546.x.View ArticleGoogle Scholar
- Waters JM, Fraser CI, Hewitt GM: Founder takes all: density-dependent processes structure biodiversity. TREE. 2013, 28 (2): 78-85.PubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.