Diet shapes the gut microbiome of pigs during nursing and weaning
© Frese et al. 2015
Received: 13 January 2015
Accepted: 6 June 2015
Published: 1 July 2015
The newborn mammal is rapidly colonized by a complex microbial community, whose importance for host health is becoming increasingly clear. Understanding the forces that shape the early community, especially during the nursing period, is critical to gain insight into how this consortium of microbes is assembled. Pigs present an attractive model for nursing humans, given physiological and compositional similarity of pig and human milk and the utility of pigs in experimental studies. However, there is a paucity of data examining the gut microbiome in nursing pigs from birth through weaning using modern molecular methods and fewer experimental studies that examine the impact of diet on these microbial communities.
We characterized the fecal microbiome of pigs from birth through 7 weeks of age, during which the animals were transitioned from an exclusive diet of sow milk to a starter diet composed of plant and animal-based components. Microbial communities were clearly distinguishable based on diet, being relatively stable absent dietary changes. Metagenomic sequencing was used to characterize a subset of animals before and after weaning, which identified glycan degradation pathways differing significantly between diets. Predicted enzymes active on milk-derived glycans that are otherwise indigestible to the host animal were enriched in the microbial metagenome of milk-fed animals. In contrast, the bacterial metagenome of weaned animals was enriched in functional pathways involved in plant glycan deconstruction and consumption.
The gut microbiome in young pigs is dramatically shaped by the composition of dietary glycans, reflected by the different functional capacities of the microbiome before and after weaning.
The mammalian gastrointestinal tract is a home to a complex and diverse microbial community that profoundly influences health and disease. The communities that assemble in this ecosystem, beginning shortly after birth, have been shown to be remarkably host-specific and broadly stable over time [1, 2]. Studies of the gut microbial community illustrate how populations of constituents are shaped by environmental exposure to microbes, diet, immunological pressures, host genetics, and ecological forces within the ecosystem itself [3–5]. The dramatic importance of this community for host health and resistance to disease has led to the evolution of elaborate mechanisms to facilitate the acquisition of this community across generations.
In shaping these populations, milk has been shown to provide an important selective advantage to some microbes, allowing them to dominate an ecosystem [6–8]. This has been best described for humans, but the secreted glycans found in human milk (human milk oligosaccharides (HMOs)) are structurally similar to those found in bovine and porcine milks [9, 10]. Understanding which microbial populations are influenced by milk and their fate at weaning provides insight into how unique, and powerful, dietary changes in the lives of mammals shape the gut microbiome. It has been hypothesized that milk, especially milk glycans, have a controlling effect on the microbiome of nursing mammals, resulting in a “milk-oriented microbiome” (MOM) . However, it is poorly understood how exposure and diet interact to shape the microbiome, that is, whether stochastic exposure or the diet, in this case milk, drives the major populations and changes observed in the microbiome early in life.
Few studies examine the microbiome of nursing pigs, and no current studies follow animals from birth past weaning longitudinally. In this study, we followed the microbial communities in nursing pigs from birth through weaning, up to 7 weeks of age. Using 16S rRNA marker gene sequencing, we found that diet, rather than litter effects or age, best described the differences between samples. Metagenomic sequencing of microbial communities confirmed this distinction, discriminating samples starkly by functional capacities related to dietary glycans, found in the animal’s respective diets. By the end of our study, the fecal microbiome of animals studied here closely resembled that of age-matched animals in other studies, suggesting that early stochasticity is replaced by later convergence and shaped by diet [12–15].
Community diversity increases over time but populations shift with diet
Functional composition of the metagenome reshaped by diet
The microbiome is rapidly assembled in pigs
Through exposure after parturition, the gut microbiome develops rapidly in a mammal’s early life. Large taxonomic shifts have been observed in humans through the first years of life [16–19], and here we report that microbial turnover in pigs early in life is also dramatic as the animals shift to an adult-like diet. Unlike humans, the weaning process for pigs is an abrupt dietary shift, from solely sow milk to a complete feed based diet, that produces similarly abrupt taxonomic and functional shifts in the gut microbiome (Fig. 3). While there is wide variation, human infants are generally introduced to novel foods gradually, the result of which is rapidly reflected in the gut microbiome . Interestingly, early colonization of these animals did not include a “succession” of taxa, as has been reported for early human colonizers, where Lactobacillus or Proteobacteria are early, transient colonizers and are supplanted by more stable populations of Bifidobacterium or Bacteroides [20–22]. Instead, in pigs, we observe dominant, stable populations from the first day after birth whose principle composition remain for the first 3 weeks of life, while minor populations begin to colonize, increasing overall diversity over time (Figs. 1 and 3). While in humans, much attention has focused on western populations and how delivery mode shapes the infant gut microbiome, some human populations show a much less dramatic effect of delivery mode . Our results in pigs suggest that exposure to these microbes may seed these early colonizers, who are present in low abundances but unable to consume available substrates, rather than a coordinated ecological succession during nursing.
In these animals, the nursing period is clearly dominated by a milk-oriented microbiome; a community whose functional genetic content reflects its focus on the consumption of milk (Fig. 4b, Additional file 2: Figure S2A). Other recent studies on pigs have found comparable trends, composed of the same taxa, suggesting the broad applicability of these results among animals in the US, South Korea, and France [13, 15]. For example, one of the most dramatic changes observed was in a small population of Prevotella that is present in nursing animals from birth. This population is very low in relative abundance and increases rapidly and dramatically when a more favorable diet is introduced at weaning (Fig. 3, Additional file 1: Figure S1). This genus is widely associated with plant polysaccharide consumption and is enriched by plant polysaccharides , which are found in the diet after weaning but not in sow milk and typifies the dramatic taxonomic shifts associated with weaning. Indeed, these Prevotellaceae appear to supplant Bacteroides populations in nursing pigs, whose ability to harvest milk glycans is well characterized [25, 26] and has been hypothesized to be linked to the dietary dichotomy we describe here .
The milk-oriented microbiome harvests dietary milk glycans
Milk glycans are a structurally and compositionally diverse and abundant carbon source in milk, best described for humans , but present in the milk of many different mammals . It has been suggested that these glycans have a profound impact on the gut microbiome, shaping a composition described as the “milk-oriented microbiome” or “MOM” which contributes to neonatal health . Porcine milk is composed primarily of N-acetylglucosamine, sialic acids (N-acetylneuraminic or N-glycolneuraminic acid), galactose and glucose monomers, and less abundantly, fucose. These monomers are linked together to form trimers, tetramers, and more complex bi-antennary structures . While structurally less diverse than human milk oligosaccharides, these compounds are compositionally similar to both human and bovine milk oligosaccharides and require a similar array of enzymatic capacities to deconstruct these structures. As the host lacks the metabolic capacity to consume these highly sialylated complex carbohydrates, they pass undigested to the distal GIT, where they shape the composition of the gut microbiome. As the host digests and absorbs simple sugars proximally to the distal gut, only glycans resistant to digestion pass through to influence the distal gut microbiota. In this study, we found that the enzymatic capacities necessary to remove and deconstruct milk sugar components were enriched, relative to the community found in weaned animals, where enzymes predicted to encode plant glycan hydrolases were increased (Fig. 4c, Additional file 2: Figure S2B). This implies also that mucus glycans, which are structurally similar to milk glycans, play a very minor role, if any, in shaping the gut microbiome of weaned pigs. Further, metabolic pathways for the catabolism of galactose, lactose, N-acetylglucosamine, and sialic acid were also enriched (Additional file 2: Figure S2A), suggesting that the gut microbiome of nursing pigs is both enzymatically and metabolically oriented to the consumption of milk oligosaccharides, as is the case for infant-associated microbes in humans .
The weaned microbiota harvest plant glycans
The abrupt weaning methods used in pig rearing allow for an excellent experimental design to compare the effects of diet. On the withdrawal of milk, we see a dramatic change in microbial populations. Specifically, in the abundances of Lactobacillus, a taxon that consumes plant-derived mono- and di-saccharides and milk sugars such as lactose but not complex milk sugars [28, 29]. Populations of saccharolytic microbes like Prevotellaceae increased from less than 0.3 % to more than 15 % of the total community, whose apparent comparative advantage for these substrates over Bacteroidaceae contributed to the concomitant decrease of Bacteroidaceae (Fig. 3). Increased relative abundances of β-glucanases, cellulases, xylanases, and cellobiases and associated downstream catabolic pathways all point to dramatic shifts in response to the compositional changes of diet, rather than exposure to new microbes from food (Fig. 4, Additional file 2: Figure S2). Of note, while the weaning feed contained antibiotics (tiamulin and chlortetracycline), the taxonomic composition of animals in this study resembled that of animals at similar ages in other studies with and without antibiotics [12, 13, 30], which indicates broad similarities in taxonomic and functional  composition, though antibiotic resistance may vary, in the microbiome of animals across different sites and production scenarios.
Recent work has cast an important light on how the gut microbiome responds to changes in diet, using mouse models and human feeding trials [5, 32]. Here, we use a simple and highly relevant diet (milk) and a relevant and highly tractable animal model that resembles human digestion (pigs), whose diet composition is changed abruptly (from milk to near-adult feed) as a matter of course in rearing. The advantages in this model over other studies are evident in our level of detail at dissecting catabolic pathways using known dietary glycan structural data as a starting point. Our observations are in agreement with previous work in both pigs [12, 13] and human dietary studies investigating differences in animal and plant-based diets . We also show the impact of abrupt and distinct changes in diet on the distal gut microbiome of pigs, at both a taxonomic and metagenomic level. Given that substrate availability was found to shape both taxonomic structure and enzymatic functional capacities, it is interesting to consider whether biogeographical variations in the gut microbiome and local variations in substrate availability may also be associated with similar taxonomic and functional differences. Future work may shed light in this respect on the role of the microbiome in digestion and whether specific pig milk oligosaccharides (PMOs) shape populations of microbes in pigs as HMOs shape the gut of human infants.
All experiments involving animals were reviewed and approved by the University of California Davis Institutional Animal Care and Use Committee prior to beginning of the experiment (approval #17776). Throughout the study, all animals were housed in a controlled-access specific pathogen free facility at the University of California Davis. Three healthy Yorkshire/Hampshire adult pregnant multiparous sows from the University of California herd were selected for this study. Upon delivery, the infant pigs (n = 27) were cohoused with sows by litter and ear notched for individual identification and dosed with 50 mg spectinomycin and Penicillin G as benzathine and procaine (75,000 IU each), following standard husbandry practices for swine. The infant pigs were allowed to nurse freely until weaning after 21 days of age without creep feed and did not, to our knowledge, consume sow feed. Nursing pigs were removed from the sow and transferred to separate housing and fed a standard starter feed (Hubbard Feeds Mankato, MN USA) after 21 days of age. The weaning feed was oat-based and also contained 444 mg/kg chlortetracycline and 39 mg/kg tiamulin, while the sow diet was unmedicated. Animals were given ad libitum access to water and feed.
Fecal samples were collected rectally from each animal using a sterile cotton swab (Puritan Medical, Guilford, ME USA) wet with sterile phosphate buffered saline (pH 7.0) 1, 3, 5, 7, 14, 21, 28, 35, and 42 days after birth.
Sequencing library construction
DNA was extracted from swabs using the Zymo Research Fecal DNA kit (ZYMO Research Irvine, CA USA) according to the manufacturer’s instructions. Extracted DNA was used as a template for PCR using barcoded primers to amplify the V4 region of the 16S rRNA gene as previously described [34, 35]. Briefly, the V4 domain of the 16S rRNA gene was amplified using primers F515 (5′-NNNNNNNN GTGTGCCAGCMGCCGCGGTAA-3′) and R806 (5′-GGACTACHVGGGTWTCTAAT-3′), where the poly-N (italicized) sequence was an 8-nt barcode unique to each sample and a 2-nt linker sequence (bold). PCR amplification was carried out in a 15-μL reaction containing 1× GoTaq Green Mastermix (Promega, Madison, WI USA), 1 mM MgCl2 and 2 pmol of each primer. The amplification conditions included an initial denaturation step of 2 min at 94 °C, followed by 25 cycles of 94 °C for 45 s, 50 °C for 60 s, and 72 °C for 90 s, followed by a single final extension step at 72 °C for 10 min. Amplicons were pooled and purified using a Qiagen PCR purification column and submitted to the UC Davis Genome Center DNA Technologies Sequencing Core for paired-end library preparation, cluster generation, and 250-bp read sequencing on an Illumina MiSeq.
Quality-filtering of de-multiplexed reads was conducted as recommended  by removing low-abundance OTUs, to remove spurious reads and OTUs which result from artifacts of sequencing, and data was analyzed using QIIME 1.8.0 . The 13_8 GreenGenes database release was used for open reference OTU picking and taxonomy assignment, and bacterial sequences were aligned using UCLUST . To ensure even sequencing depth across samples, 7,000 sequences per sample were randomly subsampled for analysis of bacterial communities. Samples with fewer than 7,000 sequences were omitted. Alpha diversity estimates were computed for phylogenetic diversity (PD) and compared  by nonparametric two-sample t test with Bonferroni correction and 999 Monte Carlo permutations. Analysis of similarities (ANOSIM) distances were calculated using compare_categories.py using a weighted UNIFRAC distance matrix calculated by beta_diversity.py, both using QIIME . Family-level taxa abundances were used in STAMP  to generate principal component analyses (PCA). All 16S rRNA sequencing data is publicly available in QIITA (http://qiita.microbio.me).
Total genomic DNA was extracted from fecal samples with the ZYMO Research Fecal DNA Extraction kit according to manufacturer’s instructions and prepared using the Illumina MiSeq v3 Reagent Chemistry for whole genome shotgun sequencing of multiplexed 150-bp libraries at the University of California Davis Genome Sequencing Core (http://dnatech.genomecenter.ucdavis.edu). Samples used for metagenomic sequencing were selected from one animal from each litter, sampled before and after weaning, on day 14 (nursing, n = 3) and day 35 (weaned, n = 3). Barcoded sample libraries were pooled and sequenced across triplicate, single-read sequencing runs to reach adequate sequencing depth for fecal samples  and prevent run-to-run variation. FASTQ files from each run were de-multiplexed, quality filtered, and trimmed to 150 bp, and then reads for each sample were pooled from the three runs, yielding 15–20 million reads per sample, and submitted to the Metagenomics Rapid Annotation using Subsystem Technology (MGRAST) pipeline for analysis , which removes host genomic DNA reads and duplicate reads, bins 16S rRNA reads, and functionally classifies remaining reads by predicted protein sequence using several different schemes. All metagenomic sequencing data is publicly available on the MGRAST website (http://metagenomics.anl.gov). Subsystem classified reads were functionally classified and normalized in MGRAST and compared between groups by multiple t test using GraphPad Prism 6 (GraphPad Software Inc, La Jolla, CA).
We thank Nicholas Bokulich, Danielle Lemay, Morgan Lee, and Karen Kalanetra for the helpful discussions and guidance in DNA sequencing methodologies and analysis. We are also thankful to the UC Davis DNA Technologies Core Facility for their sequencing services. Metagenomic sequencing data is publicly available in the MGRAST database (http://metagenomics.anl.gov), and 16S marker gene sequencing reads, metadata, and sample barcodes are publicly available in QIITA (http://qiita.microbio.me).
DAM is a co-founder of Evolve Biosystems, a company focused on diet-based manipulation of the gut microbiota.
- Martínez I, Muller CE, Walter J. Long-term temporal analysis of the human fecal microbiota revealed a stable core of dominant bacterial species. PLoS One. 2013;8, e69621.PubMed CentralPubMedView ArticleGoogle Scholar
- Faith JJ, Guruge JL, Charbonneau M, Subramanian S, Seedorf H, Goodman AL, et al. The long-term stability of the human gut microbiota. Science. 2013;341:1237439–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Benson AK, Kelly SA, Legge R, Ma F, Low SJ, Kim J, et al. Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. PNAS. 2010;107(44):18933–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Round JL, Mazmanian SK. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol. 2009;9:313–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Leamy LJ, Kelly SA, Nietfeldt J, Legge RM, Ma F, Hua K, et al. Host genetics and diet, but not immunoglobulin A expression, converge to shape compositional features of the gut microbiome in an advanced intercross population of mice. Genome Biol. 2014;15:811.View ArticleGoogle Scholar
- Marcobal A, Barboza M, Froehlich JW, Block DE, German JB, Lebrilla CB, et al. Consumption of human milk oligosaccharides by gut-related microbes. J Agric Food Chem. 2010;58:5334–40.PubMed CentralPubMedView ArticleGoogle Scholar
- Ward RE, Niñonuevo M, Mills DA, Lebrilla CB, German JB. In vitro fermentability of human milk oligosaccharides by several strains of bifidobacteria. Mol Nutr Food Res. 2007;51:1398–405.PubMedView ArticleGoogle Scholar
- Sela DA, Mills DA. Nursing our microbiota: molecular linkages between bifidobacteria and milk oligosaccharides. Trends Microbiol. 2010;18:298–307.PubMed CentralPubMedView ArticleGoogle Scholar
- Tao N, Ochonicky KL, German JB, Donovan SM, Lebrilla CB. Structural determination and daily variations of porcine milk oligosaccharides. J Agric Food Chem. 2010;58:4653–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Tao N, Wu S, Kim J, An HJ, Hinde K, Power ML, et al. Evolutionary glycomics: characterization of milk oligosaccharides in primates. J Proteome Res. 2011;10:1548–57.PubMed CentralPubMedView ArticleGoogle Scholar
- Zivkovic AM, Lewis ZT, German JB, Mills DA. Establishment of a milk-oriented microbiota (MOM) in early life: how babies meet their MOMs. Functional Food Rev. 2013;5:3–12.Google Scholar
- Kim HB, Borewicz K, White BA, Singer RS. Longitudinal investigation of the age-related bacterial diversity in the feces of commercial pigs. Vet Microbiol. 2011;153(1):124–33.PubMedView ArticleGoogle Scholar
- Pajarillo E, Chae JP, MP B, HB K, Kang DK. Assessment of fecal bacterial diversity among healthy piglets during the weaning transition. J Gen Appl Microbiol. 2014;60:140–6.View ArticleGoogle Scholar
- Poroyko V, White JR, Wang M, Donovan S, Alverdy J, Liu DC, et al. Gut microbial gene expression in mother-fed and formula-fed piglets. PLoS One. 2010;5, e12459.PubMed CentralPubMedView ArticleGoogle Scholar
- Mach N, Berri M, Estelle J, Levenez F, Lemonnier G, Denis C, et al. Early life establishment of the swine gut microbiome and impact on host phenotypes. Environ Microbiol Rep. 2015;7(3):1–33.View ArticleGoogle Scholar
- Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486(7402):222–7.PubMed CentralPubMedGoogle Scholar
- Subramanian S, Huq S, Yatsunenko T, Haque R, Mahfuz M, Alam MA, et al. Persistent gut microbiota immaturity in malnourished Bangladeshi children. Nature. 2014;510:417–21.PubMed CentralPubMedGoogle Scholar
- Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, et al. Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci. 2011;108:4578–85.PubMed CentralPubMedView ArticleGoogle Scholar
- Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO. Development of the human infant intestinal microbiota. Plos Biol. 2007;5, e177.PubMed CentralPubMedView ArticleGoogle Scholar
- Orrhage K, Nord CE. Factors controlling the bacterial colonization of the intestine in breastfed infants. Acta Paediatr Suppl. 1999;88:47–57.PubMedView ArticleGoogle Scholar
- Dominguez-Bello MG, Costello EK. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. PNAS. 2010;107(26):11971–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Favier CF, Vaughan EE, de Vos WM, Akkermans ADL. Molecular monitoring of succession of bacterial communities in human neonates. Appl Environ Microbiol. 2002;68:219–26.PubMed CentralPubMedView ArticleGoogle Scholar
- Huda MN, Lewis Z, Kalanetra KM, Rashid M, Ahmad SM, Raqib R, et al. Stool microbiota and vaccine responses of infants. Pediatrics. 2014;134:e362–72.PubMedView ArticleGoogle Scholar
- Ivarsson E, Roos S, Liu HY, Lindberg JE. Fermentable non-starch polysaccharides increases the abundance of Bacteroides–Prevotella–Porphyromonas in ileal microbial community of growing pigs. Animal. 2014;8(11):1777–87.PubMedView ArticleGoogle Scholar
- Marcobal A, Barboza M, Sonnenburg ED, Pudlo N, Martens EC, Desai P, et al. Bacteroides in the infant gut consume milk oligosaccharides via mucus-utilization pathways. Cell Host Microbe. 2011;10(5):507–14.PubMed CentralPubMedView ArticleGoogle Scholar
- Marcobal A, Sonnenburg JL. Human milk oligosaccharide consumption by intestinal microbiota. Clin Microbiol Infect. 2012;18 Suppl 4:12–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Zivkovic AM, German JB, Lebrilla CB, Mills DA. Human milk glycobiome and its impact on the infant gastrointestinal microbiota. Proc Natl Acad Sci. 2011;108:4653–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Ward RE, Niñonuevo M, Mills DA, Lebrilla CB, German JB. In vitro fermentation of breast milk oligosaccharides by Bifidobacterium infantis and Lactobacillus gasseri. Appl Environ Microbiol. 2006;72:4497–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Schwab C, Gänzle M. Lactic acid bacteria fermentation of human milk oligosaccharide components, human milk oligosaccharides and galactooligosaccharides. FEMS Microbiol Lett. 2011;315:141–8.PubMedView ArticleGoogle Scholar
- Looft T, Johnson TA, Allen HK, Bayles DO, Alt DP, Stedtfeld RD, et al. In-feed antibiotic effects on the swine intestinal microbiome. PNAS. 2012;109(5):1691–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Lamendella R, Domingo JWS, Ghosh S, Martinson J, Oerther DB. Comparative fecal metagenomics unveils unique functional capacity of the swine gut. BMC Microbiol. 2011;11:103.PubMed CentralPubMedView ArticleGoogle Scholar
- Carmody RN, Gerber GK, Luevano JM, Gatti DM, Somes L, Svenson KL, et al. Diet dominates host genotype in shaping the murine gut microbiota. Cell Host Microbe. 2015;17(1):72–84.PubMedView ArticleGoogle Scholar
- David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2013;505(7484):1–18.View ArticleGoogle Scholar
- Bokulich NA, Joseph CML, Allen G, Benson AK, Mills DA. Next-generation sequencing reveals significant bacterial diversity of botrytized wine. PLoS One. 2012;7, e36357.PubMed CentralPubMedView ArticleGoogle Scholar
- Choy YY, Quifer-Rada P, Holstege DM, Frese SA, Calvert CC, Mills DA, et al. Phenolic metabolites and substantial microbiome changes in pig feces by ingesting grape seed proanthocyanidins. Food and Function. 2014;5(9):2298–308.PubMedView ArticleGoogle Scholar
- Bokulich NA, Subramanian S, Faith JJ, Gevers D. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Meth. 2013;10:57–9.View ArticleGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Meth. 2010;7:335–6.View ArticleGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.PubMedView ArticleGoogle Scholar
- Faith DP. Conservation evaluation and phylogenetic diversity. Biol Conserv. 1992;61:1–10.View ArticleGoogle Scholar
- Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4.PubMedView ArticleGoogle Scholar
- Rodriguez-R LM, Konstantinidis KT. Estimating coverage in metagenomic data sets and why it matters. ISME J. 2014;8:2349–51.PubMedView ArticleGoogle Scholar
- Glass EM, Wilkening J, Wilke A, Antonopoulos D, Meyer F. Using the metagenomics RAST server (MG-RAST) for analyzing shotgun metagenomes. Cold Spring Harb Protoc. 2010;2010(1):pdb–prot5368.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.