BACKGROUND
the benefits of human milk compared to the use of commercial infant formulas are largely realized because of its bioactive components, including prebiotics, immune proteins and the microbiome of human milk itself.
breastfeeding is associated with a decreased incidence of gastrointestinal  tract infections
<cit> , which is corroborated by several studies that have correlated breastfeeding with a lower incidence of necrotizing enterocolitis in humans and animal models
<cit> .
breastfeeding is also associated with an altered fecal microbiome; two studies showed at two weeks of age over 90% of the total fecal bacteria of a breast-fed  infant is bifidobacteria, whereas in most formula-fed  infants bifidobacteria is non-detectable
<cit> .
because the community of gut-colonizing bacteria prevents adhesion and colonization of pathogenic bacteria whilst stimulating mucosal cell proliferation and enhancing immune development, the types of predominant bacteria in the fecal microbiome of the developing infant can affect the health outcomes of the individual, as has been discussed in a recent review article
<cit> .
human milk, the infant's first food, is a primary source of ingested microbiota.
therefore, it is paramount to fully understand the human milk microbiome and how it might influence colonization of the infant gi tract.
ingestion of viable bacteria in human milk may lead to effective colonization of the infant gi tract, but the presence of bacterial dna alone may also hold responsibility for proper infant immune development.
for example, unmethylated cytosine phosphate guanine  dinucleotides within bacterial dna are known as potent immune stimulators, acting through toll-like receptor 9
<cit> .
conversely, immune suppressive motifs including poly-guanosine or guanosine cytosine-rich sequences, such as those on the telomere region of mammalian dna, that can block immune activation induced by cpgs
<cit> .
recently, immune suppressive motifs  that are able to counter the effects of cpgs have been discovered in lactobacillus <cit> .
if immune-modulatory motifs occur in human milk derived dna, they could contribute to proper immune development by decreasing exaggerated inflammatory responses to colonizing bacteria, which are seen in infants with necrotizing enterocolitis
<cit> .
human milk bacteria have previously been analyzed by culture-dependent and -independent mechanisms, confirming the presence of a magnitude of bacterial phylotypes
<cit> .
in one study, staphylococcus and streptococcus dominated the milk microbiome of most mothers, whereas commercially well known bovine milk-associated genera, lactobacillus and bifidobacterium, contributed as minor milk microbiota members
<cit> .
another study showed that the human milk microbiome changes over time, and may be dependent on the mother’s weight and the baby’s mode of delivery
<cit> .
most recent methods for determining the milk microbiome have included amplification of 16s ribosomal rna genes  followed by pyrosequencing
<cit> .
although this technique is widely accepted as a means to determine microbial diversity, it does present limitations such as a lack of information on the functional capacity of the microbes within the milk matrix and also prevents data accumulation on the types of dna motifs to which an infant is exposed.
in this study we performed a metagenomic analysis of the bacteria in human milk using illumina sequencing and the mg-rast pipeline
<cit> .
the aims were to determine the genera of bacteria in human milk, search for immune-modulatory dna motifs, and determine the types of bacterial open reading frames  in human milk that may influence bacterial presence and stability in this complex yet foundational food matrix.
RESULTS
phyla and genera within human milk
metagenome sequencing of a pooled human milk sample resulted in  <dig> , <dig> sequenced reads of 51 bp, which were binned into those aligning to the human genome , known prokaryotic genomes  or those not aligning to either category .
using a best hit analysis of the  <dig> , <dig> 51-bp sequences, 75% aligned to staphylococcus, 15% to pseudomonas, 2% to edwardsiella, and 1% to pantoea, treponema, streptococcus and campylobacter, respectively .
the remaining 3% of the known prokaryotic sequences mapped to  <dig> bacterial genera, demonstrating the diversity of the human milk metagenome while confirming the presence of key genera like akkermansia .
sequences not aligning to prokaryotic or human genomes with a ≤2 bp mismatch were re-aligned to the human genome with decreased stringency , leaving  <dig> , <dig> sequences for contig assembly .
using ray v <dig>
<cit> ,  <dig>  contigs were assembled and submitted to the mg-rast pipeline
<cit> .
post quality control,  <dig>  sequences , with a mean length of 160 ± 55 bp, were used for further analysis .
when the contigs were analyzed using a best hit approach through mg-rast, they aligned predominantly to the phyla of proteobacteria  and firmicutes .
the contigs aligned to  <dig> known genomes at the genus level, predominantly pseudomonas , staphylococcus  and streptococcus , with the highest level of diversity at the genus level within the proteobacteria phylum .
these results are similar to the best hit analysis performed with the non-assembled sequences in that the majority of sequences are from staphylococcus and pseudomonas, but differ in their proportion .
contigs matching viral genomes were observed , including phages derived from pseudomonas and staphylococcus .
contigs also aligned to the genomes of humans, gorillas, chimps and orangutans, likely due to the 60% identity criteria used .
the observation of some of the genera, including staphylococcus, pseudomonas and pantoea, was further validated through the presence of their rrna orfs .
<dig> all sequences not matching the human genome .
open reading frames within human milk
a total of  <dig>  orfs were predicted using mg-rast, of which 82% were annotated , and 18% were unrecognized .
a total of  <dig>  orfs corresponded to a functional category .
for example, many orfs encoded proteins for basic cellular function, including those for respiration , cell signaling , rna , dna , and amino acid metabolism .
orfs encoding proteins for carbohydrate metabolism  included those for lactose metabolism , but none for human milk oligosaccharide metabolism , likely due to the lack of sequences aligning to the genome of bifidobacteria .
virulence-related orfs  included those for antibiotic resistance , adhesion , bacteriocins , as well as others .
stress-related orfs  included those for oxidative stress , osmotic stress , heat and cold shock  and many others .
human milk metagenome compared to mothers’ and infants’ feces
the metagenome of human milk was compared to that of feces from  <dig> unrelated infants  and three unrelated mothers .
using a best hit analysis at the phylum level, contigs from human milk were dissimilar from contigs from feces in regards to the lack of diversity within the human milk metagenome, as over 99% of the contigs were from just two phyla, proteobacteria and firmicutes .
bf-infants’ feces had a high proportion of actinobacteria , followed by ff-infants’ feces , mothers’ feces , and human milk .
the proportion of proteobacteria in the human milk metagenome  was most similar to that of bf-infants’ feces , but was significantly different from ff-infants’ feces and mothers’ feces .
the metagenomes of ff-infants’ feces and mothers’ feces were most similar in regards to their high proportion of bacteroidetes .
conversely, when using a lowest common ancestor approach at the phylum level in comparison to the best hit analysis, human milk appeared more similar to the fecal metagenomes in terms of an increase in diversity , but was still dominated by proteobacteria .
also, using the lowest common ancestor analysis increased the proportion of contigs aligning to actinobacteria in human milk , as well as in mothers’ feces .
the metagenomes of human milk and feces were also compared at the functional level .
the functional orf profile of the human milk metagenome is similar to that of each fecal metagenome, but two fecal profiles were even more similar, for example bf- versus ff-infants’ feces, as seen using pair-wise comparison plots .
the human milk metagenome is most dissimilar from that of ff-infants’ feces as  <dig> out of the  <dig> functional categories contain a significantly different proportion of the orfs .
the three fecal metagenomes had a significantly higher proportion of orfs encoding genes for dormancy and sporulation  than did the human milk metagenome .
both bf- and ff-infants’ fecal metagenomes had significantly higher proportions of cell division  and phosphorus metabolism related orfs  than did the human milk metagenome .
the human milk metagenome, in comparison to bf- and ff-infants’ feces, did, however, have significantly higher proportions of membrane transport , nitrogen  and rna metabolism , cell regulation , respiration , stress response  and virulence-related orfs .
immune-modulatory dna motifs in human milk and feces
when contigs were searched for the presence of immune suppressive motifs, tcaagcttga was found in  <dig> % of the human-milk assembled contigs  with an occurrence  <dig>  times that of the human genome alone .
the contigs positive for tcaagcttga aligned to the genomes of pseudomonas , nocardia , staphylococcus  and contigs of unknown origin .
when the contigs from bf-infants’ feces, ff-infants’ feces and mothers’ feces were scanned for tcaagcttga, it was found at a relative occurrence of  <dig> ,  <dig> , and  <dig>  times that in the human genome, respectively .
another immune suppressive site, ttaggg was observed  <dig>  times in the human milk metagenome , and at a relative occurrence  <dig>  times that of the human genome .
contigs containing ttaggg corresponded to genomes of staphylococcus , pseudomonas , lactobacillus ,  <dig> other known prokaryotic genomes , and contigs from unknown genomes .
when the contigs from bf-infants’ feces, ff-infants’ feces and mothers’ feces were scanned for ttaggg, this sequence was observed at a relative occurrence of  <dig> ,  <dig>  and  <dig>  times that in the human genome, respectively .
assembled contigs were also searched for the presence of synthetically-assembled immune suppressive or immune stimulatory dna motifs , such as those used in vaccine production .
no synthetically-assembled sequences were observed in the human-milk contigs, whereas three motifs were found in less than  <dig> × 10-4% of contigs from the fecal metagenomes .
dna motifs ttaggg and tcaagcttga were searched for in contigs derived from human milk, breast-fed infants’ feces , formula-fed infants’ feces  and mothers’ feces.
relative occurrence is in comparison to the human genome.
discussion
genera within human milk
determining the human milk metagenome, a bodily fluid notably absent from the human microbiome project
<cit> , is crucial for enabling better insight on the process of infant gi colonization and immune development.
by pooling dna from ten human milk samples and subjecting it to illumina sequencing we have demonstrated the large diversity of the human milk metagenome with over  <dig>  contigs aligning to  <dig> bacterial genera .
previous studies investigating the microbiome of human milk have used both culture-dependent and -independent approaches.
using 16s rrna sequencing, hunt et al.
have reported several predominant species in human milk including a core of genera found in  <dig> human milk samples across time: streptococcus, staphylococcus, serratia, pseudomonas, corynebacteria, ralstonia, propionibacteria, sphingomonas, and bradyrhizobiaceae <cit> .
other studies showed colostrum was populated mostly by weisella and leuconostoc, followed by staphylococcus, streptococcus, and lactococcus, and that akkermansia were more prevalent in overweight mothers
<cit> .
using a best hit analysis of the 51 bp illumina reads, alignments for akkermansia, propionibacteria, sphingomonas and weisella were observed , but because of the small number of base pairs used for the alignment  and the lack of assembled contigs associated with these microbes, their presence in our milk samples is a tentative identification.
using pcr-denaturing gradient gel electrophoresis and quantitative pcr, two studies from martin et al.
reported the presence of bifidobacterium breve, b. adolescentis, b. bifidum and b. dentium in human milk, which differs from our findings .
this is likely due to the method of dna extraction used in our study, as we did not incorporate bead-beating as a means to extract dna from the hard to rupture bifidobacterium <cit> .
the differences between the previously reported microbial communities and our analysis may also be due, in part, to the geographic location of the mothers, which has been shown to greatly impact the microbiome of individuals
<cit> .
furthermore, other differences between our characterization of the milk microbiome and those previously reported may be attributed to the means of milk collection.
in comparison to previous studies where human milk was expressed from an aseptic breast
<cit> , the current method determines the total microbiome  ingested by the infant , indicative of what an infant would receive from its mother during suckling.
because our samples were collected from a non-sterilized breast, it could be hypothesized the human milk metagenome reported here would be similar to that of the skin microbiome.
although no reference database was freely available within mg-rast for comparison, the metagenome of human milk is similar to previously reported skin profiles in that there is a large proportion of staphylococcus, which is found in moist areas of skin.
these moist areas, such as the antecubital fossa , also contain betaproteobacteria, such as burkholderia and bordetella, which are present in the milk metagenome .
the human milk metagenome is also similar to drier areas of the skin such as the plantar heel, which contains gamaproteobacteria such as pseudomonas <cit> .
the human milk metagenome is, however, more similar to fecal microbiomes  due to the large proportion of firmicutes bacteria within human milk, which is a very minor member of the skin microbiome .
also, the skin of adults tends to contain a high level of propionibacteria, which notably tends to colonize the skin of cesarean-section birthed babies, whereas this genus is minimally represented in our human milk sample using a best hit analysis of the 51 bp illumina reads .
this observation suggests that mother's milk may prove useful as a skin lotion, to re-balance the skin microbiome of c-section babies.
phylogenetic differences between human milk and feces
comparing the metagenome of human milk to that of publicly available infants’ and mothers’ fecal profiles provides insight as to how human milk may lead to proper colonization of the infant gut.
when comparing the human milk metagenome to the infant fecal metagenome, there are numerous differences.
for example, the metagenome of bf-infants’ feces contains a high proportion of actinobacteria , which correlates with previous studies demonstrating a high abundance of bifidobacterium in the feces of bf-infants whereas ff-infants had a more varied microbiota
<cit> .
contigs from human milk, however, aligned mostly with proteobacteria and firmicutes .
at the phylum level, the present milk metagenome was less diverse than the fecal metagenomes as over 99% of the contigs were from just two phyla, proteobacteria and firmicutes .
ff-infants’ feces and mothers’ feces were similar in that they both contained contigs aligning to the phylum bacteroidetes , whereas bacteroidetes was a very minor component of bf-infants’ feces and human milk .
also, the similar proportion of firmicutes in human milk compared to mothers’ feces  correlates with the hypothesis that mothers’ milk may be inoculated by immune cells carrying bacteria from the gi tract of the mother to her breast
<cit> .
this may be a mechanism by which the human milk microbiome is shaped by the general health of the mother, including her weight
<cit> .
functionality of the human milk metagenome
using illumina sequencing of all dna within milk samples permits the prediction of orfs within assembled contigs and allows for determination of the functional capability of the milk metagenome.
a total of  <dig>  orfs were predicted, including those for basic cell function, as well as those that may enable the bacteria to remain in human milk, such as orfs for carbohydrate metabolism .
the predominant carbohydrate in human milk, lactose, is a potential carbon source for human milk bacteria, and therefore the presence of orfs associated with its metabolism  is expected.
another carbon source for bacteria in human milk is human milk oligosaccharides , which cannot be digested by the infant
<cit> .
these oligosaccharides, which are heavily fucosylated and readily digested by bifidobacteria, are thought to be responsible for the colonization of bf-infants with high levels of bifidobacteria <cit> .
due to a lack of contigs aligning to bifidobacteria , no orfs encoding genes for hmos were observed .
recently, hmos have also been correlated with increased abundance of staphylococcus within human milk, regardless of their inability to utilize the human milk oligosaccharides as a carbon source
<cit> .
the predominance of staphylococcus-aligning contigs in our milk samples supports these findings .
furthermore, there was a significantly higher number of orfs related to nitrogen metabolism within the human milk metagenome in comparison to bf- and ff-infants’ feces .
because human milk contains  <dig> - <dig>  g of nitrogen per 100 g of milk, the bacteria within human milk may use it as a nutrient source in addition to lactose and hmos
<cit> .
human milk contains an abundance of immune cells, antibodies and antimicrobial proteins , and therefore the bacteria residing within human milk must harbor mechanisms to combat the milk-endogenous immune system
<cit> .
for example, the metagenome of human milk includes orfs for stress response and defense  including those for oxidative stress  and toxic compound resistance .
the human milk metagenome also contains orfs for both heat and cold shock , which may enable the bacteria to persist in milk post-breast pumping, refrigeration and reheating.
this may be of particular importance as human milk banks gain more popularity over time.
for example, as described in a recent review by urbaniak et al., some milk banks deem pasteurization of breast milk unnecessary, while others have an upper limit of  <dig> organisms per ml
<cit> .
in unpasteurized banked milk and in-home stored milk, if some organisms are able to survive the storage and re-heating process better than others, the bacterial profile of human milk may change to favor better surviving  bacteria.
furthermore, orfs encoding genes related to virulence and disease , are also observed in the human milk metagenome.
these orfs could allow some of the human milk microbes, such as staphylococcus aureus, to cause mastitis in humans when the balance of human milk-antimicrobials to microbes is tilted towards microbial growth
<cit> .
for example, some bacteria within human milk harbor antibiotic resistance genes  allowing them to proliferate regardless of the mother's potential antibiotic use, and some bacteria are able to produce bacteriocins , which could impact the growth of other, less virulent, microbes within the community.
immune-modulatory landscape of the human milk metagenome
because human milk contains a broad spectrum of microbes at the genus level , it likely contributes significantly towards effective colonization of the infant gi tract.
in the case of banked human milk, which is holder pasteurized , most bacteria are destroyed, but their proteins and dna remain
<cit> .
the presence of non-viable bacteria and bacterial dna in human milk, which are indistinguishable from live bacteria using our approach of dna isolation and sequencing, may be a way to prime the infant immune system and lead to tolerance of the trillions of bacteria that will inhabit the gut following birth.
for example, the immune suppressive motifs, ttaggg and tcaagcttga
<cit> , are present in  <dig> % and  <dig> % of the  <dig>  human milk-contigs, respectively .
the occurrence of the immune suppressive motifs is similar to that in the metagenomes of bf- and ff infants’ feces, as well as mothers’ feces.
this suggests that having a diverse community of microbes may lead to a similar abundance of immune suppressive motifs, regardless of the genera present in the sample.
interestingly, the immune suppressive motif ttaggg was found in higher abundance in the human genome than in bacterial contigs .
colostrum and mature human milk contain between  <dig> ×  <dig> to  <dig> ×  <dig> leukocytes/l and between  <dig> ×  <dig> to  <dig> ×  <dig> leukocytes/l, respectively, which are mostly macrophages  and neutrophils , with natural killer cells representing up to 12% of the population
<cit> .
this suggests that ingestion of the mothers’ dna, through ingestion of her immune cells and any free circulating dna may also lead to proper immune development through a balance of concomitant exposure to immune stimulatory bacterial cpgs and immune suppressive dna in the mothers’ genome and bacterial genomes.
CONCLUSIONS
current microbiome studies characterizing the microbial communities of various anatomical niches have revealed vast differences between healthy individuals
<cit> .
these differences can often be attributed to the host’s environment and diet.
as demonstrated previously by preliminary 16s rrna sequencing, the human milk microbiome is similar to other areas of the body in that its composition is unique to each individual
<cit> .
milk has evolved as the first nutrient source for mammals ex utero, with a high level of inter-mother diversity as to the proportions of bacterial genera, immune proteins and nutrients within it
<cit> .
perhaps, it is the diversity and/or sequences of dna within the milk metagenome that is beneficial for infants, as opposed to any one specific bacterial genus or species.
recent reviews on human milk outline the phylotypes of bacteria within human milk, but only speculate on the function of the human milk microbiome due to a lack of data on the functional capacity of the microbes within human milk
<cit> .
because of this, we sought to better understand the human milk metagenome on a functional level rather than a solely phylogenetic level.
the discovery of the abundance of immune suppressive dna motifs observed within bacterial and human dna from human milk, as well as orfs within the human milk metagenome that allow bacteria to persist in the biological fluid provides a first glance into the functionality of the milk metagenome.
further studies should include those determining the efficacy of milk dna to modulate the immune system in the gi tract, and a more exhaustive look at the metagenome of human milk and how it relates to infant health outcomes.
note added in proof
during revision of the manuscript, everard et al published a report suggesting akkermansia, a human mucus colonizer, helps control diet-induced obesity.
everard et al,  <dig>  proc natl acad sci usa doi/ <dig> /pnas.
<dig>

 the benefits of human milk compared to the use of commercial infant formulas are largely realized because of its bioactive components, including prebiotics, immune proteins and the microbiome of human milk itself. breastfeeding is associated with a decreased incidence of gastrointestinal tract infections, which is corroborated by several studies that have correlated breastfeeding with a lower incidence of necrotizing enterocolitis in humans and animal models. breastfeeding is also associated with an altered fecal microbiome; two studies showed at two weeks of age over 90% of the total fecal bacteria of a breast-fed infant is bifidobacteria, whereas in most formula-fed infants bifidobacteria is non-detectable. because the community of gut-colonizing bacteria prevents adhesion and colonization of pathogenic bacteria whilst stimulating mucosal cell proliferation and enhancing immune development, the types of predominant bacteria in the fecal microbiome of the developing infant can affect the health outcomes of the individual, as has been discussed in a recent review article. human milk, the infant's first food, is a primary source of ingested microbiota. therefore, it is paramount to fully understand the human milk microbiome and how it might influence colonization of the infant gi tract. ingestion of viable bacteria in human milk may lead to effective colonization of the infant gi tract, but the presence of bacterial dna alone may also hold responsibility for proper infant immune development. for example, unmethylated cytosine phosphate guanine dinucleotides within bacterial dna are known as potent immune stimulators, acting through toll-like receptor 9. conversely, immune suppressive motifs including poly-guanosine or guanosine cytosine-rich sequences, such as those on the telomere region of mammalian dna, that can block immune activation induced by cpgs. recently, immune suppressive motifs that are able to counter the effects of cpgs have been discovered in lactobacillus. if immune-modulatory motifs occur in human milk derived dna, they could contribute to proper immune development by decreasing exaggerated inflammatory responses to colonizing bacteria, which are seen in infants with necrotizing enterocolitis. human milk bacteria have previously been analyzed by culture-dependent and -independent mechanisms, confirming the presence of a magnitude of bacterial phylotypes. in one study, staphylococcus and streptococcus dominated the milk microbiome of most mothers, whereas commercially well known bovine milk-associated genera, lactobacillus and bifidobacterium, contributed as minor milk microbiota members. another study showed that the human milk microbiome changes over time, and may be dependent on the mother’s weight and the baby’s mode of delivery. most recent methods for determining the milk microbiome have included amplification of 16s ribosomal rna genes followed by pyrosequencing. although this technique is widely accepted as a means to determine microbial diversity, it does present limitations such as a lack of information on the functional capacity of the microbes within the milk matrix and also prevents data accumulation on the types of dna motifs to which an infant is exposed. in this study we performed a metagenomic analysis of the bacteria in human milk using illumina sequencing and the mg-rast pipeline. the aims were to determine the genera of bacteria in human milk, search for immune-modulatory dna motifs, and determine the types of bacterial open reading frames in human milk that may influence bacterial presence and stability in this complex yet foundational food matrix.  phyla and genera within human milk metagenome sequencing of a pooled human milk sample resulted in, sequenced reads of 51 bp, which were binned into those aligning to the human genome, known prokaryotic genomes or those not aligning to either category. using a best hit analysis of the, 51-bp sequences, 75% aligned to staphylococcus, 15% to pseudomonas, 2% to edwardsiella, and 1% to pantoea, treponema, streptococcus and campylobacter, respectively. the remaining 3% of the known prokaryotic sequences mapped to bacterial genera, demonstrating the diversity of the human milk metagenome while confirming the presence of key genera like akkermansia. sequences not aligning to prokaryotic or human genomes with a ≤2 bp mismatch were re-aligned to the human genome with decreased stringency, leaving, sequences for contig assembly. using ray v, contigs were assembled and submitted to the mg-rast pipeline. post quality control, sequences, with a mean length of 160 ± 55 bp, were used for further analysis. when the contigs were analyzed using a best hit approach through mg-rast, they aligned predominantly to the phyla of proteobacteria and firmicutes. the contigs aligned to known genomes at the genus level, predominantly pseudomonas, staphylococcus and streptococcus, with the highest level of diversity at the genus level within the proteobacteria phylum. these results are similar to the best hit analysis performed with the non-assembled sequences in that the majority of sequences are from staphylococcus and pseudomonas, but differ in their proportion. contigs matching viral genomes were observed, including phages derived from pseudomonas and staphylococcus. contigs also aligned to the genomes of humans, gorillas, chimps and orangutans, likely due to the 60% identity criteria used. the observation of some of the genera, including staphylococcus, pseudomonas and pantoea, was further validated through the presence of their rrna orfs. all sequences not matching the human genome. open reading frames within human milk a total of orfs were predicted using mg-rast, of which 82% were annotated, and 18% were unrecognized. a total of orfs corresponded to a functional category. for example, many orfs encoded proteins for basic cellular function, including those for respiration, cell signaling, rna, dna, and amino acid metabolism. orfs encoding proteins for carbohydrate metabolism included those for lactose metabolism, but none for human milk oligosaccharide metabolism, likely due to the lack of sequences aligning to the genome of bifidobacteria. virulence-related orfs included those for antibiotic resistance, adhesion, bacteriocins, as well as others. stress-related orfs included those for oxidative stress, osmotic stress, heat and cold shock and many others. human milk metagenome compared to mothers’ and infants’ feces the metagenome of human milk was compared to that of feces from unrelated infants and three unrelated mothers. using a best hit analysis at the phylum level, contigs from human milk were dissimilar from contigs from feces in regards to the lack of diversity within the human milk metagenome, as over 99% of the contigs were from just two phyla, proteobacteria and firmicutes. bf-infants’ feces had a high proportion of actinobacteria, followed by ff-infants’ feces, mothers’ feces, and human milk. the proportion of proteobacteria in the human milk metagenome was most similar to that of bf-infants’ feces, but was significantly different from ff-infants’ feces and mothers’ feces. the metagenomes of ff-infants’ feces and mothers’ feces were most similar in regards to their high proportion of bacteroidetes. conversely, when using a lowest common ancestor approach at the phylum level in comparison to the best hit analysis, human milk appeared more similar to the fecal metagenomes in terms of an increase in diversity, but was still dominated by proteobacteria. also, using the lowest common ancestor analysis increased the proportion of contigs aligning to actinobacteria in human milk, as well as in mothers’ feces. the metagenomes of human milk and feces were also compared at the functional level. the functional orf profile of the human milk metagenome is similar to that of each fecal metagenome, but two fecal profiles were even more similar, for example bf- versus ff-infants’ feces, as seen using pair-wise comparison plots. the human milk metagenome is most dissimilar from that of ff-infants’ feces as out of the functional categories contain a significantly different proportion of the orfs. the three fecal metagenomes had a significantly higher proportion of orfs encoding genes for dormancy and sporulation than did the human milk metagenome. both bf- and ff-infants’ fecal metagenomes had significantly higher proportions of cell division and phosphorus metabolism related orfs than did the human milk metagenome. the human milk metagenome, in comparison to bf- and ff-infants’ feces, did, however, have significantly higher proportions of membrane transport, nitrogen and rna metabolism, cell regulation, respiration, stress response and virulence-related orfs. immune-modulatory dna motifs in human milk and feces when contigs were searched for the presence of immune suppressive motifs, tcaagcttga was found in % of the human-milk assembled contigs with an occurrence times that of the human genome alone. the contigs positive for tcaagcttga aligned to the genomes of pseudomonas, nocardia, staphylococcus and contigs of unknown origin. when the contigs from bf-infants’ feces, ff-infants’ feces and mothers’ feces were scanned for tcaagcttga, it was found at a relative occurrence of,, and times that in the human genome, respectively. another immune suppressive site, ttaggg was observed times in the human milk metagenome, and at a relative occurrence times that of the human genome. contigs containing ttaggg corresponded to genomes of staphylococcus, pseudomonas, lactobacillus, other known prokaryotic genomes, and contigs from unknown genomes. when the contigs from bf-infants’ feces, ff-infants’ feces and mothers’ feces were scanned for ttaggg, this sequence was observed at a relative occurrence of, and times that in the human genome, respectively. assembled contigs were also searched for the presence of synthetically-assembled immune suppressive or immune stimulatory dna motifs, such as those used in vaccine production. no synthetically-assembled sequences were observed in the human-milk contigs, whereas three motifs were found in less than × 10-4% of contigs from the fecal metagenomes. dna motifs ttaggg and tcaagcttga were searched for in contigs derived from human milk, breast-fed infants’ feces, formula-fed infants’ feces and mothers’ feces. relative occurrence is in comparison to the human genome. discussion genera within human milk determining the human milk metagenome, a bodily fluid notably absent from the human microbiome project, is crucial for enabling better insight on the process of infant gi colonization and immune development. by pooling dna from ten human milk samples and subjecting it to illumina sequencing we have demonstrated the large diversity of the human milk metagenome with over contigs aligning to bacterial genera. previous studies investigating the microbiome of human milk have used both culture-dependent and -independent approaches. using 16s rrna sequencing, hunt et al. have reported several predominant species in human milk including a core of genera found in human milk samples across time: streptococcus, staphylococcus, serratia, pseudomonas, corynebacteria, ralstonia, propionibacteria, sphingomonas, and bradyrhizobiaceae. other studies showed colostrum was populated mostly by weisella and leuconostoc, followed by staphylococcus, streptococcus, and lactococcus, and that akkermansia were more prevalent in overweight mothers. using a best hit analysis of the 51 bp illumina reads, alignments for akkermansia, propionibacteria, sphingomonas and weisella were observed, but because of the small number of base pairs used for the alignment and the lack of assembled contigs associated with these microbes, their presence in our milk samples is a tentative identification. using pcr-denaturing gradient gel electrophoresis and quantitative pcr, two studies from martin et al. reported the presence of bifidobacterium breve, b. adolescentis, b. bifidum and b. dentium in human milk, which differs from our findings. this is likely due to the method of dna extraction used in our study, as we did not incorporate bead-beating as a means to extract dna from the hard to rupture bifidobacterium. the differences between the previously reported microbial communities and our analysis may also be due, in part, to the geographic location of the mothers, which has been shown to greatly impact the microbiome of individuals. furthermore, other differences between our characterization of the milk microbiome and those previously reported may be attributed to the means of milk collection. in comparison to previous studies where human milk was expressed from an aseptic breast, the current method determines the total microbiome ingested by the infant, indicative of what an infant would receive from its mother during suckling. because our samples were collected from a non-sterilized breast, it could be hypothesized the human milk metagenome reported here would be similar to that of the skin microbiome. although no reference database was freely available within mg-rast for comparison, the metagenome of human milk is similar to previously reported skin profiles in that there is a large proportion of staphylococcus, which is found in moist areas of skin. these moist areas, such as the antecubital fossa, also contain betaproteobacteria, such as burkholderia and bordetella, which are present in the milk metagenome. the human milk metagenome is also similar to drier areas of the skin such as the plantar heel, which contains gamaproteobacteria such as pseudomonas. the human milk metagenome is, however, more similar to fecal microbiomes due to the large proportion of firmicutes bacteria within human milk, which is a very minor member of the skin microbiome. also, the skin of adults tends to contain a high level of propionibacteria, which notably tends to colonize the skin of cesarean-section birthed babies, whereas this genus is minimally represented in our human milk sample using a best hit analysis of the 51 bp illumina reads. this observation suggests that mother's milk may prove useful as a skin lotion, to re-balance the skin microbiome of c-section babies. phylogenetic differences between human milk and feces comparing the metagenome of human milk to that of publicly available infants’ and mothers’ fecal profiles provides insight as to how human milk may lead to proper colonization of the infant gut. when comparing the human milk metagenome to the infant fecal metagenome, there are numerous differences. for example, the metagenome of bf-infants’ feces contains a high proportion of actinobacteria, which correlates with previous studies demonstrating a high abundance of bifidobacterium in the feces of bf-infants whereas ff-infants had a more varied microbiota. contigs from human milk, however, aligned mostly with proteobacteria and firmicutes. at the phylum level, the present milk metagenome was less diverse than the fecal metagenomes as over 99% of the contigs were from just two phyla, proteobacteria and firmicutes. ff-infants’ feces and mothers’ feces were similar in that they both contained contigs aligning to the phylum bacteroidetes, whereas bacteroidetes was a very minor component of bf-infants’ feces and human milk. also, the similar proportion of firmicutes in human milk compared to mothers’ feces correlates with the hypothesis that mothers’ milk may be inoculated by immune cells carrying bacteria from the gi tract of the mother to her breast. this may be a mechanism by which the human milk microbiome is shaped by the general health of the mother, including her weight. functionality of the human milk metagenome using illumina sequencing of all dna within milk samples permits the prediction of orfs within assembled contigs and allows for determination of the functional capability of the milk metagenome. a total of orfs were predicted, including those for basic cell function, as well as those that may enable the bacteria to remain in human milk, such as orfs for carbohydrate metabolism. the predominant carbohydrate in human milk, lactose, is a potential carbon source for human milk bacteria, and therefore the presence of orfs associated with its metabolism is expected. another carbon source for bacteria in human milk is human milk oligosaccharides, which cannot be digested by the infant. these oligosaccharides, which are heavily fucosylated and readily digested by bifidobacteria, are thought to be responsible for the colonization of bf-infants with high levels of bifidobacteria. due to a lack of contigs aligning to bifidobacteria, no orfs encoding genes for hmos were observed. recently, hmos have also been correlated with increased abundance of staphylococcus within human milk, regardless of their inability to utilize the human milk oligosaccharides as a carbon source. the predominance of staphylococcus-aligning contigs in our milk samples supports these findings. furthermore, there was a significantly higher number of orfs related to nitrogen metabolism within the human milk metagenome in comparison to bf- and ff-infants’ feces. because human milk contains -  g of nitrogen per 100 g of milk, the bacteria within human milk may use it as a nutrient source in addition to lactose and hmos. human milk contains an abundance of immune cells, antibodies and antimicrobial proteins, and therefore the bacteria residing within human milk must harbor mechanisms to combat the milk-endogenous immune system. for example, the metagenome of human milk includes orfs for stress response and defense including those for oxidative stress and toxic compound resistance. the human milk metagenome also contains orfs for both heat and cold shock, which may enable the bacteria to persist in milk post-breast pumping, refrigeration and reheating. this may be of particular importance as human milk banks gain more popularity over time. for example, as described in a recent review by urbaniak et al., some milk banks deem pasteurization of breast milk unnecessary, while others have an upper limit of organisms per ml. in unpasteurized banked milk and in-home stored milk, if some organisms are able to survive the storage and re-heating process better than others, the bacterial profile of human milk may change to favor better surviving bacteria. furthermore, orfs encoding genes related to virulence and disease, are also observed in the human milk metagenome. these orfs could allow some of the human milk microbes, such as staphylococcus aureus, to cause mastitis in humans when the balance of human milk-antimicrobials to microbes is tilted towards microbial growth. for example, some bacteria within human milk harbor antibiotic resistance genes allowing them to proliferate regardless of the mother's potential antibiotic use, and some bacteria are able to produce bacteriocins, which could impact the growth of other, less virulent, microbes within the community. immune-modulatory landscape of the human milk metagenome because human milk contains a broad spectrum of microbes at the genus level, it likely contributes significantly towards effective colonization of the infant gi tract. in the case of banked human milk, which is holder pasteurized, most bacteria are destroyed, but their proteins and dna remain. the presence of non-viable bacteria and bacterial dna in human milk, which are indistinguishable from live bacteria using our approach of dna isolation and sequencing, may be a way to prime the infant immune system and lead to tolerance of the trillions of bacteria that will inhabit the gut following birth. for example, the immune suppressive motifs, ttaggg and tcaagcttga, are present in % and % of the human milk-contigs, respectively. the occurrence of the immune suppressive motifs is similar to that in the metagenomes of bf- and ff infants’ feces, as well as mothers’ feces. this suggests that having a diverse community of microbes may lead to a similar abundance of immune suppressive motifs, regardless of the genera present in the sample. interestingly, the immune suppressive motif ttaggg was found in higher abundance in the human genome than in bacterial contigs. colostrum and mature human milk contain between × to × leukocytes/l and between × to × leukocytes/l, respectively, which are mostly macrophages and neutrophils, with natural killer cells representing up to 12% of the population. this suggests that ingestion of the mothers’ dna, through ingestion of her immune cells and any free circulating dna may also lead to proper immune development through a balance of concomitant exposure to immune stimulatory bacterial cpgs and immune suppressive dna in the mothers’ genome and bacterial genomes.  current microbiome studies characterizing the microbial communities of various anatomical niches have revealed vast differences between healthy individuals. these differences can often be attributed to the host’s environment and diet. as demonstrated previously by preliminary 16s rrna sequencing, the human milk microbiome is similar to other areas of the body in that its composition is unique to each individual. milk has evolved as the first nutrient source for mammals ex utero, with a high level of inter-mother diversity as to the proportions of bacterial genera, immune proteins and nutrients within it. perhaps, it is the diversity and/or sequences of dna within the milk metagenome that is beneficial for infants, as opposed to any one specific bacterial genus or species. recent reviews on human milk outline the phylotypes of bacteria within human milk, but only speculate on the function of the human milk microbiome due to a lack of data on the functional capacity of the microbes within human milk. because of this, we sought to better understand the human milk metagenome on a functional level rather than a solely phylogenetic level. the discovery of the abundance of immune suppressive dna motifs observed within bacterial and human dna from human milk, as well as orfs within the human milk metagenome that allow bacteria to persist in the biological fluid provides a first glance into the functionality of the milk metagenome. further studies should include those determining the efficacy of milk dna to modulate the immune system in the gi tract, and a more exhaustive look at the metagenome of human milk and how it relates to infant health outcomes. note added in proof during revision of the manuscript, everard et al published a report suggesting akkermansia, a human mucus colonizer, helps control diet-induced obesity. everard et al, proc natl acad sci usa doi/ /pnas. 
when compared to the metagenome of infants’ and mothers’ feces, the human milk metagenome was less diverse at the phylum level, and contained more open reading frames associated with nitrogen metabolism, membrane transport and stress response. the human milk metagenome was also compared to that of breast-fed and formula-fed infants’ feces and mothers’ feces at the phylum level and at a functional level using open reading frame abundance. human milk contains a diverse population of bacteria that likely influences colonization of the infant gastrointestinal tract. in the present study, a metagenomic approach using illumina sequencing of a pooled milk sample was employed to determine the genera of bacteria and the types of bacterial open reading frames in human milk that may influence bacterial establishment and stability in this primal food matrix. additionally, immune-modulatory bacterial-dna motifs were also searched for within human milk. the bacterial community in human milk contained over prokaryotic genera, with sequences aligning predominantly to the phyla of proteobacteria and firmicutes, and the genera of pseudomonas, staphylococcus and streptococcus. from assembled human milk-derived contigs, open reading frames were annotated and assigned to functional categories. the human milk metagenome also contained a similar occurrence of immune-modulatory dna motifs to that of infants’ and mothers’ fecal metagenomes. our results further expand the complexity of the human milk metagenome and enforce the benefits of human milk ingestion on the microbial colonization of the infant gut and immunity. discovery of immune-modulatory motifs in the metagenome of human milk indicates more exhaustive analyses of the functionality of the human milk metagenome are warranted. human milkmicrobiomemetagenomebacteriailluminadnaopen reading framesimmune-modulatory motifsinfant feces
 human milk contains a diverse population of bacteria that likely influences colonization of the infant gastrointestinal tract. recent studies, however, have been limited to characterization of this microbial community by 16s rrna analysis. in the present study, a metagenomic approach using illumina sequencing of a pooled milk sample was employed to determine the genera of bacteria and the types of bacterial open reading frames in human milk that may influence bacterial establishment and stability in this primal food matrix. the human milk metagenome was also compared to that of breast-fed and formula-fed infants’ feces and mothers’ feces at the phylum level and at a functional level using open reading frame abundance. additionally, immune-modulatory bacterial-dna motifs were also searched for within human milk.  the bacterial community in human milk contained over prokaryotic genera, with sequences aligning predominantly to the phyla of proteobacteria and firmicutes, and the genera of pseudomonas, staphylococcus and streptococcus. from assembled human milk-derived contigs, open reading frames were annotated and assigned to functional categories. when compared to the metagenome of infants’ and mothers’ feces, the human milk metagenome was less diverse at the phylum level, and contained more open reading frames associated with nitrogen metabolism, membrane transport and stress response. the human milk metagenome also contained a similar occurrence of immune-modulatory dna motifs to that of infants’ and mothers’ fecal metagenomes.  our results further expand the complexity of the human milk metagenome and enforce the benefits of human milk ingestion on the microbial colonization of the infant gut and immunity. discovery of immune-modulatory motifs in the metagenome of human milk indicates more exhaustive analyses of the functionality of the human milk metagenome are warranted. human milkmicrobiomemetagenomebacteriailluminadnaopen reading framesimmune-modulatory motifsinfant feces 