Effect of chemotherapy on the microbiota and metabolome of human milk, a case report
© Urbaniak et al.; licensee BioMed Central Ltd. 2014
Received: 4 March 2014
Accepted: 28 May 2014
Published: 11 July 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.
Keywords16S rRNA gene sequencing Human milk microbiome Metabolome
false discovery rate
gas chromatography-mass spectrometry
operational taxonomic unit
principal component analysis
polymerase chain reaction
principal coordinate analysis
Ribosomal Database Project.
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.
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.
Metabolites significantly altered by chemotherapy
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.
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.
- 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.PubMedView Article
- 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.PubMedView Article
- Stuebe A: The risks of not breastfeeding for mothers and infants. Rev Obstet Gynecol 2009, 2:222–231.PubMed CentralPubMed
- 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.PubMedView Article
- Guaraldi F, Salvatori G: Effect of breast and formula feeding on gut microbiota shaping in newborns. Front Cell Infect Microbiol 2012, 2:94.PubMed CentralPubMed
- 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.PubMed
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- German JB, Dillard CJ, Ward RE: Bioactive components in milk. Curr Opin Clin Nutr Metab Care 2002, 5:653–658.PubMedView Article
- Donaldson M, Goodchild JH: Pregnancy, breast-feeding and drugs used in dentistry. J Am Dent Assoc 2012, 143:858–871.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- Pistilli B, Bellettini G, Giovannetti E, Codacci-Pisanelli G, Azim HA Jr, 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.PubMedView Article
- Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26:2460–2461.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 2004, 5:113.PubMed CentralPubMedView Article
- Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 2005, 71:8228–8235.PubMed CentralPubMedView Article
- Hill TC, Walsh KA, Harris JA, Moffett BF: Using ecological diversity measures with bacterial communities. FEMS Microbiol Ecol 2003, 43:1–11.PubMedView Article
- R Development Core Team: R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2006.
- 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.PubMed CentralPubMedView Article
- Jaynes ET: Probability Theory: the Logic of Science. Cambridge, UK: Cambridge University Press; 2003.View Article
- Aitchison J, Egozcue J: Compositional data analysis: Where are we and where should we be heading? Math Geol 2005, 37:829–850.View Article
- Stein S: An integrated method for spectrum extraction and compound identification from GC/MS data. J Am Soc Mass Spectrom 1999, 10:770–781.View Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- Jost T, Lacroix C, Braegger CP, Chassard C: New insights in gut microbiota establishment in healthy breast fed neonates. PLoS One 2012, 7:e44595.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- Ward TL, Hosid S, Ioshikhes I, Altosaar I: Human milk metagenome: a functional capacity analysis. BMC Microbiol 2013, 13:116.PubMed CentralPubMedView Article
- 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.PubMedView Article
- Heikkila MP, Saris PE: Inhibition of Staphylococcus aureus by the commensal bacteria of human milk. J Appl Microbiol 2003, 95:471–478.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- Lin HV, Frassetto A, Kowalik EJ Jr, 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.PubMed CentralPubMedView Article
- Ramsey MM, Rumbaugh KP, Whiteley M: Metabolite cross-feeding enhances virulence in a model polymicrobial infection. PLoS Pathog 2011, 7:e1002012.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.