- Short report
- Open Access
Effect of chemotherapy on the microbiota and metabolome of human milk, a case report
Microbiomevolume 2, Article number: 24 (2014)
Human milk is an important source of bacteria for the developing infant and has been shown to influence the bacterial composition of the neonatal gut, which in turn can affect disease risk later in life. Human milk is also an important source of nutrients, influencing bacterial composition but also directly affecting the host. While recent studies have emphasized the adverse effects of antibiotic therapy on the infant microbiota, the effects of maternal chemotherapy have not been previously studied. Here we report the effects of drug administration on the microbiota and metabolome of human milk.
Mature milk was collected every two weeks over a four month period from a lactating woman undergoing chemotherapy for Hodgkin’s lymphoma. Mature milk was also collected from healthy lactating women for comparison. Microbial profiles were analyzed by 16S sequencing and the metabolome by gas chromatography–mass spectrometry.
Chemotherapy caused a significant deviation from a healthy microbial and metabolomic profile, with depletion of genera Bifidobacterium, Eubacterium, Staphylococcus and Cloacibacterium in favor of Acinetobacter, Xanthomonadaceae and Stenotrophomonas. The metabolites docosahexaenoic acid and inositol known for their beneficial effects were also decreased.
With milk contents being critical for shaping infant immunity and development, consideration needs to be given to the impact of drugs administered to the mother and the long-term potential consequences for the health of the infant.
Colonization of the neonatal gut plays a pivotal role in gastrointestinal, immunological and neurological development, with one of the initial sources of bacteria being the mother’s milk [1, 2]. Breastfed infants have been shown to have lower incidences of asthma, diarrhea, and necrotizing enterocolitis compared with formula-fed infants . This protective effect may be due, in part, to the types of bacteria present in milk, as infants fed formula supplemented with probiotics were better protected against these conditions compared to those just fed formula [4–6]. The bacteria acquired during infancy can influence disease risk later in life and play a major role in determining the future composition of the adult microbiome . Thus, factors that affect the milk microbiota have important health consequences for the child not only during development but also into adulthood. In addition to the microbiota, the metabolites of human milk, such as fatty acids, carbohydrates, proteins and vitamins, also play an important role in infant development and long-term health [7–10].
Post-delivery, many women are prescribed pharmaceutical agents for various reasons. While most over-the-counter drugs and antibiotics are not contraindicated during breast feeding [11, 12], when it comes to chemotherapeutics, the recommendation is that breastfeeding should be avoided until the drug has been cleared from the milk . In a case report of a 70 mg infusion of cisplatin, no detectable levels were found in milk after 66 hours , and in another case study using doxorubicin (trade name Adriamycin), no detectable levels were seen after 72 hours . In our particular study, the subject was advised that breastfeeding could resume 12 days after each chemotherapy session.Here we present the first report on the effects of chemotherapy on microbial and metabolomic profiles in human milk over a 4-month period in a breastfeeding woman undergoing treatment for Hodgkin’s lymphoma (Figure 1).
Clinical samples and study design
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.
Milk collection and processing
Mature milk was collected from a lactating woman undergoing the ABVD chemotherapy regime (Adriamycin (40 mg), Bleomycin (16 units), Vinblastine (9.6 mg), Dacarbazine (600 mg) for Hodgkin’s lymphoma at the London Health Sciences Center, London, Ontario. Mature milk was also collected from 8 healthy women recruited from London, Ontario and the surrounding area. 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 hour 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 13,000 g for 10 minutes and the supernatant 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 4,800 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 minutes. 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.
V6 16S rRNA gene sequencing
The genomic DNA isolated from the clinical samples was amplified using the barcoded primers V6-LT: 5′CCATCTCATCCCTGCGTGTCTCCGACTCAGNNNNNCWACGCGARGAACCTTACC3′ and V6-RT: 5′ CCTCTCTATGGGCAGTCGGTGATACRACACGAGCTGACGAC3′, which amplify the V6 hypervariable region of the 16S rRNA gene. The PCR was carried out in a 40 μl reaction containing 5 μl of DNA template (or nuclease-free water as a negative control), 1.5 mM MgCl2, 0.8 μM of each primer, 4 μl of 10× PCR Buffer (Invitrogen, Burlington, ON, Canada), 0.2 mM dNTPs, 0.05U Taq Polymerase (Invitrogen) and 0.15 μg/μl bovine serum albumin. Thermal cycling was carried out in an Eppendorf Mastercyler under the following conditions: initial denaturation at 95°C for 2 minutes followed by 25 cycles of 95°C for 1 minute, 55°C for 1 minute and 72°C for 1 minute. 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 pooled together and purified using the QIAquick PCR purification kit (QIAGEN). The PCR purified sample was then sent to the London Regional Genomics Center, London, Ontario, Canada for V6 16S rRNA gene sequencing using the Ion Torrent platform as per the center’s standard operating procedure.
Sequence processing and taxonomic assignment
Custom Perl and Bash scripts were used to de-multiplex the reads and assign barcoded reads to individual samples. Reads were kept if the sequence included a perfect match to the barcode and the V6 16S rRNA gene primers. Reads were clustered by 97% identity into operational taxonomic units (OTUs) using UClust v.3.0.617 . OTUs that represented ≥1% of the reads in at least one sample were kept, while those that did not meet the cutoff were discarded. Taxonomic assignments for each OTU were made by the Ribosomal Database Project (RDP) SeqMatch tool . From the top 20 matches to the RDP named isolates database, the full taxonomy was retained for matches with the highest S_ab score. For multiple top matches with equal scores, the highest common taxonomy was retained (for example, genus level if multiple species matched equally well). Since the maximum number of matches displayed per sequence is 20, the RDP taxonomic assignments were verified by BLAST against the Greengenes named isolates database with an output of 100 hits . 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. In instances where the highest percentage identity/coverage yielded a single match, if this were <90% and the S_ab score from RDP was <0.7, taxonomy was assigned at the family level instead of at the genus level. A summary of each OTU classification and its sequence is shown in Additional file 1. The raw sequencing reads generated in this study have been deposited to the NCBI Short Read Archive (SRA) database [SRA:SRP041626].
Weighted UniFrac distances were calculated in QIIME  by using a phylogenetic tree of OTU sequences built with FastTree  and based on an OTU sequence alignment with MUSCLE . The QIIME pipeline was also used to calculate Shannon’s diversity index (logarithms with base 2) and to generate principal coordinate analysis (PCoA) plots. Weighted UniFrac distances compare microbial profiles (presence/absence and abundance) between samples (i.e., beta-diversity)  while Shannon’s diversity index evaluates the microbial diversity within a sample (i.e., alpha diversity). The higher the Shannon’s diversity index, the more diverse a sample is and a value of zero indicates the presence of only one species . PCoA plots allow one to visualize the UniFrac distance matrix and plot the values on a set of orthogonal axes that capture the greatest amount of variation between all samples tested. For beta-diversity analyses, the data set was rarified to the lowest read count/sample, which was 734 reads. A summary of clinical data, including total number of sequence reads per sample, is shown in Additional file 2. Barplots, boxplots and stripcharts were all generated in R .
The ALDEx R package version 2  was used to compare genera between the non-treatment and chemotherapy treatment groups (as portrayed in the boxplots). Microbiome data represent proportional distributions and are thus not independent of each other. This means that a decrease in one organism will inevitably lead to a concomitant increase in another organism. For example, if a sample has two organisms A (50%) and B (50%) and A is completely killed by an antibiotic, the proportion of B in that sample will now be 100% even if its actual abundance has not changed. The ALDEx R package estimates the technical variation inherent in high-throughput sequencing by Monte-Carlo sampling from a Dirichlet distribution . The Monte-Carlo replicates are transformed using the centered log-ratio transformation, which takes the logarithm of the Monte-Carlo estimates of organism abundances in each sample divided by the per-sample geometric mean organism abundance . This transformation has several desirable properties that do not exist in proportional data, notably subcomposition coherence and linear sample independence. Data transformed in this way permit the use of standard statistical tests to determine significance. 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. Statistical significance for these comparisons was determined by a Mann-Whitney U test with P < 0.05 and a false discovery rate (FDR) of < 0.1 using the q values output by the fdrtool R package .
The unpaired Student’s t-test was used to compare Shannon’s diversity index (P < 0.05).
Sample preparation gas chromatography-mass spectrometry
To extract metabolites, 100 μl of milk were mixed with 400 μl pure methanol. Samples were vortexed for 15 s and centrifuged for 10 minutes at 9,000 g. Supernatants (200 μl) were transferred to gas chromatography-mass spectrometry (GC-MS) vials and 2.5 μl of ribitol solution (2 mg/ml) was added to each vial as an internal standard. Samples were dried to completeness using a SpeedVac. After drying, 100 μl of 2% methoxyamine•HCl in pyridine (MOX) was added to each sample for derivitization and samples were incubated at 50°C for 90 minutes. N-methyl-N-(trimethylsilyl) trifluoroacetamide (100 μl) was then added to each vial and incubated at 50°C for 30 minutes. Samples were transferred to microinserts before running on GC-MS (Agilent 7890A GC, 5975 inert MSD with triple axis detector, 30 m DB5-MS column with 10 m duraguard column). Samples were run using 1 μl injections on scan mode and a solvent delay of 10 minutes. Run time was 60 minutes per sample. Each sample was run twice non-consecutively to ensure consistency throughout the sequence.
GC-MS data analysis
Chromatogram files were converted to ELU format using the AMDIS Mass Spectrometry software . Chromatograms were aligned and abundance of metabolites calculated using the Spectconnect software  with the support threshold set to low. In order to determine if differences between week 0 and chemotherapy (weeks 2 to 16) existed, principle component analysis (PCA) was conducted in SIMCA (Umetrics, San Jose, CA, USA) using the relative abundance matrix output from Spectconnect. Data were mean centered and pareto scaled prior to PCA. Independent unpaired t-tests with Bonferroni correction were calculated in Excel to determine metabolites that were significantly altered by chemotherapy (P < 0.05). Compounds that also contributed to the separation of week 0 from chemotherapy samples according to the PCA loadings plot (compounds in bottom left quadrant) were chosen for further investigation.
Chemotherapy affected both bacterial diversity and bacterial profiles in human milk. Bacterial diversity within samples was lower in milk collected throughout chemotherapy compared with milk samples collected at week 0 and from healthy lactating women (Figure 2). Bacterial profiles at week 0 were similar to those from healthy women, although this changed within 2 weeks of treatment (Figure 3A). Samples collected at weeks 4 to 16 shared similar profiles and differed from week 2 and from week 0/healthy samples (Figure 3A). These differences were not due to natural changes over time, as the bacterial community in two milk samples analyzed from a healthy subject did not change over a 4-month period (Figure 3A, green samples). The bar plot in Figure 3B shows the bacterial communities present in these samples with a striking increase in abundances of Acinetobacter and Xanthomonadaceae in milk collected during chemotherapy. A comparison of relative abundances of Acinetobacter, Xanthomonadaceae and Stenotrophomonas (a genus belonging to the Xanthomonadaceae family) between the chemotherapy (weeks 4 to 16) and non-treatment (week 0 and healthy samples) groups is displayed in Figure 4 and were significantly higher during chemotherapy. We also examined the differences between three bacteria believed to confer beneficial health effects to the infant, Bifidobacterium, Eubacterium and Lactobacillus. The first two were significantly decreased during chemotherapy whereas no differences were observed for Lactobacillus (Figure 4). Overall, a total of 22 out of the 49 genera identified were differentially abundant between the two groups (Additional file 3). While the core microbiome (that is, taxa that were present in 100% of the samples) was somewhat similar between the two groups, it is interesting to note that Stenotrophomonas was present in every chemotherapy sample and Lactobacillus and Eubacterium were present in every healthy and week 0 sample (Additional file 4).
The metabolic profile also changed as a result of chemotherapy (Figure 5) and was similar between weeks 2 and 16, but different to that observed at week 0. A total of 226 metabolites were detected by our GC-MS method, 12 of which were significantly different between the week 0 and chemotherapy (weeks 2 to 16) groups (Table 1). Additional file 5 shows the relative abundances of all metabolites detected in milk and stripcharts in Additional file 6 show the distribution of key metabolites detected between the two groups.
This study shows that a course of chemotherapy has significant effects on bacterial and metabolic profiles in human milk, moving away from those of healthy lactating women. Of note, the subject did not report any additional drug use, antibiotics, illness or major changes in diet over the course of the study.
The consequences of decreased bacterial diversity in human milk and the implications on the child are still unknown; however, the decreased milk diversity could impact intestinal diversity and it has been shown that low intestinal diversity in the first weeks of life is associated with necrotizing enterocolitis  and an increased risk of allergy and atopy in school-age children . Lower intestinal diversity has also been observed in children with type 1 diabetes compared to age-matched controls .
In addition to overall changes in microbial profiles, we observed a significant decrease in the relative abundance of Bifidobacterium in the chemotherapy group compared to the non-treatment group. Bifidobacterium is the predominant organism in the gut of breastfed infants, attributed to its ability to metabolize the human milk oligosaccharides present in large amounts in milk [33, 34]. Maternal levels can also impact gut Bifidobacterium abundance, with low levels in milk correlating with low levels in the neonatal gut . The potential consequences of decreased numbers of Bifidobacterium being passed on from the mother to the neonate could be an increased risk of asthma or obesity later in life. High levels of Bacteroides have been reported in the gut of infants with low levels of Bifidobacterium and early colonization with high Bacteroides counts has been associated with an increased risk of developing asthma and obesity [37–39]. In addition, depleted levels of Bifidobacterium have been shown to promote colonization of opportunistic pathogens such as Klebsiella and Citrobacter.
Staphylococcus, including coagulase negative species, are one of the predominant organisms in human milk [40–42], and were also significantly reduced as a result of chemotherapy. It has been shown that numerous human milk isolates of Staphylococcus epidermidis can inhibit the growth of Staphylococcus aureus, the main causative agent of mastitis, which is a painful inflammatory condition of the breast that often leads to premature cessation of breastfeeding in many women. While we were not able to identify the Staphylococcus in our samples down to the species level with 16S rRNA gene sequencing, culture analysis on mannitol salt agar plates did show that the Staphylococcus isolates were not S. aureus and the select few that were tested were coagulase negative. This reduction of Staphylococcus (likely coagulase negative species) as a result of chemotherapy could make lactating women more prone to infections, affecting both themselves and their infants. Like Bifidobacterium, Staphylococcus is passed from the milk to the neonate, with higher numbers in the intestine of breastfed compared to formula-fed infants . Interestingly, a metagenome analysis revealed the presence of immunosuppressive motifs in bacterial DNA from human milk, with the majority of these belonging to Staphylococcus. The exposure of the neonate to this DNA, either ingested from the milk or through live bacteria that have released their DNA once in the gut, could help to regulate the infant’s immune response against a variety of innocuous bacterial, environmental and food antigens.
The utilization of bacterial products by other bacteria is termed metabolic cross-feeding and plays an important role in bacterial selection. For example, the byproducts of bacterial metabolism, such as lactate and acetate production, are utilized as an energy source by many butyrate-producing bacteria [44–47] such as Eubacterium, which was decreased in our chemotherapy group. Butyrate is important for health, as it reduces inflammation and metabolic diseases, promotes colonic repair and protects against colon cancer [48, 49]. On the other hand, some pathogens persist and cause disease only in the presence of certain commensal bacteria , likely due to the metabolites produced. Thus, changes in bacterial communities in human milk will inevitably alter the metabolic milieu, selecting for bacteria able to utilize those metabolites. As a result, potential shifts from a healthy and balanced intestinal microbiota can occur, having important consequences on health.
Docosahexaenoic acid (DHA), inositol and an unknown polyunsaturated fatty acid were among the 12 metabolites that differed between week 0 and weeks 2 to 16, with reduced levels detected during chemotherapy (weeks 2 to 16). DHA is the most abundant long chain polyunsaturated fatty acid in the brain, retina and nerve cells and is supplied mainly through breast milk . DHA deficiencies lead to reduced brain, eye and neuronal development  and it has been observed that breastfed infants have better visual acuity and neuronal development compared to those fed formula . Reduced levels of DHA and alpha-linolenic acid (a precursor of DHA) have been reported in milk of mothers with atopic children compared to milk from mothers with non-atopic children [53, 54]. Another principal metabolite in the neonatal brain is inositol, which is important for osmoregulation, cellular nutrition and detoxification . Palmitic acid levels were also reduced during chemotherapy, though the results were not significant. Palmitic acid is the most abundant lipid in human milk and has been shown to increase bone strength in infants  and limit intestinal damage and pro-inflammatory immune responses in mice . While changes in metabolite concentrations do occur over the course of lactation, especially between colostrum and mature milk, no changes in the above metabolites over the course of the first year of life have been observed in mature milk [58–61]. Due to the high variability in milk metabolites between individuals, we did not have enough samples from healthy women to make substantive claims as to how the Hodgkin’s patient compared, although there were no obvious differences in the metabolic profiles of control samples taken at early compared to later stages of lactation.
We recognize that the main limitation of the study is its single case study content of a patient undergoing chemotherapy. However, the findings were revealing. Of note, even with the disease, her milk microbiota before treatment was similar to that of healthy lactating women and only after intervention did the microbiota patterns alter. It has been reported by Hunt et al. that the milk microbiota of a specific individual is stable over time, consistent with our study from the healthy sample collected 4 months apart. We believe that the changes in the microbiota are a result of therapeutic agents and not specific to just one patient.
Bacterial and metabolic compositions in human milk, so critical for immunity and infant development, can change significantly after maternal exposure to chemotherapeutic agents. Further larger cohort studies are warranted to examine microbiota and metabolomic changes associated with chemotherapy and other medications prescribed to lactating mothers and the consequences for the microbiome, the metabolome and long-term health of infants.
false discovery rate
gas chromatography-mass spectrometry
operational taxonomic unit
principal component analysis
polymerase chain reaction
principal coordinate analysis
Ribosomal Database Project.
Martin V, Maldonado-Barragan A, Moles L, Rodriguez-Banos M, Campo RD, Fernandez L, Rodriguez JM, Jimenez E: Sharing of bacterial strains between breast milk and infant feces. J Hum Lact. 2012, 28: 36-44.
Martin R, Langa S, Reviriego C, Jiminez E, Marin ML, Xaus J, Fernandez L, Rodriguez JM: Human milk is a source of lactic acid bacteria for the infant gut. J Pediatr. 2003, 143: 754-758.
Stuebe A: The risks of not breastfeeding for mothers and infants. Rev Obstet Gynecol. 2009, 2: 222-231.
Silvers KM, Frampton CM, Wickens K, Pattemore PK, Ingham T, Fishwick D, Crane J, Town GI, Epton MJ, New Zealand Asthma and Allergy Cohort Study Group: Breastfeeding protects against current asthma up to 6 years of age. J Pediatr. 2012, 160: 991–6.e1-
Guaraldi F, Salvatori G: Effect of breast and formula feeding on gut microbiota shaping in newborns. Front Cell Infect Microbiol. 2012, 2: 94-
Chouraqui JP, Grathwohl D, Labaune JM, Hascoet JM, de Montgolfier I, Leclaire M, Giarre M, Steenhout P: 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-1373.
Boudry G, David ES, Douard V, Monteiro IM, Le Huerou-Luron I, Ferraris RP: Role of intestinal transporters in neonatal nutrition: carbohydrates, proteins, lipids, minerals, and vitamins. J Pediatr Gastroenterol Nutr. 2010, 51: 380-401.
Kunz C, Rudloff S, Baier W, Klein N, Strobel S: Oligosaccharides in human milk: structural, functional, and metabolic aspects. Annu Rev Nutr. 2000, 20: 699-722.
Agostoni C: Role of long-chain polyunsaturated fatty acids in the first year of life. J Pediatr Gastroenterol Nutr. 2008, 47 (Suppl 2): S41-S44.
German JB, Dillard CJ, Ward RE: Bioactive components in milk. Curr Opin Clin Nutr Metab Care. 2002, 5: 653-658.
Donaldson M, Goodchild JH: Pregnancy, breast-feeding and drugs used in dentistry. J Am Dent Assoc. 2012, 143: 858-871.
Drugs and Lactation Database. [http://toxnet.nlm.nih.gov/cgi-bin/sis/htmlgen?LACT]
Gartner LM, Morton J, Lawrence RA, Naylor AJ, O’Hare D, Schanler RJ, Eidelman AI, American Academy of Pediatrics Section on Breastfeeding: Breastfeeding and the use of human milk. Pediatrics. 2005, 115: 496-506.
Hays KE, Ryu RJ, Swisher EM, Reed E, McManus T, Rybeck B, Petros WP, Hebert MF: Duration of cisplatin excretion in breast milk. J Hum Lact. 2013, 29: 469-472.
Pistilli B, Bellettini G, Giovannetti E, Codacci-Pisanelli G, Azim HA, Benedetti G, Sarno MA, Peccatori FA: Chemotherapy, targeted agents, antiemetics and growth-factors in human milk: how should we counsel cancer patients about breastfeeding?. Cancer Treat Rev. 2013, 39: 207-211.
Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461.
Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM: The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37 (Database issue): D141-D145.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006, 72: 5069-5072.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336.
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-1650.
Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 113-
Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005, 71: 8228-8235.
Hill TC, Walsh KA, Harris JA, Moffett BF: Using ecological diversity measures with bacterial communities. FEMS Microbiol Ecol. 2003, 43: 1-11.
R Development Core Team: R: A Language and Environment for Statistical Computing. 2006, Vienna, Austria: R Foundation for Statistical Computing
Fernandes AD, Reid J, 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-
Jaynes ET: Probability Theory: the Logic of Science. 2003, Cambridge, UK: Cambridge University Press
Aitchison J, Egozcue J: Compositional data analysis: Where are we and where should we be heading?. Math Geol. 2005, 37: 829-850.
Stein S: An integrated method for spectrum extraction and compound identification from GC/MS data. J Am Soc Mass Spectrom. 1999, 10: 770-781.
Styczynski MP, Moxley JF, Tong LV, Walther JL, Jensen KL, Stephanopoulos GN: Systematic identification of conserved metabolites in GC/MS data for metabolomics and biomarker discovery. Anal Chem. 2007, 79: 966-973.
Mshvildadze M, Neu J, Shuster J, Theriaque D, Li N, Mai V: Intestinal microbial ecology in premature infants assessed with non-culture-based techniques. J Pediatr. 2010, 156: 20-25.
Madan JC, Farzan SF, Hibberd PL, Karagas MR: Normal neonatal microbiome variation in relation to environmental factors, infection and allergy. Curr Opin Pediatr. 2012, 24: 753-759.
Giongo A, Gano KA, Crabb DB, Mukherjee N, Novelo LL, Casella G, Drew JC, Ilonen J, Knip M, Hyoty H, Veijola R, Simell T, Simell O, Neu J, Wasserfall CH, Schatz D, Atkinson MA, Triplett EW: Toward defining the autoimmune microbiome for type 1 diabetes. ISME J. 2011, 5: 82-91.
Sela DA, Chapman J, Adeuya A, Kim JH, Chen F, Whitehead TR, Lapidus A, Rokhsar DS, Lebrilla CB, German JB, Price NP, Richardson PM, Mills DA: The genome sequence of Bifidobacterium longum subsp. infantis reveals adaptations for milk utilization within the infant microbiome. Proc Natl Acad Sci U S A. 2008, 105: 18964-18969.
Asakuma S, Hatakeyama E, Urashima T, Yoshida E, Katayama T, Yamamoto K, Kumagai H, Ashida H, Hirose J, Kitaoka M: Physiology of consumption of human milk oligosaccharides by infant gut-associated bifidobacteria. J Biol Chem. 2011, 286: 34583-34592.
Gronlund MM, Gueimonde M, Laitinen K, Kociubinski G, Gronroos T, Salminen S, Isolauri E: 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-1772.
Jost T, Lacroix C, Braegger CP, Chassard C: New insights in gut microbiota establishment in healthy breast fed neonates. PLoS One. 2012, 7: e44595-
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-
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-
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-
Hunt KM, Foster JA, Forney LJ, Schutte UM, Beck DL, Abdo Z, Fox LK, Williams JE, McGuire MK, McGuire MA: Characterization of the diversity and temporal stability of bacterial communities in human milk. PLoS One. 2011, 6: e21313-
Ward TL, Hosid S, Ioshikhes I, Altosaar I: Human milk metagenome: a functional capacity analysis. BMC Microbiol. 2013, 13: 116-
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-1262.
Heikkila MP, Saris PE: Inhibition of Staphylococcus aureus by the commensal bacteria of human milk. J Appl Microbiol. 2003, 95: 471-478.
Belenguer A, Duncan SH, Calder AG, Holtrop G, Louis P, Lobley GE, Flint HJ: Two routes of metabolic cross-feeding between Bifidobacterium adolescentis and butyrate-producing anaerobes from the human gut. Appl Environ Microbiol. 2006, 72: 3593-3599.
Barcenilla A, Pryde SE, Martin JC, Duncan SH, Stewart CS, Henderson C, Flint HJ: Phylogenetic relationships of butyrate-producing bacteria from the human gut. Appl Environ Microbiol. 2000, 66: 1654-1661.
Duncan SH, Holtrop G, Lobley GE, Calder AG, Stewart CS, Flint HJ: Contribution of acetate to butyrate formation by human faecal bacteria. Br J Nutr. 2004, 91: 915-923.
Duncan SH, Louis P, Flint HJ: Lactate-utilizing bacteria, isolated from human feces, that produce butyrate as a major fermentation product. Appl Environ Microbiol. 2004, 70: 5810-5817.
Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost FJ, Brummer RJ: Review article: the role of butyrate on colonic function. Aliment Pharmacol Ther. 2008, 27: 104-119.
Lin HV, Frassetto A, Kowalik EJ, Nawrocki AR, Lu MM, Kosinski JR, Hubert JA, Szeto D, Yao X, Forrest G, Marsh DJ: Butyrate and propionate protect against diet-induced obesity and regulate gut hormones via free fatty acid receptor 3-independent mechanisms. PLoS One. 2012, 7: e35240-
Ramsey MM, Rumbaugh KP, Whiteley M: Metabolite cross-feeding enhances virulence in a model polymicrobial infection. PLoS Pathog. 2011, 7: e1002012-
Morse NL: Benefits of docosahexaenoic acid, folic acid, vitamin D and iodine on foetal and infant brain development and function following maternal supplementation during pregnancy and lactation. Nutrients. 2012, 4: 799-840.
Hoffman DR, Boettcher JA, Diersen-Schade DA: Toward optimizing vision and cognition in term infants by dietary docosahexaenoic and arachidonic acid supplementation: a review of randomized controlled trials. Prostaglandins Leukot Essent Fatty Acids. 2009, 81: 151-158.
Yu G, Duchen K, Bjorksten B: Fatty acid composition in colostrum and mature milk from non-atopic and atopic mothers during the first 6 months of lactation. Acta Paediatr. 1998, 87: 729-736.
Duchen K, Casas R, Fageras-Bottcher M, Yu G, Bjorksten B: Human milk polyunsaturated long-chain fatty acids and secretory immunoglobulin A antibodies and early childhood allergy. Pediatr Allergy Immunol. 2000, 11: 29-39.
Kok RD, van den Bergh AJ, Heerschap A, Nijland R, van den Berg PP: Metabolic information from the human fetal brain obtained with proton magnetic resonance spectroscopy. Am J Obstet Gynecol. 2001, 185: 1011-1015.
Litmanovitz I, Davidson K, Eliakim A, Regev RH, Dolfin T, Arnon S, Bar-Yoseph F, Goren A, Lifshitz Y, Nemet D: High Beta-palmitate formula and bone strength in term infants: a randomized, double-blind, controlled trial. Calcif Tissue Int. 2013, 92: 35-41.
Lu P, Bar-Yoseph F, Levi L, Lifshitz Y, Witte-Bouma J, de Bruijn AC, Korteland-van Male AM, van Goudoever JB, Renes IB: High beta-palmitate fat controls the intestinal inflammatory response and limits intestinal damage in mucin Muc2 deficient mice. PLoS One. 2013, 8: e65878-
Marangoni F, Agostoni C, Lammardo AM, Giovannini M, Galli C, Riva E: Polyunsaturated fatty acid concentrations in human hindmilk are stable throughout 12-months of lactation and provide a sustained intake to the infant during exclusive breastfeeding: an Italian study. Br J Nutr. 2000, 84: 103-109.
Jozwik M, Jozwik M, Teng C, Jozwik M, Battaglia FC: Human breast milk sugars and polyols over the first 10 puerperium days. Am J Hum Biol. 2013, 25: 198-204.
Ribeiro M, Balcao V, Guimaraes H, Rocha G, Moutinho C, Matos C, Almeida C, Casal S, Guerra A: Fatty acid profile of human milk of Portuguese lactating women: prospective study from the 1st to the 16th week of lactation. Ann Nutr Metab. 2008, 53: 50-56.
Lopez-Lopez A, Lopez-Sabater MC, Campoy-Folgoso C, Rivero-Urgell M, Castellote-Bargallo AI: Fatty acid and sn-2 fatty acid composition in human milk from Granada (Spain) and in infant formulas. Eur J Clin Nutr. 2002, 56: 1242-1254.
We would like to thank Jean Macklaim for her advice and guidance in analyzing the 16S rRNA sequencing data and Tim McDowell at Agriculture and Agri-Food Canada for technical support with GC-MS. Funding sources to individuals: CU - CIHR Strategic Training Program in Cancer Research and Technology Transfer; and the Translational Breast Cancer Studentship; AM - Ontario Graduate Scholarship. Written informed consent was obtained from the patients for publication of this manuscript and accompanying images.
The authors declare that they have no competing interests.
CU designed the study, recruited subjects, collected and prepared the samples for microbial analysis, analyzed 16S rRNA sequencing data and wrote the manuscript. AM performed GC-MS analysis on the milk samples and reviewed the manuscript. MA helped with study design and recruitment and reviewed the manuscript. JB provided input into study design and reviewed the manuscript. MS contributed substantially to the design and acquisition of the metabolomic data, analysis and manuscript review. GG provided co-supervision, instruction and input into microbiome acquisition and analysis and manuscript review. GR conceptualized the study, helped with study design and manuscript writing, supervised data collection and analysis and provided financial support. All authors read and approved the final manuscript.