- Research article
- Open Access
Longitudinal study of the early-life fecal and nasal microbiotas of the domestic pig
© Slifierz et al. 2015
- Received: 11 March 2015
- Accepted: 27 August 2015
- Published: 21 September 2015
The mammalian microbiota plays a key role in host health and disease susceptibility. However, knowledge of the early-age microbiota of pigs is lacking. The purpose of this study was to use high-throughput next-generation sequencing to characterize the fecal and nasal microbiotas of pigs during early life.
Ten commercially-raised pigs were randomly enrolled at birth and sampled throughout the first 7 weeks of life. DNA was extracted from fecal and nasal samples and the hypervariable region V4 of the 16S rRNA gene was amplified. The product was sequenced using the Illumina MiSeq platform and 2 × 250 chemistry. Sequencing data was processed and analyzed with the mothur algorithms using an operational taxonomic unit approach. In total, 4.7 million and 5.4 million high-quality sequences were recovered from fecal and nasal samples, respectively. Analysis revealed that these microbiotas contain a very rich and diverse population of bacteria that display a remarkable evolution during the first 7 weeks of life. During this developmental period, a pig was exposed to an average of 1,976 and 6,257 species of bacteria by way of the gastrointestinal and respiratory tracts, respectively. Aging was significantly associated with an increasing measure of richness and diversity as well as with distinct changes to the core microbiota. At 2–3 weeks post-weaning, the rapidly developing microbiotas appeared to reach a developmental milestone as a relative degree of stability was evident.
Pigs are exposed to an incredibly rich and diverse mixture of bacteria during early-life as demonstrated by next-generation sequencing methodology. These findings expand the knowledge of the developing porcine microbiota which is important for understanding susceptibility to disease, particularly for vulnerable neonatal pigs.
- Fecal Microbiota
- Bacterial Phyla
- Community Membership
- Nasal Sample
- Clostridium Cluster XIVa
The microorganisms that populate the body, collectively known as the microbiota, are important drivers of host health and metabolism. Characterization of the porcine microbiota using previous (primarily culture-dependent) techniques provides only a narrow understanding of the complexity of these ecosystems due to methodological limitations [1, 2]. Advancements in next-generation sequencing and bioinformatics have only recently given researchers the opportunity to examine the composition and diversity of these microbial populations in significant detail. Presently, high-throughput sequencing of the hypervariable regions of the 16S rRNA gene can provide a depth of coverage that is far greater than any previous method .
The microbiotas of the gastrointestinal and respiratory tracts of pigs are of particular interest due to the association of these body sites with common swine diseases or pathogens. These microbiotas provide a first line of defense against foreign invaders as competition and interaction between bacteria can protect a host from becoming colonized with particular pathogens . However, the developing microbiota of young pigs is particularly vulnerable to disruptions which can result in long-term impacts that may affect disease susceptibility and growth performance . Such variation in the porcine microbiota has been associated with stress, diet, management practices, and antimicrobial compounds [5–7].
Previous research exploring the gut microbiota in adult pigs using next-generation sequencing methods has revealed an incredibly diverse population of bacteria . The adult porcine gut microbiota has been found to contain at least 7 identifiable bacterial phyla, the predominant being Firmicutes and Bacteroidetes, and at least 171 genera of bacteria . Similarly, the nasal cavity of adult pigs also harbours an incredibly diverse and rich microbial ecosystem which contains an estimated 1,749 species of bacteria from 124 different genera. The porcine nasal microbiota was found to contain 9 phyla of bacteria, of which Proteobacteria, Firmicutes and Spirochaetes predominated .
However, despite the emerging understanding of the microbiota of adult pigs, characterization of the development of the early-age microbiota using next-generation techniques is lacking. Further understanding of the fecal and nasal microbiotas can provide significant insight into swine health and susceptibility to disease, especially among young pigs who are particularly vulnerable. Therefore, the objective of this longitudinal investigation was to characterize the transformation of the fecal and nasal microbiotas of conventionally-raised pigs throughout the first 7 weeks of life using high-throughput next-generation sequencing.
Relative abundance of the top five predominant classes and genera during each period of development
Age Range (days)
Dominant classes (Mean relative abundance, %)
Dominant genera (Mean relative abundance, %)
Clostridium XIVa (4.3)
Clostridium XIVa (7.4)
unclassified Firmicutes (3.4)
unclassified Firmicutes (5.2)
unclassified Firmicutes (4.2)
Erysipelotrichaceae incertae sedis (2.9)
unclassified Firmicutes (3.9)
unclassified Firmicutes (2.9)
Erysipelotrichaceae incertae sedis (3.9)
The genus-level taxonomy of abundant core OTUs of the porcine fecal and nasal microbiotas
Genera of the core OTUs* of the porcine microbiotas (% of pigs with OTU)
Escherichia (Shigella) (100), Clostridium sensu stricto (100)
Globicatella (100), Haemophilus (100), Clostridium sensu stricto (100), Rothia (80), Staphylococcus (80)
Clostridium sensu stricto (100), Clostridium XIVa (80), Haemophilus (80), Lachnospiracea incertae sedis (80)
Globicatella (100), Rothia (100), Haemophilus (80), Clostridium sensu stricto (80)
Clostridium sensu stricto (80), Clostridium XIVa (80), Desulfovibrio (80)
Clostridium sensu stricto (80)
Escherichia (Shigella) (100), Clostridium XIVa (80)
Clostridium sensu stricto (80), Kingella (80), unclassified Firmicutes (80)
Clostridium sensu stricto (100), Kingella (80), unclassified Firmicutes (80)
Megasphaera (80), Lactobacillus (80), Acidaminococcus (80)
Lactobacillus (100), Erysipelotrichaceae incertae sedis (100), Enterococcus (100), Haemophilus (80), Blautia (80)
Megasphaera (100), Lactobacillus (100), Butyricicoccus (100), Erysipelotrichaceae incertae sedis (80), Selenomonas (80), Roseburia (80), Acidaminococcus (80), Faecalibacterium (80), unclassified Firmicutes (80)
Megasphaera (100), Butyricicoccus (100), Erysipelotrichaceae incertae sedis (100), Roseburia (100), Acidaminococcus (100), Faecalibacterium (100), Lactobacillus (80), Selenomonas (80), Streptococcus (80), Prevotella (80), Clostridium XI, unclassified Firmicutes (80)
Lactobacillus (100), Sarcina (100), Erysipelotrichaceae incertae sedis (80), Kingella (80), unclassified Firmicutes (80)
Megasphaera (100), Butyricicoccus (100), Erysipelotrichaceae incertae sedis (100), Roseburia (100), Lactobacillus (80), Streptococcus (80), Prevotella (80), Gemmiger(80)
Lactobacillus (80), unclassified Firmicutes (80)
The most common core OTUs (>1 % relative abundance in ≥80 % of pigs) prior to weaning were from Clostridia (Clostridium sensu stricto and Clostridium cluster XIVa). Interestingly, no core OTUs with >1 % relative abundance could be identified at 21 days of age. However, after weaning, 3 core OTUs were present throughout the entire nursery phase at >1 % relative abundance in at least 80 % of pigs: Megasphaera, Lactobacillus, and Erysipelotrichaceae incertae sedis (Table 2).
A random subsample of 5,381 sequences was used for subsequent analysis. Median sample richness was 1,152 OTUs (Range: 467–5,793), and no age-related association was observed (ρ = 0.15, P = 0.164). The median inverse Simpson’s diversity index was 16.7 (Range: 6.7–72.7) and, unlike richness, there was a positive correlation with age (ρ = 0.22, P = 0.035), although aging only accounted for a small proportion of the variation in diversity. A summary of the number of observed and core OTUs is shown in Fig. 2, and the age-related change in richness and diversity is presented in Fig. 3.
The evolving core nasal microbiota is described in Table 2. During the first 3 days of life, the genera Haemophilus, Globicatella, and Rothia formed the core microbiota in addition to Clostridium sensu stricto which was also present throughout the entire pre-weaning period. Following weaning, Lactobacillus was the only genus present throughout the entire nursery phase at >1 % abundance in at least 80 % of pigs.
The fecal and nasal microbiotas of young pigs display rich and diverse bacterial communities which undergo a very rapid and profound evolution during the first 7 weeks of life. During this developmental period a herd of pigs is exposed to thousands of bacterial species which play an important pioneering role in establishing a stable microbiota that is evident at 2–3 weeks post-weaning, particularly for the fecal microbiota. However, the instability of the microbiota during early life means it may be more susceptible to dysbiosis which can create predispositions to other illnesses as has been reported in human studies . This has created the framework for the role of microorganisms in the ‘programming hypothesis’ which claims that manifestation of diseases later in life can be predetermined by exposure to environmental stimuli during early life . Some research with swine have demonstrated that antimicrobials, excessive hygiene, stress, and management practices have a long-term impact on the gut microbiota [5–7, 11, 12] which may contribute to disease susceptibility or resistance. For example, it has been suggested that the abundance of lactobacilli may influence host immunity and disease susceptibility due to the immunomodulatory activity of these bacteria . Further research exploring the relationship between the microbiota and the early-life programming of health in swine is warranted.
The fecal microbiota described in the present study is similar to previous studies which have found Firmicutes and Proteobacteria to be the predominant phyla pre-weaning [11, 12]. However, studies of the fecal microbiota at 3-months of age demonstrate that a taxonomical shift results in the predominance of Firmicutes and Bacteroidetes [8, 14]. The beginning of this shift was also observed in the present study; Proteobacteria are very dominant during early life (39 % relative abundance) but undergo a considerable decline into adulthood (5 % relative abundance), and the transition into adulthood is marked by an establishment of core OTUs belonging to Bacteroidetes (particularly Prevotella spp.). Proteobacteria, mainly comprising of Escherichia, are common early colonizers of the colon and the early decline in this taxa is possibly due to increased competition from obligate anaerobes (Clostridium, Megasphaera) and colonization resistance due to the immunological properties of colostrum .
Some notable OTUs that were abundant and formed the core microbiota after weaning belonged to butyrate-producing genera including Megasphaera, Butyricicoccus, and Roseburia. Butyrate is an essential energy source for enterocytes and it supports gut health by inhibiting inflammation, reducing oxidative stress, and promoting gut-barrier defense . Additionally, the predominance of Megasphaera after weaning may be related to its extensive metabolic capabilities which include the production of vitamins, amino acids, and short-chain fatty acids .
The microbiota of the anterior nasal passage of young pigs described in the current investigation demonstrates similarities to the nasal microbiota of adult pigs that has previously been studied using next-generation sequencing methods . During early-life, there was an apparent co-dominance between Proteobacteria and Firmicutes, but after weaning there was a noticeable shift in dominance towards Proteobacteria which is also the predominant phylum observed in adult pigs . This transition appeared to coincide with weaning and was marked by an increase in Gammaproteobacteria (particularly Moraxella) and a decrease in Clostridium and Bacilli. The genus Moraxella is of particular interest due to its high abundance in pigs throughout life. Previous research has found that Moraxella is inversely associated with Haemophilus colonization in humans . It is possible that Moraxella plays a protective role in the nasal cavity of pigs by preventing colonization with certain swine pathogens (such as Haemophilus parasuis); although this requires much further research.
Aging appeared to be the most significant driver of development of the fecal and nasal microbiotas, although weaning also played an important role. The reported changes in richness, diversity, relative abundance, and community structure and membership all appeared to follow trends that were independent of weaning, but that is not to say weaning was insignificant. The microbiota did respond to this event with noticeable changes in core OTUs and the detection of specific biomarker OTUs. However, this study was unable to determine whether these changes at weaning were associated with the change in diet, mixing of pigs, stress, or other variables encountered during this event. Furthermore, the rapidly evolving fecal microbiota of young pigs appears to reach a developmental milestone at approximately 2–3 weeks post-weaning. This is demonstrated by the decelerating change in relative abundance of taxa, the clustering of community membership and structure, the plateauing of richness and diversity, and the change in core OTUs. This is less evident for the nasal microbiota, but there are still indications of stability or convergence towards stability during this period.
The present characterization of the early-life microbiota of conventionally-raised pigs using next-generation sequencing provides a broader view of the microbial landscape than previous studies that have employed culture-based methods, fingerprinting, Sanger sequencing, microarrays, and 454-sequencing [1, 2, 5, 11, 12, 14] and is novel in using next-generation sequencing technology to characterize the early-life nasal microbiota, which is poorly described in the literature. However, one of the limitations of the present study is over- or under-representation of particular taxa due to amplification bias. This study also involved only healthy pigs with no clinical presentation of disease and the study was conducted at only one farm, so inter-farm variability and the impact of farm environment and management could not be assessed. Additionally, the fecal microbiota is only representative of the distal gastrointestinal tract (colon, rectum) and characterization of important bacteria populations from the stomach and small intestine may be absent or under-represented.
The fecal and nasal microbiotas of pigs undergo a remarkable evolution during the first 7 weeks of life. The gastrointestinal and respiratory tracts are exposed to thousands of bacteria species with significant turn-over in community membership and structure until a relative degree of stability is evident at 2–3 weeks post-weaning. Future research should explore the significance of variation in the microbial population ecology during early-life and its implication for pathogen shedding and disease susceptibility in both young and adult pigs.
The use of animals in this study was approved by the Animal Care Committee at the University of Guelph. Ten conventionally-raised Yorkshire-Landrace pigs (4 males, 6 females) at a 300-sow batch-farrowing facility were enrolled on the day of birth from ten different litters (1 pig per litter). Each pig in the herd had a unique identifier assigned at birth and the pigs enrolled in the present study were randomly selected (random number generator) from a stratified sampling frame of pigs that weighed ≥1 kg at birth and did not demonstrate clinical signs of disease. Each pig was raised with their mother and litter-mates for 21 days without creep feed and, at weaning, randomized into 4 pens (located in the same nursery room) containing 11–13 pen-mates. Pigs were given a first phase ration for 14 days followed by a second phase ration for the remainder of the nursery period (Additional file 3: Table S1). The rations did not contain antimicrobial compounds and pigs were not administered any antimicrobials during the trial.
Samples were collected on days 1, 3, 7, 14, 21 (prior to weaning), 28, 35, 42, and 49. Samples were collected at the same time of day by the same investigator. At each sampling, a swab of the rectum was collected, as well as whole feces when possible, and a swab of the interior of both nares (approximately 10 mm penetration). Samples were transported at 4 °C and stored at −80 °C before being processed.
Extraction of DNA from fecal and nasal samples was completed using a commercially prepared kit and following the manufacturer’s protocol for stool DNA extraction for pathogen detection (E.Z.N.A. Stool DNA Kit, Omega Bio-Tek Inc., Doraville, Georgia, USA). For nasal samples, the whole tip of the swab was processed through the lysis stage of extraction.
Amplification and sequencing of bacterial 16S rRNA gene
The V4 region of the 16S rRNA gene was amplified using forward (5′-AYTGGGYDTAAAGNG-3′) and reverse (5′-TACNVGGGTATCTAATCC-3′) primers previously designed . The 16S primers contained adapter regions (Forward: TCGTCGGCAGCGTCAGATGTGTATAA-GAGACAG, Reverse: GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG) for annealing to Illumina universal index sequencing adaptors that were added in a later PCR.
Amplification of the 16S rRNA V4 region was completed in a 25 μl reaction consisting of 12.5 μl of KAPA 2G Fast HotStart ReadyMix 2X (KapaBiosystems), 9.0 μl of molecular-grade water, 2.5 μl template DNA, and 0.5 μl each of both the forward and reverse 16S rRNA V4 primers (10.0 μM). The reaction conditions for PCR were 94 °C for 10 min, and 27 cycles of 94 °C for 45 s, 53 °C for 60 s, and 72 °C for 90 s, followed by a final period of 72 °C for 10 min. The product was then purified using AMPure X (Beckman Coulter Inc, Mississauga, Ontario, Canada). The AMPure (20 μl) was mixed with the amplicon (25 μl) and incubated at room temperature for 2 min. After applying a magnetic field, the supernatant was discarded and the beads were washed twice with 80 % ethanol. The beads were then incubated at room temperature for 10 min before eluting with 50 μl of 10 mM Tris buffer.
Illumina universal adapters (Forward: AATGATACGGCGACCACCGAGATCTACAC-index-TCGTCGGCAGCGTC, Reverse: CAAGCAGAAGACGGCATACGAGAT-index-GTCTCGTGGGCTCGG) were then added to the purified 16S rRNA gene product by PCR using a 25 μl reaction consisting of 12.5 μl KAPA 2G Fast HotStart ReadyMix 2X (KapaBiosystems), 8.0 μl of molecular-grade water, 2.5 μl template DNA, and 1.0 μl each of the forward and reverse sample-specific Illumina universal adapters. The PCR conditions were as follows; a single cycle of 94 °C for 3 min, 8 cycles of 94 °C for 45 s, 50 °C for 60 s, and 72 °C for 90 s, followed by a final cycle of 72 °C for 10 min. The product was purified with AMPure X (Beckman Coulter Inc, Mississauga, Ontario, Canada) as previously described and DNA was eluted into 30 μl of 10 mM Tris buffer. The samples were then quantified by spectrophotometry (Nanodrop, Roche, Mississauga, Canada) and normalized to a final concentration of 2 nM. Sequencing of the library pool was performed using an Illumina MiSeq (San Diego, USA) and 2 × 250 chemistry at the University of Guelph’s Advanced Analysis Centre.
Analysis of sequencing data
Sequencing data were analyzed using the mothur software package v.1.33.0 . The paired-end reads were aligned and screened to remove sequences with the following irregularities: contiguous sequence lengths >245 bp or <239 bp, ambiguous base calls, stretches of homopolymers >8 bp, and misalignment with the target region. Sequences were then screened for chimeras using the UCHIME tool  and sequences belonging to non-bacterial domains, including chloroplasts, mitochondria, Archaea and Eukaryotes, were removed. Operational taxonomic units (OTUs) were created using a 3.0 % dissimilarity threshold and average neighbour algorithm, then assigned taxonomy using the RDP reference database (http://rdp.cme.msu.edu). The observed relative abundance of taxa was analyzed prior to subsampling and the mean relative abundance of the 10 pigs was plotted graphically and in a chart to demonstrate changes in phyla, classes, and genera over time. Random subsampling of sequences from each individual sample was completed to normalize the sequence count. Community diversity (inverse Simpson index) and richness (CatchAll; ) were computed and sampling coverage was assessed by Good’s coverage value. The core microbiota (OTUs with >1 % relative abundance and shared amongst ≥80 % of pigs) was also explored for each sample type at different ages. The Jaccard index of dissimilarity was used as an assessment of community membership and the Yue and Clayton measure of dissimilarity was used to as an assessment of community structure. These indices were visualized by principle coordinate analysis, and compared between groups using analysis of molecular variance (AMOVA). Probabilistic modelling using Dirichlet multinomial mixtures (DMM) was utilized in examining community types associated with phase of production . Additionally, the linear discriminant analysis (LDA) effect size (LEfSe) algorithm was used to identify OTUs associated with pre-weaning and post-weaning microbiotas  and cladograms were produced using the online LEfSe tool (http://huttenhower.sph.harvard.edu/galaxy/). Non-parametric statistical analysis (Spearman’s rank correlation) was performed in STATA 10.0 I/C (Stata Corporation, College Station, TX) to determine age-related differences in taxonomical abundance, diversity, and richness. The null hypothesis for all statistical tests was rejected at P < 0.05.
Raw sequences were uploaded to the Sequence Read Archive under study number SRP059335 (experiment number SRX1055388).
We greatly appreciate the assistance of Joyce Rousseau, Amy Sturgeon, Bryan Bloomfield, Glen Cassar, and Karen Richardson. Funding for this project was provided by the Ontario Ministry of Agriculture, Food, and Rural Affairs and Ontario Pork.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Pryde SE, Richardson AJ, Stewart CS, Flint HJ. Molecular analysis of the microbial diversity present in the colonic wall, colonic lumen, and cecal lumen of a pig. Appl Environ Microbiol. 1999;65:5372–7.PubMed CentralPubMedGoogle Scholar
- Leser TD, Amenuvor JZ, Jensen TK, Lindecrona RH, Boye M, Møller K. Culture-independent analysis of gut bacteria: the pig gastrointestinal tract microbiota revisited. Appl Environ Microbiol. 2002;68:673–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Shokralla S, Spall JL, Gibson JF, Hajibabaei M. Next‐generation sequencing technologies for environmental DNA research. Mol Ecol. 2012;21:1794–805.View ArticlePubMedGoogle Scholar
- Pettigrew MM, Gent JF, Revai K, Patel JA, Chonmaitree T. Microbial interactions during upper respiratory tract infections. Emerg Infect Dis. 2008;14:1584–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Thompson CL, Wang B, Holmes AJ. The immediate environment during postnatal development has long-term impact on gut community structure in pigs. ISME J. 2008;2:739–48.View ArticlePubMedGoogle Scholar
- Isaacson R, Kim HB. The intestinal microbiome of the pig. Anim Health Res Rev. 2012;13:100–9.View ArticlePubMedGoogle Scholar
- Weese JS, Slifierz M, Jalali M, Friendship R. Evaluation of the nasal microbiota in slaughter-age pigs and the impact on nasal methicillin-resistant Staphylococcus aureus (MRSA) carriage. BMC Vet Res. 2014;15:69.View ArticleGoogle Scholar
- Kim HB, Borewicz K, White BA, Singer RS, Sreevatsan S, Tu ZJ, et al. Longitudinal investigation of the age-related bacterial diversity in the feces of commercial pigs. Vet Microbiol. 2011;153:124–33.View ArticlePubMedGoogle Scholar
- Munyaka PM, Khafipour E, Ghia JE. External influence of early childhood establishment of gut microbiota and subsequent health implications. Front Pediatr. 2014;2:109.View ArticlePubMed CentralPubMedGoogle Scholar
- Azad MB, Kozyrskyj AL. Perinatal programming of asthma: the role of gut microbiota. Clin Dev Immunol. 2012;2012:932072.View ArticlePubMed CentralPubMedGoogle Scholar
- Schmidt B, Mulder IE, Musk CC, Aminov RI, Lewis M, Stokes CR, et al. Establishment of normal gut microbiota is compromised under excessive hygiene conditions. PloS One. 2011;6, e28284.View ArticlePubMed CentralPubMedGoogle Scholar
- Schokker D, Zhang J, Vastenhouw SA, Heilig HG, Smidt H, Rebel JMJ, et al. Long-lasting effects of early-life antibiotic treatment and routine animal handling on gut microbiota composition and immune system in pigs. PloS One. 2015;10, e0116523.View ArticlePubMed CentralPubMedGoogle Scholar
- Schokker D, Zhang J, Zhang L, Vastenhouw SA, Heilig HGHJ, Smidt H, et al. Early-life environmental variation affects intestinal microbiota and immune development in new-born piglets. PloS One. 2014;9, e100040.View ArticlePubMed CentralPubMedGoogle Scholar
- Looft T, Allen HK, Cantarel BL, Levine UY, Bayles DO, Alt DP, et al. Bacteria, phages and pigs: the effects of in-feed antibiotics on the microbiome at different gut locations. ISME J. 2014;8:1566–76.View ArticlePubMedGoogle Scholar
- Fairbrother JM, Nadeau E, Gyles CL. Escherichia coli in postweaning diarrhea in pigs: an update on bacterial types, pathogenesis, and prevention strategies. Anim Health Res Rev. 2005;6:17–39.View ArticlePubMedGoogle Scholar
- Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost FJ, Brummer RJ. The role of butyrate on colonic function. Aliment Pharmacol Ther. 2008;27:104e19.Google Scholar
- Shetty SA, Marathe NP, Lanjekar V, Ranade D, Shouche YS. Comparative genome analysis of Megasphaera sp. reveals niche specialization and its potential role in the human gut. PloS One. 2013;8:e79353.View ArticlePubMed CentralPubMedGoogle 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 Methods. 2010;7:335–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.View ArticlePubMed CentralPubMedGoogle Scholar
- Bunge J, Woodard L, Bohning D, Foster JA, Connolly S, Allen HK. Estimating population diversity with CatchAll. Bioinformatics. 2012;28:1045–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Holmes I, Harris K, Quince C. Dirichlet multinomial mixtures: generative models for microbial metagenomics. PLoS One. 2012;7, e30126.View ArticlePubMed CentralPubMedGoogle Scholar
- Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.View ArticlePubMed CentralPubMedGoogle Scholar