Background: The newly assembled Bos taurus genome sequence enables the linkage of bovine milk and lactation data with other mammalian genomes. Results: Using publicly available milk proteome data and mammary expressed sequence tags, 197 milk protein genes and over 6,000 mammary genes were identified in the bovine genome. Intersection of these genes with 238 milk production quantitative trait loci curated from the literature decreased the search space for milk trait effectors by more than an order of magnitude. Genome location analysis revealed a tendency for milk protein genes to be clustered with other mammary genes. Using the genomes of a monotreme (platypus), a marsupial (opossum), and five placental mammals (bovine, human, dog, mice, rat), gene loss and duplication, phylogeny, sequence conservation, and evolution were examined. Compared with other genes in the bovine genome, milk and mammary genes are: more likely to be present in all mammals; more likely to be duplicated in therians; more highly conserved across Mammalia; and evolving more slowly along the bovine lineage. The most divergent proteins in milk were associated with nutritional and immunological components of milk, whereas highly conserved proteins were associated with secretory processes. Conclusions: Although both copy number and sequence variation contribute to the diversity of milk protein composition across species, our results suggest that this diversity is primarily due to other mechanisms. Our findings support the essentiality of milk to the survival of mammalian neonates and the establishment of milk secretory mechanisms more than 160 million years ago.

Lactation persistency and milk production are among the most economically important traits in the dairy industry. In this study, we explored the association of over 6.1 million imputed whole-genome sequence variants with lactation persistency (LP), milk yield (MILK), fat yield (FAT), fat percentage (FAT%), protein yield (PROT), and protein percentage (PROT%) in North American Holstein cattle. We identified 49, 3991, 2607, 4459, 805, and 5519 SNPs significantly associated with LP, MILK, FAT, FAT%, PROT, and PROT%, respectively. Various known associations were confirmed while several novel candidate genes were also revealed, including ARHGAP35, NPAS1, TMEM160, ZC3H4, SAE1, ZMIZ1, PPIF, LDB2, ABI3, SERPINB6, and SERPINB9 for LP; NIM1K, ZNF131, GABRG1, GABRA2, DCHS1, and SPIDR for MILK; NR6A1, OLFML2A, EXT2, POLD1, GOT1, and ETV6 for FAT; DPP6, LRRC26, and the KCN gene family for FAT%; CDC14A, RTCA, HSTN, and ODAM for PROT; and HERC3, HERC5, LALBA, CCL28, and NEURL1 for PROT%. Most of these genes are involved in relevant gene ontology (GO) terms such as fatty acid homeostasis, transporter regulator activity, response to progesterone and estradiol, response to steroid hormones, and lactation. The significant genomic regions found contribute to a better understanding of the molecular mechanisms related to LP and milk production in North American Holstein cattle.

High-fat diet (HFD) during lactation alters milk composition and is associated with development of metabolic diseases in the offspring. We hypothesized that HFD affects milk microRNA (miRNA) and mRNA content, which potentially impact offspring development. Our objective was to determine the effect of maternal HFD on secreted milk transcriptome. To meet this objective, 4 wk old female ICR mice were divided into two treatments: control diet containing 10% kcal fat and HFD containing 60% kcal fat. After 4 wk on CD or HFD, mice were bred while continuously fed the same diets. On postnatal day 2 (P2), litters were normalized to 10 pups, and half the pups in each litter were cross-fostered between treatments. Milk was collected from dams on P10 and P12. Total RNA was isolated from milk fat fraction of P10 samples and used for mRNA-Seq and small RNA-Seq. P12 milk was used to determine macronutrient composition. After 4 wk of prepregnancy feeding HFD mice weighed significantly more than did the control mice. Lactose and fat concentration were significantly ( P < 0.05) higher in milk of HFD dams. Pup weight was significantly greater ( P < 0.05) in groups suckled by HFD vs. control dams. There were 25 miRNA and over 1,500 mRNA differentially expressed (DE) in milk of HFD vs. control dams. DE mRNA and target genes of DE miRNA enriched categories that were primarily related to multicellular organismal development. Maternal HFD impacts mRNA and miRNA content of milk, if bioactive nucleic acids are absorbed by neonate differences may affect development.