Human milk microbiota profiles in relation to birthing method, gestation and infant gender
© Urbaniak et al. 2015
Received: 22 June 2015
Accepted: 27 November 2015
Published: 6 January 2016
Human milk is an important source of bacteria for the developing infant and has been shown to influence the bacterial composition of the neonate, which in turn can affect disease risk later in life. Very little is known about what factors shape the human milk microbiome. The goal of the present study was to examine the milk microbiota from a range of women who delivered vaginally or by caesarean (C) section, who gave birth to males or females, at term or preterm.
Milk was collected from 39 Caucasian Canadian women, and microbial profiles were analyzed by 16S ribosomal RNA (rRNA) sequencing using the Illumina platform.
A diverse community of milk bacteria was found with the most dominant phyla being Proteobacteria and Firmicutes and at the genus level, Staphylococcus, Pseudomonas, Streptococcus and Lactobacillus. Comparison of bacterial profiles between preterm and term births, C section (elective and non-elective) and vaginal deliveries, and male and female infants showed no statistically significant differences.
The study revealed the diverse bacterial types transferred to newborns. We postulate that there may be a fail-safe mechanism whereby the mother is “ready” to pass along her bacterial imprint irrespective of when and how the baby is born.
With the incidence of various non-infectious diseases on the rise, there is much interest in the developmental origins of health and disease and the potential role of early life feeding practices in modulating these outcomes. Breast-fed infants have been shown to be better protected than formula-fed infants against necrotizing enterocolitis and diarrhoea, allergy and asthma, inflammatory bowel disease, type 1 and type II diabetes, obesity and cardiovascular disease [1, 2]. In addition to immune protection and bioactive compounds being conveyed through maternal milk, a possible protective role of bacteria has been suggested. Lower than average levels of Bifidobacterium in human milk correlate with low levels of Bifidobacterium in the neonatal gut , allowing for higher than normal levels of Bacteroides to be established . These high levels of Bacteroides early in life have been associated with an increased risk of asthma and obesity later in life [5–7]. Indeed, efforts to manipulate the microbiota of formula-fed infants through probiotic supplementation have resulted in protection against some of the above diseases, comparable to that observed for breast-fed infants [8–10].
Differences exist in bioactive components, macronutrients, cytokines, enzymes, proteins and immunological factors between preterm and term milk and milk from mothers giving birth by vaginal and caesarean deliveries [11–16]. As well, the energy content differs in milk depending on gender of the newborn, with breast milk from mothers who give birth to sons having more fat content than that of daughters [17, 18].
We hypothesized that physiological or hormonal triggers that influence milk composition might also support different bacterial genera. Thus, we studied human milk samples from women giving birth at different stages of gestation, by vaginal or caesarean delivery, and examined whether gender of the newborn also affected the microbiota profiles.
The dataset was also analyzed using ALDEx2  to examine whether specific taxa were differentially expressed based on gestation, mode of delivery or gender. Again, no differences were detected at the genus (Additional file 6: Table S3), family, class or phyla levels (data not shown).
Of interest, the abundance profiles clearly showed that mothers are transmitting very different bacterial profiles to newborns. For example, the milk of subject 40 contained over 80 % abundance of staphylococci, whereas that of subject 25 comprised mostly Gram-negative organisms, especially Pseudomonas.
This study revealed a range of bacterial genera in human milk, consistent with previous studies [20–22]. Interestingly, even when a baby was born extremely prematurely (subject S24), the mother’s milk was similar in composition to that of a woman giving birth at full term (subject S30). It can be speculated that this might be a fail-safe mechanism whereby the mother is “ready” to pass along her bacterial imprint irrespective of when the baby is born. If so, the microbiota would appear to be recalcitrant to late pregnancy hormonal and inflammatory changes, which could indicate an evolutionary pressure protecting this niche for the baby’s benefit. Further studies on this concept are warranted, and cases where the milk profiles are very different or the outcome for the baby negative could be particularly insightful.
It was not surprising that milk from emergency caesarean (C) section deliveries (i.e. non-elective) did not differ from women who gave birth vaginally, as this decision is made at the time of labour, when the hormones and timing for birth have been initiated. We also did not see differences between non-elective C section deliveries and vaginal births which is in contrast to a study published by Cabrera-Rubio et al. in which they concluded that the human milk microbiome is shaped by mode of delivery . However, we do not believe that the authors adequately proved this claim. In their analysis, the milk samples collected were from both obese and normal weight women with no indication of the proportion in each group. Since it was shown in that same paper that body mass index influences the milk microbiota, the subject’s weight could have been a confounding factor responsible for the observed differences . In addition, the author’s claim that the milk microbiota is influenced by mode of delivery was based solely on visual observations of bar plots (which were not very distinct between the two groups), with no principal coordinate analysis (PCoA) plots or statistical analyses to support this claim. These conflicting results between our group and that of Cabrera-Rubio highlight the need for future studies with larger sample sizes and inclusion of women from various demographics.
Because of the multivariate nature of the data, differences, if they are present, may have been masked by the different variables confounding each other. Thus, a larger sample size allowing for linear regression analysis would strengthen our claims. However, there is the possibility that no matter the sample size, differences will only be apparent when examined at the level of the individual. Schwarzberg et al.  showed that treatment for periodontitis had no effects on microbial profiles when compared to controls. However, when pre- and post-treatment effects for each individual were compared, significant changes in bacterial composition were observed. The authors thus emphasized the importance of comparing shifts from a personalized healthy state to a personalized disease state in order to truly understand biological changes .
Proteobacteria and Firmicutes were the dominant phyla, consistent with a previous high-throughput study performed in Switzerland and analyzed with a different sequencing platform (454 sequencing) .
The detection of Staphylococcus, as the most abundant organism, is likewise consistent with other studies [20–22]. Staphylococcus is more abundant in the gut of breast-fed newborns compared to formula-fed ones [25, 26], but their numbers start to decrease after the first week of life when oxygen has been consumed and an environment favourable for anaerobes is created [4, 27]. Unlike Staphylococcus, which is present in high abundance in at least the first week of life, Proteobacteria are never found in high numbers in the faeces of newborns, infants or Western adults [4, 28, 29]. Thus, an important question arises, as to what the evolutionary significance is of having such a diverse population of bacteria in milk, if only a select few seem to colonize during the development of the neonate. There are a few possibilities; firstly, persistent colonization is not always needed for beneficial effects, as transient exposure can be just as effective [30, 31]. Secondly, bacteria in milk may not have to be passed on from the mother to the child to exert their beneficial effects. Many of the protective factors in milk such as antibodies, immune cells, lactoferrin and beta defensins originate from the mammary gland and not from the blood . With the ability of bacteria to regulate host gene expression, such as anti-microbials, and their ability to stimulate the immune system, the plethora of bacteria in breast milk could be inducing up-regulation of these protective factors in the breast that then get passed on in high concentrations to the neonate via milk.
From another perspective, there is the possibility that some of these milk microbes have either limited or no effects on the offspring but are present for the benefit of the mother, such as in the protection against mastitis. Mastitis is a painful inflammatory condition of the breast with the main causative agent being Staphylococcus aureus, and it has been shown that some milk commensals have the ability to prevent S .aureus growth and infection . The same is true for a skin derived strain of Propionibacterium acnes, which prevents growth of S. aureus via its by-products of glycerol fermentation . With the high abundance of glycerol present in milk, milk-derived strains of Propionibacterium may have similar growth-inhibiting properties, which would account for its presence in every milk sample collected.
Lactobacillus was present in high abundance in milk, so for those women who deliver by C section and thus do not transfer lactobacilli from the vagina to the infant, the milk could provide a means for these beneficial organisms to reach the infant gut.
A diverse population of bacteria is present in breast milk dominated by the phyla Proteobacteria and Firmicutes and the taxa Staphylococcus, Pseudomonas, Enterobacteriaceae, Streptococcus and Lactobacillus. While no differences in microbial profiles were apparent based on gestation, mode of delivery or gender, more studies are still needed on what factors do influence bacterial communities in milk and how these changes impact neonatal and maternal health.
Milk collection and processing
A single milk sample (day 6 and onwards post-partum) was collected from 39 Caucasian Canadian women recruited from London, Ontario, and the surrounding area, representing a homogenous community. Even though the samples collected were from different days after birth, with some considered “transitional” and others “mature,” a Kendall’s tau correlation test showed no statistically significant correlation between operational taxonomic unit (OTU) relative abundances over time after a Benjamini-Hochberg correction for multiple hypothesis tests. This shows that it is acceptable to use all milk samples, regardless of when it was collected, in our analyses. Ethical approval was obtained from Western Research Ethics Board and Lawson Health Research Institute, London, Ontario, Canada. Subjects provided written consent for sample collection and subsequent analyses. Participants were excluded if they were suffering from mastitis and were/had been on antibiotics during lactation. Caesarean deliveries were classified as either (i) “non-elective”, if there were complications during labour, or (ii) “elective”, if they were planned in advance or if the health of the foetus and/or mother was at risk prior to labour.
Wearing sterile gloves, the women cleaned their nipple and surrounding area with sterile saline swabs to reduce the presence of skin bacteria. Milk was collected using a sterile HygieniKit Milk Collection System (Ameda, Buffalo Grove, IL, USA) attached to an electric breast pump. Between 5 and 15 ml of milk was pumped into a sterile tube and kept on ice until transfer to the laboratory, which occurred within 1 h of collection. Samples were aliquoted and stored at −20 °C until DNA extraction.
After thawing on ice, 2 ml of milk were spun down at 20,000×g for 10 min and the supernatant was discarded. The pellet was then homogenized in 1.4 ml of ASL buffer (QIAamp® DNA Stool Kit, QIAGEN: Valencia, CA, USA) and 400 mg of 0.1-mm diameter zirconium-glass beads (BioSpec Products, Bartlesville, OK, USA). Mechanical and chemical lyses were performed by bead beading at 4800 rpm for 60 s, then 60 s on ice (repeated twice) using a mini-beadbeater-1 (BioSpec Products) and then incubated at 95 °C for 5 min. Subsequent procedures were performed using the QIAGEN QIAamp® DNA Stool Kit according to the manufacturer’s protocol, with the exception of the last step in which the column was eluted with 120 μl of elution buffer. DNA was stored at −20 °C until further use.
A no template PCR control and a DNA extraction kit reagent control were sequenced alongside the samples. We observed that the taxon abundances in the controls were uncorrelated with the abundances in the experimental samples, and the distance between the controls and samples was large. Thus, we conclude that the controls had different profiles than that of the milk samples (Additional file 7: Figure S4).
V6 16S rRNA gene sequencing
V6-Forward: 5′ACACTCTTTCCCTACACGACGCTCTTCCGATCTnnnn(8)CWACGCGARGAACCTTACC 3′
V6-Reverse: 5′CGGTCTCGGCATTCCTGCTGAAC CGCTCTTCCGATCTnnnn(8)ACRACACGAGCTGACGAC 3′
nnnn indicates four randomly incorporated nucleotides and “8” represents a specific sample barcode sequence. The PCR was carried out in a 42 μl reaction containing 2 μl of DNA template (or nuclease-free water as a negative control), 0.15 μg/μl bovine serum albumin, 20 μl of 2X GoTaq hot start colourless master mix (Promega) and 10 μl of each primer (initial concentration 3.2 pmol/μl). Thermal cycling was carried out in an Eppendorf Mastercycler under the following conditions: initial denaturation at 95 °C for 2 min followed by 25 cycles of 95 °C for 1 min, 55 °C for 1 min, and 72 °C for 1 min. After amplification, the DNA concentration was measured with the Qubit® 2.0 Fluorometer (Invitrogen) using the broad range assay. Equimolar amounts of each PCR product were then pooled together and purified using the QIAquick PCR purification kit (QIAGEN). The pooled PCR purified sample was then paired end sequenced on the Illumina MiSeq platform using a 150-cycle kit with a 2 × 80 run at the London Regional Genomics Centre, London, ON, Canada, following standard operating procedures.
Sequence processing and taxonomic assignment
Custom Perl and Bash scripts were used to de-multiplex the reads and assign barcoded reads to individual samples. Multiple layers of filtering were employed: (i) Paired end sequences were overlapped with Pandaseq, allowing zero mismatch in the overlapped reads; (ii) Reads were kept if the sequence included a perfect match to the V6 16S rRNA gene primers; (iii) Barcodes were 8mers with an edit distance of >4, and reads were kept if the sequence were a perfect match to the barcode; (iv) Reads were clustered by 97 % identity into OTUs using the Uclust algorithm of USEARCH v7  which has a de novo chimera filter built into it; and (v) All singleton OTUs were discarded, and those that represented ≥1 % of the reads in at least one sample were kept.
Taxonomic assignments for each OTU were made by extracting the best hits from the SILVA database  and then manually verified using the Ribosomal Database Project (RDP) SeqMatch tool (rdp.cme.msu.edu) and by BLAST against the Green genes database (greengenes.lbl.gov). Taxonomy was assigned based on hits with the highest percentage identities and coverage. If multiple hits fulfilled this criterion, classification was re-assigned to a higher common taxonomy. A summary of each OTU classification and its sequence is shown in Additional file 8: Table S4. The raw sequencing reads generated in this study have been deposited to the NCBI Short Read Archive (SRA) database accession # SRP 064311.
Weighted and un-weighted UniFrac distances  were calculated in QIIME  by using a tree of OTU sequences built with FASTTREE  based on an OTU sequence alignment made with MUSCLE . The QIIME pipeline was also used to generate PCoA plots to visualize the Bray-Curtis dissimilarity. Changes in the microbial community composition were also analyzed by calculating the generalized UniFrac distance (alpha = 0.5, rooted phylogenetic tree) using the GUniFrac package in R, version 3.1.2 . PERMANOVA was used to test for statistical significance between the groups using 10,000 permutations (QIIME package). Bar plots, box plots, stripcharts, PCoA plots and k means clustering analysis were all generated in R (http://www.R-project.org/).
The ALDEx R package version 2  was used to compare genera, class, family and phyla between preterm and term milk; caesarean and vaginal deliveries, and male and female children. Values reported in this manuscript represent the expected values of 128 Dirichlet Monte Carlo instances. A value of zero indicates that organism abundance is equal to the geometric mean abundance. Thus, organisms more abundant than the mean will have positive values, and those less abundant than the mean will have negative values. Base 2 was used for the logarithm, so differences between values represent fold changes. Significance (p < 0.05) was based on the Benjamini-Hochberg corrected p value of the Welch’s t test and the Wilcoxon test.
Availability of supporting data
The data set supporting the results of this article are included within Additional file 1: Table S1; Additional file 2: Figure S1; Additional file 3: Table S2; Additional file 4: Figure S2; Additional file 5: Figure S3; Additional file 6: Table S3; Additional file 7: Figure S4; Additional file 8: Table S4.
We would like to thank all the women who provided the samples. In addition, we would like to thank Dr. Jean Macklaim and Ms. Ruth Grace Wong for fruitful discussions on 16S rRNA sequencing data and David Carter at the London Regional Genomics Centre. CU was supported by a CIHR Strategic Training Program in Cancer Research and Technology Transfer scholarship. The Translational Breast Cancer Research Studentship Written informed consent was obtained from the subjects for publication of this manuscript and accompanying images.
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.
- Guaraldi F, Salvatori G. Effect of breast and formula feeding on gut microbiota shaping in newborns. Front Cell Infect Microbiol. 2012;2:94.View ArticlePubMedPubMed CentralGoogle Scholar
- Stuebe A. The risks of not breastfeeding for mothers and infants. Rev Obstet Gynecol. 2009;2:222–31.PubMedPubMed CentralGoogle Scholar
- Grönlund M-M, Gueimonde M, Laitinen K, Kociubinski G, Grönroos T, Salminen S, et al. Maternal breast-milk and intestinal bifidobacteria guide the compositional development of the Bifidobacterium microbiota in infants at risk of allergic disease. Clin Exp Allergy. 2007;37:1764–72.View ArticlePubMedGoogle Scholar
- Jost T, Lacroix C, Braegger CP, Chassard C. New insights in gut microbiota establishment in healthy breast fed neonates. PLoS One. 2012;7:e44595.View ArticlePubMedPubMed CentralGoogle Scholar
- Vael C, Vanheirstraeten L, Desager KN, Goossens H. Denaturing gradient gel electrophoresis of neonatal intestinal microbiota in relation to the development of asthma. BMC Microbiol. 2011;11:68.View ArticlePubMedPubMed CentralGoogle Scholar
- Vael C, Nelen V, Verhulst SL, Goossens H, Desager KN. Early intestinal Bacteroides fragilis colonisation and development of asthma. BMC Pulm Med. 2008;8:19.View ArticlePubMedPubMed CentralGoogle Scholar
- Vael C, Verhulst SL, Nelen V, Goossens H, Desager KN. Intestinal microflora and body mass index during the first three years of life: an observational study. Gut Pathog. 2011;3:8.View ArticlePubMedPubMed CentralGoogle Scholar
- Chouraqui JP, Grathwohl D, Labaune JM, Hascoet JM, de Montgolfier I, Leclaire M, et al. Assessment of the safety, tolerance, and protective effect against diarrhea of infant formulas containing mixtures of probiotics or probiotics and prebiotics in a randomized controlled trial. Am J Clin Nutr. 2008;87:1365–73.PubMedGoogle Scholar
- Dilli D, Aydin B, Fettah ND, Özyazıcı E, Beken S, Zenciroğlu A, et al. The ProPre-Save study: effects of probiotics and prebiotics alone or combined on necrotizing enterocolitis in very low birth weight infants. J Pediatr. 2015;166:545–51. e1.View ArticlePubMedGoogle Scholar
- Kankaanpää PE, Yang B, Kallio HP, Isolauri E, Salminen SJ. Influence of probiotic supplemented infant formula on composition of plasma lipids in atopic infants. J Nutr Biochem. 2002;13:364–9.View ArticlePubMedGoogle Scholar
- Mehta R, Petrova A. Biologically active breast milk proteins in association with very preterm delivery and stage of lactation. J Perinatol. 2011;31:58–62.View ArticlePubMedGoogle Scholar
- Dizdar EA, Sari FN, Degirmencioglu H, Canpolat FE, Oguz SS, Uras N, et al. Effect of mode of delivery on macronutrient content of breast milk. J Matern Fetal Neonatal Med. 2014;27:1099–102.View ArticlePubMedGoogle Scholar
- Wang X-F, Cao R-M, Li J, Wu J, Wu S-M, Chen T-X. Identification of sociodemographic and clinical factors associated with the levels of human β-defensin-1 and human β-defensin-2 in the human milk of Han Chinese. Br J Nutr. 2014;111:867–74.View ArticlePubMedGoogle Scholar
- Mehta R, Petrova A. Very preterm gestation and breastmilk cytokine content during the first month of lactation. Breastfeed Med. 2011;6:21–4.View ArticlePubMedGoogle Scholar
- Koenig A, de Albuquerque Diniz EM, Barbosa SFC, Vaz FAC. Immunologic factors in human milk: the effects of gestational age and pasteurization. J Hum Lact. 2005;21:439–43.View ArticlePubMedGoogle Scholar
- Zanardo V, Nicolussi S, Carlo G, Marzari F, Faggian D, Favaro F, et al. Beta endorphin concentrations in human milk. J Pediatr Gastroenterol Nutr. 2001;33:160–4.View ArticlePubMedGoogle Scholar
- Powe CE, Knott CD, Conklin-Brittain N. Infant sex predicts breast milk energy content. Am J Hum Biol. 2010;22:50–4.View ArticlePubMedGoogle Scholar
- Fujita M, Roth E, Lo Y-J, Hurst C, Vollner J, Kendell A. In poor families, mothers’ milk is richer for daughters than sons: a test of Trivers-Willard hypothesis in agropastoral settlements in Northern Kenya. Am J Phys Anthropol. 2012;149:52–9.View ArticlePubMedGoogle Scholar
- Chen J, Bittinger K, Charlson ES, Hoffmann C, Lewis J, Wu GD, et al. Associating microbiome composition with environmental covariates using generalized UniFrac distances. Bioinformatics. 2012;28:2106–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Hunt KM, Foster JA, Forney LJ, Schütte UME, Beck DL, Abdo Z, et al. Characterization of the diversity and temporal stability of bacterial communities in human milk. PLoS One. 2011;6:e21313.View ArticlePubMedPubMed CentralGoogle Scholar
- Jost T, Lacroix C, Braegger C, Chassard C. Assessment of bacterial diversity in breast milk using culture-dependent and culture-independent approaches. Br J Nutr. 2013;110:1253–62.View ArticlePubMedGoogle Scholar
- Urbaniak C, Burton JP, Reid G. Breast, milk and microbes: a complex relationship that does not end with lactation. Womens Health (Lond Engl). 2012;8:385–98.View ArticleGoogle Scholar
- Cabrera-Rubio R, Collado MC, Laitinen K, Salminen S, Isolauri E, Mira A. The human milk microbiome changes over lactation and is shaped by maternal weight and mode of delivery. Am J Clin Nutr. 2012;96:544–51.View ArticlePubMedGoogle Scholar
- Schwarzberg K, Le R, Bharti B, Lindsay S, Casaburi G, Salvatore F, et al. The personal human oral microbiome obscures the effects of treatment on periodontal disease. PLoS One. 2014;9:e86708.View ArticlePubMedPubMed CentralGoogle Scholar
- Lundequist B, Nord CE, Winberg J. The composition of the faecal microflora in breastfed and bottle fed infants from birth to eight weeks. Acta Paediatr Scand. 1985;74:45–51.View ArticlePubMedGoogle Scholar
- Balmer SE, Wharton BA. Diet and faecal flora in the newborn: breast milk and infant formula. Arch Dis Child. 1989;64:1672–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Adlerberth I. Factors influencing the establishment of the intestinal microbiota in infancy. Nestle Nutr Workshop Ser Pediatr Program. 2008;62:13–29. Discussion 29–33.View ArticlePubMedGoogle 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 U S A. 2011;108 Suppl:4578–85.View ArticleGoogle Scholar
- Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.View ArticleGoogle Scholar
- Hapfelmeier S, Lawson MAE, Slack E, Kirundi JK, Stoel M, Heikenwalder M, et al. Reversible microbial colonization of germ-free mice reveals the dynamics of IgA immune responses. Science. 2010;328:1705–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Gan XT, Ettinger G, Huang CX, Burton JP, Haist JV, Rajapurohitam V, et al. Probiotic administration attenuates myocardial hypertrophy and heart failure after myocardial infarction in the rat. Circ Heart Fail. 2014;7:491–9.View ArticlePubMedGoogle Scholar
- Tunzi CR, Harper PA, Bar-Oz B, Valore EV, Semple JL, Watson-MacDonell J, et al. Beta-defensin expression in human mammary gland epithelia. Pediatr Res. 2000;48:30–5.View ArticlePubMedGoogle Scholar
- Heikkila MP, Saris PEJ. Inhibition of Staphylococcus aureus by the commensal bacteria of human milk. J Appl Microbiol. 2003;95:471–8.View ArticlePubMedGoogle Scholar
- Shu M, Wang Y, Yu J, Kuo S, Coda A, Jiang Y, et al. Fermentation of Propionibacterium acnes, a commensal bacterium in the human skin microbiome, as skin probiotics against methicillin-resistant Staphylococcus aureus. PLoS One. 2013;8:e55380.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticlePubMedGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35:7188–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35.View ArticlePubMedPubMed CentralGoogle 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 ArticlePubMedPubMed CentralGoogle Scholar
- Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26:1641–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.View ArticlePubMedPubMed CentralGoogle Scholar
- Fernandes AD, Reid JN, Macklaim JM, McMurrough TA, Edgell DR, Gloor GB. Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome. 2014;2:15.View ArticlePubMedPubMed CentralGoogle Scholar
- Aitchson J, Greenacre M. Biplots of compositional data. J R Stat Soc Ser C. 2002;51:375–92.View ArticleGoogle Scholar