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
Skip to main content
© 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.
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.
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.
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.
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.
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).
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.
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.
false discovery rate
gas chromatography-mass spectrometry
operational taxonomic unit
principal component analysis
polymerase chain reaction
principal coordinate analysis
Ribosomal Database Project.
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.
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.