Human Milk Oligosaccharides (HMOs) are abundant carbohydrates fundamental to infant health and development. Although these oligosaccharides were discovered more than half a century ago, their biosynthesis in the mammary gland remains largely uncharacterized. Here, we used a systems biology framework that integrated glycan and RNA expression data to construct an HMO biosynthetic network and predict glycosyltransferases involved. To accomplish this, we constructed models describing the most likely pathways for the synthesis of the oligosaccharides accounting for >95% of the HMO content in human milk. Through our models, we propose candidate genes for elongation, branching, fucosylation, and sialylation of HMOs. We further explored selected enzyme activities through kinetic assay and their co-regulation through transcription factor analysis. These results provide the molecular basis of HMO biosynthesis necessary to guide progress in HMO research and application with the ultimate goal of understanding and improving infant health and development.Significance statementWith the HMO biosynthesis network resolved, we can begin to connect genotypes with milk types and thereby connect clinical infant, child and even adult outcomes to specific HMOs and HMO modifications. Knowledge of these pathways can simplify the work of synthetic reproduction of these HMOs providing a roadmap for improving infant, child, and overall human health with the specific application of a newly limitless source of nutraceuticals for infants and people of all ages.