Skip to main content

Lipidome analysis of milk composition in humans, monkeys, bovids, and pigs



Lipids contained in milk are an essential source of energy and structural materials for a growing neonate. Furthermore, lipids’ long-chain unsaturated fatty acid residues can directly participate in neonatal tissue formation. Here, we used untargeted mass spectrometric measurements to assess milk lipid composition in seven mammalian species: humans, two macaque species, cows, goats, yaks, and pigs.


Analysis of the main milk lipid class, triacylglycerides (TAGs), revealed species-specific quantitative differences in the composition of fatty acid residues for each of seven species. Overall, differences in milk lipid composition reflect evolutionary distances among species, with each species group demonstrating specific lipidome features. Among them, human milk contained more medium and long-chain unsaturated fatty acids compared to other species, while pig milk was the most distinct, featuring the highest proportion of long-chain polyunsaturated fatty acids.


We show that milk lipidome composition is dynamic across mammalian species, changed extensively in pigs, and contains features particular to humans.


Milk is the sole source of nutrients to support growth and development in neonates [1, 2]. The nutritional requirements differ among species due to ecological and biological variability of the newborns’ development, mirrored by the variation in milk composition. For example, the composition of marsupial milk adjusts according to the developmental stage of neonates born altricial [3, 4]. In placental mammals, milking strategy and milk composition depend on the species lifestyle and the ability to accompany offspring after birth. If the lactation period is short, then milk will be dense, allowing to transfer all the necessary nutrition to the baby within this short period. Those species that have the opportunity to accompany their offspring for a longer time after birth, such as bovids and primates, produce more diluted milk [5]. In addition to temporal milk composition changes along with the baby’s growth, there is evidence that milk’s microelements and caloricity vary due to baby’s sex [6,7,8].

Besides nutritional components, milk contains steroid and peptide hormones, including leptin, ghrelin, and adiponectin [9], immune and growth factors [2, 10, 11], as well as poly-unsaturated fatty acids directly participating in the neonatal tissue development and growth [12, 13]. The most variable element of milk composition is total fat content ranging from 60% in the hooded seals [14] to less than 1% in white rhinoceroses and ring-tailed lemur [15] – the phenomenon mainly linked to the duration of breastfeeding and ecological conditions. In most species, however, a single lipid class, TAGs, composes on average, 98% of the milk fat [16]. TAGs have simple chemical composition, consisting of glycerol and three fatty acid residues. Several studies analyzed milk TAG composition in mammalian species, including humans [17,18,19,20,21], demonstrating differences in the distribution of docosahexaenoic acid (DHA)-containing TAGs in comparison with the baby formula [22]. However, systematic characterization of the milk lipidome across mammalian species is yet to be performed.

Maternal diet can play an essential role in the milk fatty acid composition - thus, dietary enrichment with particular fatty acids, including the essential ones, will lead to an increase of these fatty acids in the milk TAGs [23,24,25]. Similarly, seasonal changes in diet altering gut microbial community were shown to affect milk fatty acid composition in bovids [26]. Maternal parity was also associated with differential fatty acid content in bovids and humans, favoring primiparous mothers [26, 27]. Another factor that can influence milk composition, as well as the milk yield, is genetic background and variation in the genes associated with fatty acid biosynthesis [28,29,30].

While most lipids contained in milk get broken down to energy sources and simple building blocks, there is evidence that some of TAG fatty acid residues could directly participate in neonatal tissue development and growth. For example, DHA derived from milk TAGs is transported across the blood-brain barrier in the form of lysophosphatidylcholine (LPC) with the help of a specific protein transporter [31]. The direct use of milk fatty acid residues in neonatal tissue development indicates that the milk lipidome composition of a species might reflect specific requirements of its neonatal tissues.

If milk lipids indeed represent the best match for neonate’s nutritional demand, we expect to find corresponding species-specific differences in milk lipidome composition. More specifically, since human neonates are distinct even from non-human primates, we anticipate seeing human-specific milk lipidome features. To address these hypotheses, we used liquid chromatography-mass spectrometry to assess milk TAG composition across seven mammalian species, including humans.


We examined the lipid composition of milk samples in two orders of mammals: artiodactyla, or cloven-hoofed mammals, and primate. Within artiodactyla, our study contains representatives of two families: suidae (pig) and bovidae (cow, domestic yak, and goat). Primate species included cercopithecidae (crab-eating and rhesus monkeys) and hominidae (human) (Fig. 1a). We sampled 23 primate individuals (19 humans and four macaques), 10 bovids (four cows, four goats, and two yaks), and four pigs (Additional file 5: Table S1).

Fig. 1

Milk lipidome evolution. a Dendrogram showing the phylogenetic relationship among seven mammalian species used in the study. Colored dots indicate species, with colors used consistently throughout figures. Colored frames indicate three groups: primate (red), artiodactyla (blue), and bovidae (green). b The relationship among milk samples based on the signal intensities of 472 detected lipid features plotted in two dimensions using a multidimensional scaling algorithm. Colors represent species. Small dots represent individual milk samples. Large dots represent pools of samples from one species. c The relationship between lipid-intensity-based distances between species pairs and phylogenetic distances

Randomized milk samples were extracted, separated using liquid chromatography, and measured using untargeted mass spectrometry in positive ionization mode. The measurements yielded a total of 472 mass spectrometry features representing distinct hydrophobic compounds (lipids) with molecular weights below 1200 Da (Da).

Visualization of the relationship among samples based on the abundance levels of these 472 detected lipids using multidimensional scaling (MDS) revealed good separation of species and phylogenetic groups (Fig. 1b). Furthermore, distances between species calculated using the normalized intensities of mass spectrometry signals generally agreed with the phylogenetic distances (Fig. 1c). Pig milk was the only obvious exception from this linear relationship, showing a greater difference to the bovidae species than expected from the phylogeny.

Of the 472 detected lipids, 403 (85%) showed significant intensity differences among species (Analysis of Variance (ANOVA), Benjamini-Hochberg (BH)-corrected P < 0.05). Unsupervised clustering of these 403 species-dependent lipids based on their intensity profiles across samples yielded four clusters (Fig. 2a). Lipid intensities within these clusters differed within the mammalian orders as much as between them. Notably, the lipid composition of the pig milk stood out in three of the four clusters, while cluster 4 contained milk lipidome composition features shared between monkeys and bovids (Fig. 2b).

Fig. 2

Patterns of lipid composition differences among species. a Unsupervised hierarchical clustering dendrogram and heatmap based on pairwise correlation-based distances between lipid intensity vectors. Color bar indicates four major lipid clusters: cluster 1 – yellow, cluster 2 – pink, cluster 3 – brown, cluster 4 – gray. b Distribution of centered lipid intensities in each cluster. Each box represents the interquartile distribution of lipid intensities in a sample. The number of lipids contained in a cluster is shown at the top of the panels. Samples colored and ordered by species. c Mass spectrometry output showing 76 computationally annotated TAGs colored by clusters. The x-axis shows the compounds’ mass-by-charge ratio (m/z). The y-axis shows retention time (RT) of the compound on the liquid chromatography preceding mass spectrometry

In agreement with previous knowledge, most of the lipids detected in milk belonged to a specific lipid class – TAGs. Interestingly, 76 lipid features computationally annotated as TAGs were present in all four clusters, covering the entire spectrum of lipid variation patterns (Additional file 6: Table S2). Furthermore, the average length and unsaturation extent of the fatty acid chains of these TAGs differed among the clusters (Fig. 2c). For instance, cluster 2 TAGs contained long-chain polyunsaturated fatty acid residues, while cluster 4 TAGs preferentially contained medium-chain fatty acid residues (Fig. 2c; Additional file 1: Figure S1).

Notably, relative abundance analysis of detected TAGs revealed apparent intensity differences characteristic of each species. Specifically, cow milk contained more TAGs composed of long monounsaturated fatty acids. By contrast, goat milk contained more TAGs composed of medium-chain saturated and monounsaturated fatty acids. Pig milk stood out from the rest of the species by having TAGs composed of long- and very-long-chain polyunsaturated fatty acids. Monkey milk had more TAGs composed of medium-chain monounsaturated and polyunsaturated fatty acids. Finally, human milk tended to contain more TAGs with long-chain polyunsaturated fatty acids (Fig. 3a; Additional file 1: Figure S1). The magnitude of normalized intensity differences for the annotated TAGs across species ranged from 1.1 to 1.4 fold (Additional file 7: Table S3).

Fig. 3

Milk TAG repertoire in each species. a Distribution of milk TAGs in each species and two human populations. The x-axis represents the total chain length of three TAG fatty acid residues. The y-axis represents the number of double bonds in TAG fatty acid residues. The size of the circles represents TAG intensity calculated as the mean intensity of the TAG in a species divided by the sum of the mean intensities of all TAGs in this species. Colors represent the number of double bonds with a gradient from orange – zero double bonds to purple – seven double bonds. b Cumulative distribution of species-specific TAGs features. The axes as in panel A. Size of the circles represents the maximum intensity of a TAG across species. Colors indicate species with the maximum intensity of a TAG

We next specifically searched for lipids showing statistically significant intensity differences between humans and the other three species’ groups represented in our study. Of the 472 detected lipids, 94 differed in intensity between humans and macaques, 23 of them annotated as TAGs (ANOVA, BH-corrected P < 0.05). Clustering of these lipids revealed a notable intensity increase in the human milk for a group of lipids, including nine TAGs with long- and very-long-chain fatty acids (Fig. 4a; Additional file 2: Figure S2).

Fig. 4

Comparison of milk TAG profiles between humans and other species. a Left: Distribution of centered intensities of 94 lipids showing differences between humans and two macaque species separated into two clusters. Each box represents the interquartile distribution of lipid intensities in a sample. The number of lipids contained in a cluster is shown at the top of the panels. Samples colored and ordered by species. Right: Mass spectrometry output with 76 computationally annotated TAGs showing the difference between humans and two macaque species colored by clusters. b As panel A, for 76 lipids showing intensity differences between humans and three bovidae species. c As panel A, for 76 lipids showing intensity differences between humans and pigs

Comparison between humans and bovids yielded significant intensity differences for 269 lipids, 61 of them annotated as TAGs (ANOVA, BH-corrected P < 0.05). Among them, 23 TAGs with long- and very-long-chain polyunsaturated fatty acids showed increased intensities in human milk (Fig. 4b; Additional file 3: Figure S3).

Comparison between humans and pigs revealed significant intensity differences for 343 lipids, 64 of them annotated as TAGs (ANOVA, BH-corrected P < 0.05). Among them, 15 TAGs containing medium-chain saturated, monounsaturated, and polyunsaturated fatty acids were increased in the human milk, while the rest of TAGs containing long- and very-long-chain polyunsaturated fatty acids were elevated in the pig milk (Fig. 4c; Additional file 4: Figure S4).


Our study provides the first systematic description of the milk lipid composition of several mammalian species, including primates and, in particular, humans. We show that TAGs, the main components of the milk fat, differ widely among mammalian species in terms of their average fatty acid length and unsaturation degree.

A comparison of the human milk with bovid milk shows that TAGs containing saturated fatty acids are particularly abundant in bovids. In contrast, human milk has a significantly higher content of long-chain polyunsaturated fatty acids with the cumulative number of double bonds ranging from 3 to 8 (Fig. 4b). Pig milk lipid composition is more similar to the human milk lipids when compared to bovids (Fig. 2b, c clusters 2 and 3). Yet, the comparison of human milk lipids to pig milk reveals a higher abundance of long-chain polyunsaturated fatty acids in pigs (Fig. 4c). Of the seven species examined in our study, pig milk shows the most rapid evolutionary divergence, resulting in an unusually high proportion of polyunsaturated fatty acids.

Importantly, we detect specific features of the milk lipidome composition in each of the seven species, including humans. While we were unable to obtain milk samples from apes, the comparison to two old world monkey species - rhesus and crab-eating macaques - revealed milk lipidome features potentially unique to humans. Specifically, we show that after the human-monkey species’ divergence approximately 30 million years ago, human ancestors started to produce milk with a higher abundance of TAGs containing long-chain fatty acids with high levels of unsaturation. This observation is intriguing, given the reports of particular long-chain polyunsaturated fatty acids accumulating in the human brain during the last trimester of pregnancy and after birth [19] potentially influencing brain development and functionality [8]. In addition to TAGs, we detect human-specific intensity differences for lipids representing other lipid classes present in milk, such as phospholipids. However, given the small sample size of the study, we did not consider low-abundance lipid classes in our analysis.

Both pig and human milk display elevated amounts of long-chain polyunsaturated fatty acids. Pigs have a substantially shorter lactation period compared to humans and many other mammals. Thus, their milk composition might reflect the need to provide the most critical nutrition within a short period [5, 6]. Our study supports this hypothesis by showing that pig milk indeed stands out from the rest of the assessed species in terms of the long-chain polyunsaturated fatty acids content. In humans and non-human primates, the lactation strategy is different from the artiodactyls and pigs in particular. Since primates can provide care to their youth for an extended time after birth, reaching to several years, they tend to produce more dilute milk rich in carbohydrates and water, and low in fats.

Nonetheless, we detect an increase in long-chain unsaturated fatty acid content in human milk. Given that the human lactation period is not shorter than that of monkey species, this human milk feature should not reflect the shortening of the lactation period. Instead, it might represent an adaptation for the increased demand in the particular fatty acids raised by growing organs, such as the brain.


Our study revealed substantial differences in TAG composition among seven mammalian species: three primates, three bovids, and pigs. While for most species, changes in milk lipidome composition fit the general evolutionary pattern, with distances proportional to the phylogenetic times, there is an exception. Specifically, pig milk stood out by containing unusually high amounts of long-chain polyunsaturated fatty acids. Notably, human milk was second in terms of long-chain polyunsaturated fatty acids abundance, followed by two macaque species, and then by the bovids. While both pig and human milk contain more long-chain polyunsaturated fatty acids, this increase could represent an adaptation to a shorter lactation period in pigs, but reflect the enhanced demand of the growing brain in humans. These results indicate the need for further studies of milk lipidome evolution in conjunction with other developing tissues, especially the brain.


Sample collection

We collected milk samples from healthy human volunteers representing Russian (n = 9) and Chinese (n = 10), cows (n = 4), goats (n = 4), pigs (n = 4), yaks (n = 2), rhesus monkeys (n = 2), and crab-eating monkeys (n = 2) (Additional file 5: Table S1). Informed consent for the use of milk in this study was obtained from each of the human volunteers. In all species, milk samples were collected at matched lactation stage, several weeks after neonate birth. In each case, milk was sampled at the end of breastfeeding or milking event into the same type of 10-ml plastic container and immediately placed into − 20 °C freezer for short-term storage not exceeding two weeks. Samples were then transported on dry ice without de-freezing into − 80 °C freezer for long-term storage, not exceeding four months.

Ethic statement

All animals were treated in accordance with good animal practice as defined by the local welfare authorities; all human volunteers have signed an Informed Consent Form, confirming that they understand the purpose of the research and their participation is entirely gratuitous.

Lipid extraction

The milk aliquots were thawed at 0 °C mixed, and 16 μl of milk were transferred to a 2.0 ml Eppendorf safe-lock tube and resuspended in 34 μl of LC-MS grade water. Prior to extraction, samples were randomized with regard to species’ identity. Furthermore, for each species, we prepared a pooled sample containing equal volumes of milk from each individual. For lipid extraction, a modified two-phase protocol was used as described in [32]. All manipulations with samples were performed on ice. Briefly, 750 μl of MeOH:MTBE (1:3) solution containing internal standards in concentration of 1 mg/L were added to each sample, vortexed for 1 min, sonicated for 15 min in an ice-cooled sonication bath, incubated for 30 min at 4 °C, and sonicated for the second time in a pre-cooled sonication bath. Then, 560 μl of MeOH:H2O (1:3) solution was added to each sample, vortexed for 10 s and centrifuged for 10 min at 14.000 x g at 4 °C. The 400 μl of the upper-phase, containing organic fraction, was transferred to a new 2.0 ml Eppendorf tube and dried in a Speedvac for 1 h at 30 °C.


Dried lipid pellets were resuspended in 400 μl of acetonitrile:isopropanol (1:3) solution. Samples were rigorously vortexed for 10 s, shaken for 10 min at 4 °C and sonicated in an ice bath. Then, 5 μl of each sample was transferred to the autosampler glass vial and diluted 1:20 with 95 μl of acetonitrile:isopropanol (1:3) solution. A pool of all samples was prepared by mixing 5 μl from each sample in a separate Eppendorf tube, transferred to glass vials and diluted 1:20 with acetonitrile:isopropanol (1:3) solution to get quality control (QC) samples. Sample pools for each species were made by mixing 10 μl of each sample of the corresponding species and diluted 1:20 with acetonitrile:isopropanol (1:3) solution. From each diluted sample, 3 μl were injected to a reversed-phase Bridged Ethylene Hybrid (BEH) C8 reverse column (100 mm × 2.1 mm, containing 1.7 μm diameter particles, Waters) coupled to a Vanguard pre-column with the same dimensions, using a Waters Acquity UPLC system (Waters, Manchester, UK). The mobile phases used for the chromatographic separation were: water, containing 10 mM ammonium acetate, 0.1% formic acid (Buffer A) and acetonitrile:isopropanol (7:3 (v:v)), containing 10 mM ammonium acetate, 0.1% formic acid (Buffer B). The gradient separation was: 1 min 55% B, 3 min linear gradient from 55 to 80% B, 8 min linear gradient from 80% B to 85% B, and 3 min linear gradient from 85% A to 100% A. After 4.5 min washing with 100% B the column was re-equilibrated with 55% B. The flow rate was set to 400 μl/min. The mass spectra were acquired in a positive mode using a heated electrospray ionization source in combination with a Bruker Impact II QTOF (quadrupole-Time-of-Flight) mass spectrometer (Bruker Daltonics, Bremen, Germany).

Four blank samples were run at the beginning of the queue, followed by four QC samples to equilibrate the column. After them, 38 samples were queued in the same random order used for extraction with all samples randomized by species, interleaved with seven pooled samples, one QC preceding the first sample and then a QCs after every 9th sample. At the end of the queue, we performed two injections containing 100% acetonitrile to wash the column, followed by blank samples. Blank samples were prepared as usual samples, but contained only extraction buffers to reveal all contaminants that could come from the extraction and other technical steps, and not from the sample itself.

Data preprocessing

After the acquisition, Bruker raw data .d files were automatically calibrated using the internal calibration and converted into mzXML format using a custom DataAnalysis script (Bruker, Version 4.3). The mzXML files were then subjected to the standard alignment and peak picking procedure using xcms software [33]. We then filtered from the output table all lipid features with the coefficient of variation (CoV, calculated as the standard deviation over the mean across QC samples) > 30%, and peaks with zero values in > 50% of the individual samples. Lipid features’ intensities were then normalized using the upper quartile normalization and base-two log-transformed. Raw data is uploaded to the metabolomics study data repository MetaboLights [34].

Phylogenetic distances calculation

The lipid intensity-based distances between species were calculated as the Euclidean distance between the vectors containing the intensities of 472 lipids detected in each pair of species. The phylogenetic distances between the species’ pairs were obtained from the TimeTree database [35].

Statistical analysis

Species-dependent lipids were defined with ANOVA and BH-corrected p-value cutoff of 0.05. The correlation matrix of species-dependent lipids was calculated as (1-cor) Pearson’s distances between all lipids. Unsupervised clustering of the species-dependent lipids was performed using hierarchical clustering with complete linkage in the R statistical environment.

TAGs annotation

All detected lipid features were annotated against the theoretical list with all possible masses of NH4+ adducts of TAGs. The theoretical list of masses was generated using ALEX Target List Generator with the ALEX lipid database (5.2) [36]. Among the detected lipids, 76 matched the theoretical masses with < 10 ppm and were considered for further analysis.

Species-dependent TAG intensity differences

We calculated the mean intensities of TAGs in each species using the raw intensities of the annotated TAGs. To assess TAGs concentrations in each species, the mean intensity of a particular TAG was divided by the sum of the mean intensities of all TAGs in that species. For comparison between species, the 100% intensity value of the TAG was defined as the maximal intensity of the TAG across the species.

Availability of data and materials

The datasets generated and analysed during the current study are available in the MetaboLights repository,



Analysis of Variance


Bridged Ethylene Hybrid




Coefficient of Variation




Docosahexaenoic acid


Fold change


Liquid chromatography–mass spectrometry




Multidimensional scaling




Methyl tert-butyl ether


Million years ago


Quality control




Retention time




Ultra-performance liquid chromatography


  1. 1.

    Wu X, Jackson RT, Khan SA, Ahuja J, Pehrsson PR. Human Milk Nutrient Composition in the United States: Current Knowledge, Challenges, and Research Needs. Curr Dev Nutr. 2018;2:nzy025.

    PubMed  PubMed Central  Google Scholar 

  2. 2.

    Ballard O, Morrow AL. Human Milk composition. Pediatr Clin N Am. 2013;60:49–74.

    Article  Google Scholar 

  3. 3.

    Nicholas KR. Asynchronous dual lactation in a marsupial, the tammar wallaby (Macropus eugenii). Biochem Biophys Res Commun. 1988;154:529–36.

    CAS  Article  Google Scholar 

  4. 4.

    Trott JF, et al. Maternal regulation of milk composition, milk production, and pouch young development during lactation in the tammar wallaby (Macropus eugenii). Biol Reprod. 2003;68:929–36.

    CAS  Article  Google Scholar 

  5. 5.

    Hinde K, Milligan LA. Primate milk: proximate mechanisms and ultimate perspectives. Evol Anthropol. 2011;20:9–23.

    Article  Google Scholar 

  6. 6.

    Powe CE, Knott CD, Conklin-Brittain N. Infant sex predicts breast milk energy content. Am J Hum Biol. 2010;22:50–4.

    Article  Google Scholar 

  7. 7.

    Galante L, et al. Sex-specific human Milk composition: the role of infant sex in determining early life nutrition. Nutrients. 2018;10.

  8. 8.

    Hinde K. Richer milk for sons but more milk for daughters: sex-biased investment during lactation varies with maternal life history in rhesus macaques. Am J Hum Biol. 2009.

  9. 9.

    Savino F, Petrucci E, Lupica MM, Nanni GE, Oggero R. Assay of ghrelin concentration in infant formulas and breast milk. World J Gastroenterol. 2011;17:1971–5.

    CAS  Article  Google Scholar 

  10. 10.

    Bhinder G, et al. Milk fat globule membrane supplementation in formula modulates the neonatal gut microbiome and normalizes intestinal development. Sci Rep. 2017;7:45274.

    CAS  Article  Google Scholar 

  11. 11.

    Hampel D, Shahab-Ferdows S, Islam MM, Peerson JM, Allen LH. Vitamin concentrations in human Milk vary with time within feed, circadian rhythm, and single-dose supplementation. J Nutr. 2017;147:603–11.

    CAS  Article  Google Scholar 

  12. 12.

    Innis SM. Dietary (n-3) fatty acids and brain development. J Nutr. 2007;137:855–9.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Mulder KA, Elango R, Innis SM. Fetal DHA inadequacy and the impact on child neurodevelopment: a follow-up of a randomised trial of maternal DHA supplementation in pregnancy. Br J Nutr. 2018;119:271–9.

    CAS  Article  Google Scholar 

  14. 14.

    Oftedal OT, Bowen WD, Boness DJ. Energy transfer by lactating hooded seals and nutrient deposition in their pups during the four days from birth to weaning. Physiol Zool. 1993;66:412–36.

    Article  Google Scholar 

  15. 15.

    Power LM, Shulkin J. Milk: the biology of lactation. JHU Press. 2016:120.

  16. 16.

    Walstra P. Dairy technology: principles of Milk properties and processes: CRC Press; 1999. p. 50.

  17. 17.

    Morera S, Castellote AI, Jauregui O, Casals I, López-Sabater MC. Triacylglycerol markers of mature human milk. Eur J Clin Nutr. 2003;57:1621–6.

    CAS  Article  Google Scholar 

  18. 18.

    Haddad I, Mozzon M, Strabbioli R, Frega NG. A comparative study of the composition of triacylglycerol molecular species in equine and human milks. Dairy Sci Technol. 2012;92:37–56.

    CAS  Article  Google Scholar 

  19. 19.

    Beccaria M, et al. High performance characterization of triacylglycerols in milk and milk-related samples by liquid chromatography and mass spectrometry. J Chromatogr A. 2014;1360:172–87.

    CAS  Article  Google Scholar 

  20. 20.

    Linderborg KM, et al. Tandem mass spectrometric analysis of human milk triacylglycerols from normal weight and overweight mothers on different diets. Food Chem. 2014;146:583–90.

    CAS  Article  Google Scholar 

  21. 21.

    Ten-Doménech I, Beltrán-Iturat E, Herrero-Martínez JM, Sancho-Llopis JV, Simó-Alfonso EF. Triacylglycerol analysis in human Milk and other mammalian species: small-scale sample preparation, characterization, and statistical classification using HPLC-ELSD profiles. J Agric Food Chem. 2015;63:5761–70.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Liu Z, Cocks BG, Rochfort S. Comparison of molecular species distribution of DHA-containing Triacylglycerols in Milk and different infant formulas by liquid chromatography–mass spectrometry. J Agric Food Chem. 2016;64:2134–44.

    CAS  Article  Google Scholar 

  23. 23.

    Cimatti AG, et al. Maternal supplementation with krill oil during breastfeeding and long-chain polyunsaturated fatty acids (LCPUFAs) composition of human Milk: a feasibility study. Front Pediatr. 2018;6:407.

    Article  Google Scholar 

  24. 24.

    Kuipers RS, et al. High contents of both docosahexaenoic and arachidonic acids in milk of women consuming fish from lake Kitangiri (Tanzania): targets for infant formulae close to our ancient diet? Prostaglandins Leukot Essent Fatty Acids. 2005;72:279–88.

    CAS  Article  Google Scholar 

  25. 25.

    Martin MA, et al. Fatty acid composition in the mature milk of Bolivian forager-horticulturalists: controlled comparisons with a US sample. Matern Child Nutr. 2012;8:404–18.

    Article  Google Scholar 

  26. 26.

    Hanuš O, Samková E, Křížová L, Hasoňová L, Kala R. Role of fatty acids in Milk fat and the influence of selected factors on their variability-a review. Molecules. 2018;23.

  27. 27.

    Prentice A, Jarjou LM, Drury PJ, Dewit O, Crawford MA. Breast-milk fatty acids of rural Gambian mothers: effects of diet and maternal parity. J Pediatr Gastroenterol Nutr. 1989;8:486–90.

    CAS  Article  Google Scholar 

  28. 28.

    Miliku K, et al. Human milk fatty acid composition is associated with dietary, genetic, sociodemographic, and environmental factors in the CHILD cohort study. Am J Clin Nutr. 2019.

  29. 29.

    Hlusko LJ, et al. Environmental selection during the last ice age on the mother-to-infant transmission of vitamin D and fatty acids through breast milk. Proc Natl Acad Sci U S A. 2018;115:E4426–32.

    CAS  Article  Google Scholar 

  30. 30.

    Mida K, Shamay A, Argov-Argaman N. Elongation and desaturation pathways in mammary gland epithelial cells are associated with modulation of fat and membrane composition. J Agric Food Chem. 2012;60:10657–65.

    CAS  Article  Google Scholar 

  31. 31.

    Nguyen LN, et al. Mfsd2a is a transporter for the essential omega-3 fatty acid docosahexaenoic acid. Nature. 2014;509:503–6.

    CAS  Article  Google Scholar 

  32. 32.

    Sarafian MH, et al. Objective set of criteria for optimization of sample preparation procedures for ultra-high throughput untargeted blood plasma lipid profiling by ultra performance liquid chromatography–mass spectrometry. Anal Chem. 2014;86:5766–74.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Smith CA, Want EJ, O’Maille G, Abagyan R, Siuzdak G. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem. 2006;78:779–87.

    CAS  Article  Google Scholar 

  34. 34.

    Haug K, et al. MetaboLights - an open-access general-purpose repository for metabolomics studies and associated meta-data. Nucleic Acids Res. 2013;41:D781–6.

    CAS  Article  Google Scholar 

  35. 35.

    Hedges B, Kumar S. TimeTree knowledge-base for information on the tree-of-life and its evolutionary timescale. Accessed 9 Apr 2019.

  36. 36.

    Husen P, et al. Analysis of lipid experiments (ALEX): a software framework for analysis of high-resolution shotgun Lipidomics data. PLoS One. 2013;8:e79736.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank Peng Shi for providing yak milk samples; Vadim Bakushkin and Bulatovo Farm for their generous help in goat milk sample collection; Ksenia Fominykh and Svetlana Gribova for organising human milk samples collection; Markus R. Wenk and his team at the National University of Singapore for valuable comments and suggestions.


This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (grant XDB13010200); the National Natural Science Foundation of China (grant 31420103920); the National One Thousand Foreign Experts Plan (grant WQ20123100078); the National Key R&D Program of China (grant 2017YFA0505700); and the Russian Science Foundation (grant 16–14-00220). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




PK and WM conceived and designed the study. SG and AM organized samples collection. AV, NA, WM and AM performed samples preparation, lipid extraction and mass-spectrometry. PM and AM performed statistical analysis and contributed figures. PK and AM wrote the manuscript. All authors read and approved the final version of manuscript.

Corresponding authors

Correspondence to Pavel Mazin or Philipp Khaitovich.

Ethics declarations

Ethics approval and consent to participate

The study was carried out with the ethical approval of the Institutional Review Board of the Skolkovo Institute of Science and Technology. All participants have signed an Informed Consent form.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1. Mass spectrometry output showing 76 lipid features computationally annotated as TAGs. The x-axis shows the compounds’ mass-by-charge ratio (m/z). The y-axis shows retention time (RT) of the compound on the liquid chromatography preceding mass spectrometry. Colors indicate four clusters of lipids showing intensity differences among species. Each point represents annotated TAG feature; labels include the cumulative length of the carbon chains and the total number of double bonds of the fatty acid residues.

Additional file 2: Figure S2. Mass spectrometry output showing 76 lipid features computationally annotated as TAGs. The x-axis shows the compounds’ mass-by-charge ratio (m/z). The y-axis shows retention time (RT) of the compound on the liquid chromatography preceding mass spectrometry. Colors indicate two clusters of lipids showing intensity differences between humans and two macaque species. Each point represents annotated TAG feature; labels include the cumulative length of the carbon chains and the total number of double bonds of the fatty acid residues.

Additional file 3: Figure S3. Mass spectrometry output showing 76 lipid features computationally annotated as TAGs. The x-axis shows the compounds’ mass-by-charge ratio (m/z). The y-axis shows retention time (RT) of the compound on the liquid chromatography preceding mass spectrometry. Colors indicate two clusters of lipids showing intensity differences between humans and bovid species (cows, yaks, goats). Each point represents annotated TAG feature; labels include the cumulative length of the carbon chains and the total number of double bonds of the fatty acid residues.

Additional file 4: Figure S4. Mass spectrometry output showing 76 lipid features computationally annotated as TAGs. The x-axis shows the compounds’ mass-by-charge ratio (m/z). The y-axis shows retention time (RT) of the compound on the liquid chromatography preceding mass spectrometry. Colors indicate two clusters of lipids showing intensity differences between humans and pigs. Each point represents annotated TAG feature; labels include the cumulative length of the carbon chains and the total number of double bonds of the fatty acid residues.

Additional file 5: Table S1. Data overview. Number of samples used in this study.

Additional file 6: Table S2. Computationally annotated TAG features. The feature annotation was based on m/z values of the intact ions and, therefore, informed about the cumulative length and unsaturation degree of three fatty acid residues.

Additional file 7: Table S3. Signal intensities and fold change for annotated TAG features. Average values for upper-quartile normalized, log2 transformed signal intensities for seven species, with minimal (Min) and maximal (Max) values across species, difference between the minimum and maximum (log2FC) and fold change (FC).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mitina, A., Mazin, P., Vanyushkina, A. et al. Lipidome analysis of milk composition in humans, monkeys, bovids, and pigs. BMC Evol Biol 20, 70 (2020).

Download citation


  • Lipidome
  • Milk
  • Evolution
  • Human
  • Fatty acids