Variation in Human Milk Composition Is Related to Differences in Milk and Infant Fecal Microbial Communities.

Previously published data from our group and others demonstrate that human milk oligosaccharide (HMOs), as well as milk and infant fecal microbial profiles, vary by geography. However, little is known about the geographical variation of other milk-borne factors, such as lactose and protein, as well as the associations among these factors and microbial community structures in milk and infant feces. Here, we characterized and contrasted concentrations of milk-borne lactose, protein, and HMOs, and examined their associations with milk and infant fecal microbiomes in samples collected in 11 geographically diverse sites. Although geographical site was strongly associated with milk and infant fecal microbiomes, both sample types assorted into a smaller number of community state types based on shared microbial profiles. Similar to HMOs, concentrations of lactose and protein also varied by geography. Concentrations of HMOs, lactose, and protein were associated with differences in the microbial community structures of milk and infant feces and in the abundance of specific taxa. Taken together, these data suggest that the composition of human milk, even when produced by relatively healthy women, differs based on geographical boundaries and that concentrations of HMOs, lactose, and protein in milk are related to variation in milk and infant fecal microbial communities.

Key genetic variants associated with variation of milk oligosaccharides from diverse human populations.

Human milk oligosaccharides (HMO), the third most abundant component of human milk, are thought to be important contributors to infant health. Studies have provided evidence that geography, stage of lactation, and Lewis and secretor blood groups are associated with HMO profile. However, little is known about how variation across the genome may influence HMO composition among women in various populations. In this study, we performed genome-wide association analyses of 395 women from 8 countries to identify genetic regions associated with 19 different HMO. Our data support FUT2 as the most significantly associated (P < 4.23-9 to P < 4.5-70) gene with seven HMO and provide evidence of balancing selection for FUT2. Although polymorphisms in FUT3 were also associated with variation in lacto-N-fucopentaose II and difucosyllacto-N-tetrose, we found little evidence of selection on FUT3. To our knowledge, this is the first report of the use of genome-wide association analyses on HMO.