Integrated Analysis of Human Milk Microbiota With Oligosaccharides and Fatty Acids in the CHILD Cohort

Human milk contains many bioactive components that are typically studied in isolation, including bacteria. We performed an integrated analysis of human milk oligosaccharides and fatty acids to explore their associations with milk microbiota. We studied a sub-sample of 393 mothers in the CHILD birth cohort. Milk was collected at 3-4 months postpartum. Microbiota was analyzed by 16S rRNA gene V4 sequencing. Oligosaccharides and fatty acids were analyzed by rapid high-throughput high performance and gas liquid chromatography, respectively. Dimension reduction was performed with principal component analysis for oligosaccharides and fatty acids. Center log-ratio transformation was applied to all three components. Associations between components were assessed using Spearman rank correlation, network visualization, multivariable linear regression, redundancy analysis, and structural equation modeling. -values were adjusted for multiple comparisons. Key covariates were considered, including fucosyltransferase-2 (FUT2) secretor status of mother and infant, method of feeding (direct breastfeeding or pumped breast milk), and maternal fish oil supplement use. Overall, correlations were strongest between milk components of the same type. For example, FUT2-dependent HMOs were positively correlated with each other, and was negatively correlated with other core taxa. Some associations were also observed between components of different types. Using redundancy analysis and structural equation modeling, the overall milk fatty acid profile was significantly associated with milk microbiota composition. In addition, some individual fatty acids [22:6n3 (docosahexaenoic acid), 22:5n3, 20:5n3, 17:0, 18:0] and oligosaccharides (fucosyl-lacto-N-hexaose, lacto-N-hexaose, lacto-N-fucopentaose I) were associated with microbiota α diversity, while others (C18:0, 3'-sialyllactose, disialyl-lacto-N-tetraose) were associated with overall microbiota composition. Only a few significant associations between individual HMOs and microbiota were observed; notably, among mothers using breast pumps prevalence was associated with lower abundances of disialyl-lacto-N-hexaose. Additionally, among non-secretor mothers, was positively correlated with sialylated HMOs. Using multiple approaches to integrate and analyse milk microbiota, oligosaccharides, and fatty acids, we observed several associations between different milk components and microbiota, some of which were modified by secretor status and/or