Submitted:
12 February 2026
Posted:
13 February 2026
You are already at the latest version
Abstract
Background/Objectives: Human milk composition is shaped by gestational age at de-livery and stage of lactation; however, proteomic differences between milk from mothers of preterm and term infants and their temporal patterns remain incompletely characterised. Methods: This prospective study enrolled 40 lactating mothers: 20 who delivered preterm infants (< 32 weeks’ gestation) and 20 who delivered at term (37–42 weeks). Each provided milk samples during first 10 days postpartum and at the fifth week. Milk serum was analysed using quantitative data-independent acquisition mass spectrometry. Differential protein abundance was assessed separately at each time point; functional annotation was performed using Gene Ontology biological process analysis. Results: Eighty samples were analysed. During early lactation, 10 proteins differed significantly, most showing lower abundance in preterm milk. At week five, 19 pro-teins were differentially abundant, predominantly higher in preterm samples. Im-mune-related proteins constituted the largest functional category at both stages. Im-munoglobulin heavy constant gamma 4 remained consistently downregulated in pre-term milk (fold change −1.6). Ferritin heavy chain (1.5) and HLA class II histocompati-bility antigen gamma chain (1.8) were elevated only early, whereas calprotectin subu-nits S100A8 (5.6) and S100A9 (5.2) were markedly upregulated later. Conclusions: Proteomic differences vary across lactation stages, identifying lactation timing as a critical contextual factor in comparative human milk proteome studies.
Keywords:
human milk
; preterm birth
; milk proteome
; lactation stage
; immune-related proteins
; data-independent acquisition (DIA)
; quantitative proteomics
; preterm vs term infants
1. Introduction
Preterm birth is associated with increased nutritional vulnerability resulting from the immaturity of metabolic, gastrointestinal, and immune systems. Optimal postnatal nutrition is therefore a critical component of early care in preterm infants, with human milk consistently recognized as the preferred source of enteral feeding due to its unique nutritional and bioactive properties [1].
Beyond its macronutrient content, human milk contains a complex and dynamic array of proteins involved in immune regulation, inflammatory signaling, tissue development, and host defense [2,3]. The composition of human milk is not static but changes across the course of lactation, reflecting both physiological maturation and adaptive responses to infant needs. These temporal changes are particularly relevant in the context of preterm delivery, where early postnatal requirements differ substantially from those of term infants [4].
Previous studies have reported compositional differences between milk from mothers of preterm and term infants, especially during early lactation. However, most investigations have focused on selected proteins or predefined immunological components, providing limited insight into coordinated changes within the broader milk proteome. As a result, it remains unclear whether preterm delivery is associated with global quantitative alterations in milk protein composition or with selective, time-dependent modulation of specific functional protein groups [5].
High-throughput proteomic approaches offer the opportunity to comprehensively characterize the milk proteome and to detect subtle but biologically coherent changes in protein abundance that may not be captured by targeted analyses. Quantitative proteomic approaches based on data-independent acquisition enable reproducible profiling of complex biological fluids, including human milk, and are well suited for assessing time-dependent changes across lactation [6,7,8,9]. Despite this potential, comparative proteomic studies that simultaneously address gestational age at delivery and changes across lactation remain scarce. In parallel, recent proteomic and multi-omics studies have demonstrated substantial interindividual variability and temporal dynamics in human milk composition [7], underscoring the need for carefully designed comparative analyses.
Understanding the temporal dynamics of the human milk proteome in mothers of preterm infants is essential for interpreting observed compositional differences within a nutritional framework and for distinguishing early, transient adaptations from more persistent features of milk composition. Such knowledge is particularly relevant in clinical nutrition, where human milk is increasingly viewed as an active, adaptive component of neonatal nutritional support rather than a uniform dietary substrate.
Therefore, the aim of this study was to perform a comparative quantitative proteomic analysis of human milk from mothers of preterm and term infants collected at two predefined stages of lactation, to assess whether proteomic differences are consistent or stage-specific.
2. Materials and Methods
2.1. Study Population
The study prospectively included 40 lactating mothers whose infants were admitted to the Neonatal Intensive Care Unit of the Department of Paediatrics, Jagiellonian University Medical College, Krakow, Poland, between October 2020 and November 2021. Two groups were enrolled: mothers of preterm infants born before 32 weeks of gestation and mothers of term infants born between 37 and 42 weeks of gestation. Exclusion criteria for both groups included maternal diabetes mellitus (gestational or pregestational), neonatal phenylketonuria, and other severe maternal diseases potentially affecting lactation, as assessed by the investigators. Each participating mother provided two paired milk samples, collected at two predefined time points: within the first 10 days of lactation and during the fifth week of lactation.
2.2. Sample Preparation
2.2.1. Milk Collection
All milk samples were collected at noon after a minimum fasting period of four hours. Milk expression was performed using an electric breast pump (Medela Symphony) with disposable Symphony one-day pump sets for a duration of 10 min. From each expressed sample, 5 mL of milk was collected and immediately stored at -80 °C until further analysis.
2.2.2. Milk Serum Separation
To minimize interference from milk fat globules and casein micelles during proteomic analysis, milk serum was isolated prior to protein extraction. Samples were thawed at room temperature and centrifuged at 1,500 × g for 10 min to remove the lipid layer. The resulting supernatant was subsequently ultracentrifuged at 100,000 × g for 90 minutes at 4 °C using a Beckman L-60 ultracentrifuge. The clarified milk serum fraction was collected and used for further proteomic analysis.
2.2.3. Sample Preparation for LC-MS/MS Analysis
Milk samples (20 μl) were mixed with 80 µl of lysis buffer (2.5% SDS, 60 mM DTT in 0.1M Tris-HCl pH 7.6), vortexed, incubated in 95 °C for 5 min and clarified by centrifugation at 14 000 g for 30 min. Before protein digestion the total protein concentration in collected lysates was determined by WF-assay [10]. Next, a volume containing 70 µg of total protein were transferred to Microcon-30kDa centrifugal filter units (Merck, Darmstadt, Germany), denaturated with 8 M urea in 0.1 M Tris-HCl pH 8.5 and digested to peptides with a use of filter-aided sample preparation (FASP) protocol [11]. Briefly, proteins were alkylated with iodoacetamide and cleaved with LysC-trypsin mix (Thermo Scientific, Waltham, MA) with the enzyme to protein ratio 1:50. Digestions were carried out overnight in 50 mM Tris-HCl pH 8.5 at 37 °C. After digestion the peptide yields were determined by WF-assay and the aliquots containing equal amount of total peptides were desalted on 96-Well MiniSpin C18 columns (Harvard Apparatus, Holliston, MA). Samples were then concentrated to a volume of ~ 5 µL and stored at -80 °C. For project-specific spectral libraries preparation equal amount of peptides from all samples were combined and subjected to fractionation protocol. HpH-fractionation on C18 Micro SpinColumns (Harvard Apparatus, Holliston, MA) was performed in 50 mM ammonium formate buffer (pH 10) with 13 consecutive injections of the eluent buffer, comprising 5, 10, 12.5, 15, 17.5, 20, 22.5, 25, 27.5, 30, 35 and 50% acetonitrile in 50 mM ammonium formate buffer (pH 10), collected by centrifugation (300 × g, 2 min) and dried in a speedvac concentrator (Eppendorf, Hamburg, Germany). In this way peptides were distributed across 12 HpH fractions and analyzed by LC-MS/MS in DDA acquisition mode for library generation. Prior the analysis all samples and library peptide fractions were solubilized in 0.1% formic acid in a concentration of 0.5 µg/µl and spiked with the iRT peptide mix (Biognosys, Schlieren, Switzerland) for normalization of the retention time.
2.2.4. Liquid Chromatography – Tandem Mass Spectrometry
Peptides (1 µg) were injected onto a nanoEase M/Z Peptide BEH C18 75 µm i.d. × 25 cm column (Waters, Milford, MA) via a trap column nanoEase M/Z Symmetry C18 180 µm i.d. × 2 cm column (Waters, Milford, MA). For library generation, each peptide fraction was separated using a 98 min 1% to 40% B phase linear gradient (A phase - 0.1% FA; B phase - 80% ACN and 0.1% FA) operating at a flow rate of 300 nL/min on an UltiMate 3000 HPLC system (Thermo Scientific, Waltham, MA) and applied to a TripleTOF 6600+ (Sciex, Framingham, MA) mass spectrometer. The main working nano-electrospray ion source (Optiflow, Sciex, Framingham, MA) parameters were as follows: ion spray voltage 3.2 kV, interface heater temperature (IHT) 200 °C, ion source gas 1 (GS1) 10 and curtain gas (CUR) 25. For DDA acquisition, spectra were collected in full scan mode (350–1400 Da), followed by one hundred CID MS/MS scans of one hundred most intense precursor ions from the preceding survey full scan exceeding 100 cps intensity under dynamic exclusion criteria. Samples analyzed in DIA acquisition mode were separated using a 63 min 1% to 40% B phase linear gradient at a flow rate of 300 nL/min. For DIA acquisition, spectra were collected in full scan mode (400–1250 Da), followed by one hundred DIA MS/MS scans using a variable precursor isolation window approach, with m/z windows ranging from 6 to 90 Da.
2.2.5. Mass Spectrometric Raw Data Analysis, Spectral Library Generation and DIA Quantitation
DDA data were searched against the human UniProt database (release 2021_01_04, 17 056 entries) using the Pulsar search engine implemented in Spectronaut software (Biognosys, Schlieren, Switzerland) [12] with default parameters (± 40 ppm mass tolerance on MS1 and MS2 level, mutated decoy generation method, trypsin enzyme specificity). Deep Learning Assisted iRT Regression was set as iRT reference strategy for RT to iRT calibration with minimum R2 set to 0.8. Peptide, protein and PSM FDR were set to 1%. Library was generated using 3-6 fragment ions per precursor.
Project-specyfic library was be then used to analyze the DIA data in Spectronaut (Biognosys, Schlieren, Switzerland). Data were filtered by 1% FDR on peptide and protein level, while quantitation and interference correction were done on the MS2 level. Protein grouping was performed based on the ID picker algorithm [13]. Protein quantities were calculated by averaging the respective peptide intensities, while the latter were obtained as mean precursor quantities. Data were normalized by global regression strategy, while statistical testing for differential protein abundance was done using t-tests with multiple testing correction after Storey [14]. The Partial Least Squares-Discriminant Analysis (PLS-DA) was done in RStudio 2026.01.0 (Build 392) with the use of mixOmics package [15].
2.3. Statistical Analysis
Categorical variables are presented as numbers and percentages. Continuous variables are expressed as mean and standard deviation (SD) or as median with first and third quartiles (Q1-Q3), as appropriate. Normality of distribution was assessed using the Shapiro-Wilk test.
Differences between groups were analyzed using Student’s t-test for normally distributed variables and the Wilcoxon test for non-normally distributed variables. Paired analyses were applied where appropriate to account for the longitudinal, within-subject design. Pairing refers to sample structure, not time-trend testing. Categorical variables were compared using Pearson’s chi-squared test or Fisher’s exact test when the assumptions of the chi-squared test were not met. Differential protein abundance analyses were conducted separately for each lactation time point. Proteins were considered significantly different between groups if they met the predefined statistical threshold (q-value < 0.05) and an absolute fold change value greater than 1.5, as reported in the corresponding result tables.
3. Results
3.1. Study Population and Maternal Characteristics
Milk samples were collected from 40 lactating mothers of enrolled infants, including 20 mothers of preterm infants (group A) and 20 mothers of term infants (group B). Baseline maternal, pregnancy, and neonatal characteristics are summarized in Table 1.
No significant differences were observed between mothers of preterm and term infants with respect to gravidity, parity, maternal age, height, educational level, prevalence of allergic diseases, smoking status, alcohol use, or previous breastfeeding history. Maternal weight gain during pregnancy was significantly lower in mothers of preterm infants compared with mothers of term infants (median 9.5 vs. 15.0 kg, p = 0.0003). No significant differences were found in maternal maximal weight before delivery or maternal weight one week postpartum.
Pregnancy-related conditions, including thyroid disease, hypertension, anemia, proteinuria, swelling, and bacterial or viral infections, did not differ significantly between groups. Antenatal corticosteroid prophylaxis was significantly more frequent in the preterm group (p = 0.0036).
As expected, gestational age and birth weight were significantly lower in the preterm group (both p < 0.0001). Caesarean section was more common among preterm deliveries (p = 0.0084). No significant differences were observed in newborn sex distribution.
3.2. Differential Milk Proteome in Early Lactation
To facilitate functional interpretation, differentially abundant proteins were annotated using Gene Ontology (GO) terms with the DAVID bioinformatics resource [16,17].
Global proteomic differences between milk from mothers of preterm and term infants at each lactation stage are summarized in Figure 1. Volcano plot analysis illustrates the magnitude and statistical significance of protein abundance differences during early and later lactation, while multivariate PLS-DA based on whole-proteome data demonstrates clear separation between preterm and term milk samples at both time points.
Quantitative proteomic analysis revealed moderate differences in protein abundance between milk samples from mothers of preterm and term infants at two predefined stages of lactation. Volcano plots illustrate differentially abundant proteins during the first 10 days postpartum (A) and at the fifth week of lactation (B). Multivariate analysis based on whole-proteome data using partial least squares–discriminant analysis (PLS-DA) demonstrated clear separation between preterm and term milk samples at both time points (C and D, respectively).
During early lactation, only a limited number of proteins differed between groups, with changes predominantly involving immune-related proteins and modest effect sizes. Specifically, quantitative proteomic analysis of milk samples collected within the first 10 days postpartum identified 10 proteins that differed significantly between mothers of preterm and term infants (Table 2).
Most differences reflected lower protein abundance in preterm milk, with only a limited number of proteins showing higher levels compared with term controls.
Proteins increased in preterm milk during early lactation were predominantly related to iron handling and immune regulation, including ferritin heavy chain and HLA class II histocompatibility antigen gamma chain, both showing approximately 1.5-1.8-fold higher abundance. In contrast, several immunoglobulin-related proteins, including immunoglobulin heavy constant gamma 4 and immunoglobulin lambda variable 3-25, were present at lower levels in preterm milk, with fold changes ranging from approximately 1.5- to 2.3-fold lower abundance.
Overall, the magnitude of early-lactation differences was moderate, and the differentially abundant proteins were mainly annotated to immune system–related biological processes, with additional representation of proteins involved in signaling and digestive pathways.
3.3. Differential Milk Proteome in Later Lactation
By the fifth week of lactation, the number and magnitude of proteomic differences increased, with a clear predominance of higher protein abundance in preterm milk (Table 3). At this time point, the pattern of differences shifted, with most proteins exhibiting higher abundance in preterm milk.
The most pronounced changes were observed for calprotectin subunits S100A8 and S100A9, which showed more than five-fold higher abundance in preterm milk compared with term milk. Additional immune-related proteins, including S100A6 and multiple immunoglobulin variable region components, also demonstrated increased abundance in preterm milk, typically in the range of 1.6- to 2.7-fold. In contrast, only a small number of proteins, including immunoglobulin heavy constant gamma 4, remained lower in preterm milk at this later stage.
As in early lactation, most differentially abundant proteins identified at week five were associated with immune-related biological processes, with a smaller subset linked to signaling and structural functions.
4. Discussion
4.1. Summary of Principal Findings
In this paired proteomic analysis of human milk from mothers delivering preterm and term infants, we identified a limited but biologically coherent set of differentially abundant proteins between groups. The number of differentially abundant proteins differed between predefined lactation stages, with fewer proteins detected at the earlier stage and a higher number at the later stage. Importantly, the direction of differences between preterm and term milk was not uniform across lactation stages, suggesting stage-specific modulation rather than a stable group-defining proteomic signature. Proteins related to immune function constituted the dominant functional category among the differentially abundant proteins at both stages. Several proteins exhibited clear stage-dependent patterns, including IGHG4, ferritin heavy chain, and the calprotectin subunits S100A8 and S100A9. Accordingly, the present analyses were restricted to predefined stage-based comparisons, and the findings should be interpreted within this analytical framework.
4.2. Biological Interpretation of Key Proteins
4.2.1. Immune-Related Proteins and Immunoglobulins
Immunoglobulin-related proteins represented a prominent component of the observed proteomic differences between preterm and term milk, consistent with the established immunological role of human milk [3,5]. Among these, IGHG4 showed a consistent decrease in preterm milk compared with term milk across lactation stages. IgG subclasses are known to differ substantially in effector function and immunological relevance, and IgG4 in particular has been associated with regulatory or tolerogenic immune responses rather than classical pro-inflammatory activity [18]. Therefore, the observed reduction in IGHG4 should not be interpreted as indicating a global impairment of immune protection but rather as a selective alteration within the immunoglobulin repertoire. Marked interindividual variability in immunoglobulin abundance has been repeatedly reported in human milk proteomic studies, even within homogeneous populations [8,18], supporting cautious interpretation of subclass-specific differences. In early lactation, the detection of HLA class II–related proteins aligns with previous observations that antigen presentation–associated components may be more prominent shortly after delivery, potentially reflecting maternal immune signaling rather than direct infant immune activation [4].
4.2.2. Proteins Related to Innate Immunity and Inflammatory Response
Proteins involved in innate immune defense constituted a second major group among the differentially abundant proteins. The calprotectin complex, composed of S100A8 and S100A9, showed stage-dependent differences between preterm and term milk. Calprotectin is well recognized for its antimicrobial properties and its role in nutritional immunity through metal chelation, but it is also widely used as a biomarker of inflammatory activity [19,20]. Calprotectin, a key component of innate immune responses, has been proposed to contribute to mucosal immune protection in early life; however, evidence in the context of human milk remains indirect and largely inferential [21]. However, increased abundance may also reflect maternal inflammatory status rather than an adaptive milk-specific mechanism, and these interpretations are not mutually exclusive. Accordingly, the present findings are best interpreted as indicating differential representation of innate immune–related proteins rather than direct evidence of altered inflammatory signaling.
4.2.3. Iron Metabolism and Adaptive Functions
Ferritin heavy chain emerged as a protein with clear stage-dependent differences between preterm and term milk. Ferritin in human milk has been proposed to play a dual role, acting both as an iron-binding protein and as a potential modulator of oxidative stress [22,23]. Previous studies have reported variability in milk ferritin concentrations across lactation and between individuals, suggesting tight regulation rather than passive diffusion from maternal circulation [22,24,25]. In the context of preterm birth, altered ferritin abundance may reflect adaptive modulation of iron handling in milk rather than differences in total iron supply. Importantly, the present study was not designed to assess iron bioavailability or clinical iron status, and no clinical inferences regarding infant iron nutrition can be drawn from these findings.
4.3. Lactation Stage as Contextual Framework Rather than Primary Driver
The interpretation of proteomic differences in human milk must account for the strong influence of lactation stage on milk composition. Rather than assuming that lactation stage itself drives the observed differences, our findings suggest that differences observed at predefined stages reflect temporal modulation of group-specific patterns. This distinction is important, as the direction and magnitude of differences between preterm and term milk varied between stages, arguing against a single, stable proteomic phenotype. Formal longitudinal analyses within individuals across lactation stages were not performed, and therefore no causal inferences regarding temporal changes can be made. Similar stage-dependent variability has been reported in previous proteomic and compositional studies of human milk, underscoring the importance of contextualizing group comparisons within clearly defined lactation windows [4,9].
4.4. Novelty and Significance of the Study
The primary strength of this study lies in its design and analytical approach rather than in the identification of individual biomarkers. The use of paired sampling, with two samples obtained from the same mother at predefined lactation stages, minimizes interindividual variability that is known to be substantial in human milk proteomics [7]. Quantitative DIA-based proteomics enabled robust comparison at the proteome level rather than reliance on selected candidate proteins. Importantly, comparisons between preterm and term milk were conducted separately for each lactation stage, avoiding implicit assumptions of temporal continuity. Together, these design elements allow for a more nuanced interpretation of proteomic differences than is typically possible in cross-sectional milk studies.
4.5. Limitations
Several limitations of the present study should be acknowledged. The sample size was limited, reflecting the complexity of paired proteomic analyses in human milk. Functional validation of the identified proteins was not performed, and conclusions are restricted to differences in protein abundance. Within-group longitudinal changes across lactation stages were not analyzed, precluding direct assessment of temporal trajectories. Finally, unmeasured maternal and perinatal factors may have contributed to the observed variability. Potential confounders were not adjusted due to sample size.
4.6. Implications and Future Directions
The present findings contribute to a growing body of evidence indicating that differences in human milk composition between preterm and term deliveries are protein-specific and context-dependent rather than uniform. These results underscore the importance of considering lactation stage and individual variability when interpreting milk proteomic data. Future studies combining longitudinal designs with functional assays will be required to clarify the biological significance of stage-dependent protein modulation. At present, the implications of these findings are primarily relevant for research into human milk biology rather than for clinical practice.
5. Conclusions
In conclusion, this comparative proteomic analysis demonstrates that differences between human milk from mothers of preterm and term infants are not static but vary according to the stage of lactation. Rather than revealing a single, stable proteomic phenotype associated with preterm delivery, the findings indicate stage-specific modulation of a limited set of proteins, predominantly related to immune function. Importantly, the direction and magnitude of these differences changed between early and later lactation, underscoring lactation timing as a critical contextual factor in the interpretation of human milk proteomic data. These results suggest that observed compositional differences should be viewed as dynamic and context-dependent rather than as fixed characteristics of preterm milk. While the functional implications of these protein-level differences remain to be established, the present study highlights the importance of temporally resolved, paired study designs in advancing our understanding of human milk biology and cautions against oversimplified comparisons between preterm and term milk across undefined lactation stages.
Author Contributions
Conceptualization, N.M., M.Z, M.S.; methodology, N.M, M.S; software, M.S.; validation, M.S.; formal analysis, W.Z, M.S.; investigation, M.Z., M.S.; resources. N.M, M.S. data curation, M.S.; writing—original draft preparation, N.M.; writing—review and editing, N.M.; visualization, N.M, M.S.; supervision, N.M., P.K.; project administration, N.M.; funding acquisition, N.M. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by Jagiellonian University Medical College, subsidy for research activities (DBS), grant number N41/DBS/000422.
Institutional Review Board Statement
The study was conducted in accordance with the Declaration of Helsinki. The study protocol was approved by the Ethical Committee of the Jagiellonian University Medical College (approval No. KBET/1072.6120.196.2020, 25 June 2020).
Informed Consent Statement
Written informed consent was obtained from all participating women.
Data Availability Statement
Proteomic data generated in this study are available from the corresponding author upon reasonable request. Clinical and demographic data are not publicly available due to ethical and privacy restrictions.
Acknowledgments
We would like to express our gratitude to all members of the medical staff of our department who contributed to the conduct of this study, as well as to all the mothers who participated. We also thank Magdalena Oliwa, clinical dietitian, for her assistance with the technical aspects of human milk sample collection.
Conflicts of Interest
The authors declare no conflicts of interests.
References
- Victora, C.G.; Bahl, R.; Barros, A.J.; França, G.V.; Horton, S.; Krasevec, J.; et al. Breastfeeding in the 21st century: Epidemiology, mechanisms, and lifelong effect. Lancet 2016, 387(10017), 475–490. [Google Scholar] [CrossRef]
- Ballard, O.; Morrow, A.L. Human milk composition: Nutrients and bioactive factors. Pediatr Clin North Am. 2013, 60(1), 49–74. [Google Scholar] [CrossRef]
- Andreas, N.J.; Kampmann, B.; Mehring Le-Doare, K. Human breast milk: A review on its composition and bioactivity. Early Hum Dev. 2015, 91(11), 629–635. [Google Scholar] [CrossRef]
- Castellote, C.; Casillas, R.; Ramírez-Santana, C.; Pérez-Cano, F.J.; Castell, M.; Moretones, M.G.; et al. Premature delivery influences the immunological composition of colostrum and transitional and mature human milk. J Nutr. 2011, 141(6), 1181–1187. [Google Scholar] [CrossRef]
- Lönnerdal, B. Bioactive proteins in human milk: Mechanisms of action. J Pediatr. 2010, 156((2) Suppl, S26–S30. [Google Scholar] [CrossRef]
- Zasada, M.; Suski, M.; Bokiniec, R.; Szwarc-Duma, M.; Borszewska-Kornacka, M.K.; Madej, J.; et al. Comparative two time-point proteome analysis of the plasma from preterm infants with and without bronchopulmonary dysplasia. Ital J Pediatr. 2019, 45(1), 112. [Google Scholar] [CrossRef] [PubMed]
- Dekker, P.M.; Boeren, S.; van Goudoever, J.B.; Vervoort, J.J.M.; Hettinga, K.A. Exploring Human Milk Dynamics: Interindividual Variation in Milk Proteome, Peptidome, and Metabolome. J Proteome Res. 2022, 21(4), 1002–1016. [Google Scholar] [CrossRef] [PubMed]
- Molinari, C.E.; Casadio, Y.S.; Hartmann, B.T.; Livk, A.; Bringans, S.; Arthur, P.G.; et al. Proteome mapping of human skim milk proteins in term and preterm milk. J Proteome Res. 2012, 11(3), 1696–1714. [Google Scholar] [CrossRef] [PubMed]
- Zhang, Y.; Zhang, X.; Mi, L.; Li, C.; Zhang, Y.; Bi, R.; et al. Comparative Proteomic Analysis of Proteins in Breast Milk during Different Lactation Periods. Nutrients. 2022, 14(17), 3648. [Google Scholar] [CrossRef]
- Wiśniewski, J.R.; Gaugaz, F.Z. Fast and sensitive total protein and Peptide assays for proteomic analysis. Anal Chem. 2015, 87(8), 4110–4116. [Google Scholar] [CrossRef]
- Wiśniewski, J.R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat Methods. 2009, 6(5), 359–362. [Google Scholar] [CrossRef] [PubMed]
- Bruderer, R.; Bernhardt, O.M.; Gandhi, T.; Miladinović, S.M.; Cheng, L.Y.; Messner, S.; et al. Extending the limits of quantitative proteome profiling with data-independent acquisition and application to acetaminophen-treated three-dimensional liver microtissues. Mol Cell Proteomics. 2015, 14(5), 1400–1410. [Google Scholar] [CrossRef] [PubMed]
- Zhang, B.; Chambers, M.C.; Tabb, D.L. Proteomic parsimony through bipartite graph analysis improves accuracy and transparency. J Proteome Res. 2007, 6(9), 3549–3557. [Google Scholar] [CrossRef] [PubMed]
- Storey, J.D. A Direct Approach to False Discovery Rates. Journal of the Royal Statistical Society Series B: Statistical Methodology. 2002, 64(3), 479–498. [Google Scholar] [CrossRef]
- Rohart, F.; Gautier, B.; Singh, A.; Le Cao, K.A. mixOmics: An R package for ‘omics feature selection and multiple data integration. PLoS Comput Biol. 2017, 13(11), e1005752. [Google Scholar] [CrossRef]
- Huang da, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4(1), 44–57. [Google Scholar] [CrossRef]
- Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; et al. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 50(W1), W216–W221. [Google Scholar] [CrossRef]
- Vidarsson, G.; Dekkers, G.; Rispens, T. IgG subclasses and allotypes: From structure to effector functions. Front Immunol. 2014, 5, 520. [Google Scholar] [CrossRef]
- Stríz, I.; Trebichavský, I. Calprotectin - a pleiotropic molecule in acute and chronic inflammation. Physiol Res. 2004, 53(3), 245–253. [Google Scholar] [CrossRef]
- Burri, E.; Beglinger, C. Faecal calprotectin -- a useful tool in the management of inflammatory bowel disease. Swiss Med Wkly. 2012, 142, w13557. [Google Scholar] [CrossRef]
- Wang, S.; Song, R.; Wang, Z.; Jing, Z.; Ma, J. S100A8/A9 in Inflammation. Front Immunol. 2018, 9, 1298. [Google Scholar] [CrossRef]
- Lönnerdal, B.; Keen, C.L.; Hurley, L.S. Iron, copper, zinc, and manganese in milk. Annu Rev Nutr. 1981, 1, 149–174. [Google Scholar] [CrossRef]
- Lönnerdal, B. Excess iron intake as a factor in growth, infections, and development of infants and young children. Am J Clin Nutr. 2017, 106 (Suppl 6), 1681s–1687s. [Google Scholar] [CrossRef]
- Dempster, W.S.; Heese, H.D.; Pocock, F.H.; Kirsten, G.; Watermeyer, S. Ferritin levels in human milk. Ann Trop Paediatr. 1986, 6(3), 209–212. [Google Scholar] [CrossRef]
- Siimes, M.A.; Vuori, E.; Kuitunen, P. Breast milk iron--a declining concentration during the course of lactation. Acta Paediatr Scand. 1979, 68(1), 29–31. [Google Scholar] [CrossRef]
Figure 1.
Differential proteomic profiles of human milk.
Table 1.
Maternal, pregnancy, and neonatal characteristics of the study population.
| Mothers of Full-Term Babies (n=20) | Mothers of Preterm Babies (n=20) | p-Value | |
|---|---|---|---|
| Characteristics of mothers | |||
| Gravidity; Me (Q1-Q3) | 1 (1-2) | 2 (1-3) | 0.1645W |
| Parity; Me (Q1-Q3) | 1 (1-2) | 2 (1-2) | 0.2057W |
| Maternal age [years]; mean (SD) | 30.2 (5.5) | 29.4 (6.1) | 0.7346T |
| Level of education; n (%) elementary secondary university vocational |
2 (10%) 9 (45%) 9 (45%) 0 (0%) |
1 (5%) 14 (70%) 4 (20%) 1 (5%) |
0.2356F |
| Maternal height [centimetres]; mean (SD) | 165.6 (6.0) | 162.4 (5.2) | 0.0839T |
| Maternal weight gain during pregnancy [kilograms]; Me (Q1-Q3) | 15.0 (12.6 – 17.0) | 9.5 (7.3 – 11.8) | 0.0003W |
| Maternal maximal weight before labour [kilograms]; Me (Q1-Q3) | 76.0 (60.8 – 85.8) | 70.0 (55.5 – 75.0) | 0.0933W |
| Maternal weight one week after labour [kilograms]; Me (Q1-Q3) | 72 (65 – 84) | 64 (58 – 85) | 0.0983W |
| Previously breastfeeding (whole group; n (%) Breastfeeding a previous child (only if there was a previous child); n (%) |
6 (30%) 6 (86%) |
10 (50%) 10 (91%) |
0.1967P 0.7324P |
| Interval between end of previous lactation [months]; mean (SD) | 74 (45) | 52 (32) | 0.3299T |
| Characteristics of the current pregnancy | |||
| Thyroid disease (any); n (%) | 5 (25%) | 2 (10%) | 0.2119P |
| Diabetes (any); n (%) | 0 (0%) | 0 (0%) | - |
| Heart problems; n (%) | 0 (0%) | 0 (0%) | - |
| Hypertension; n (%) | 3 (15%) | 4 (20%) | 0.6773P |
| Proteinuria; n (%) | 0 (0%) | 2 (10%) | 0.1468P |
| Swelling; n (%) | 2 (10%) | 5 (25%) | 0.2119P |
| Anaemia; n (%) | 4 (20%) | 3 (15%) | 0.6773P |
| Bacterial infection (any); n (%) | 4 (20%) | 4 (20%) | 0.9999P |
| Viral infection (any); n (%) | 6 (30%) | 2 (10%) | 0.1138P |
| Urinary tract infection; n (%) | 3 (15%) | 1 (5%) | 0.2918P |
| Vaginal rectal screening for GBS colonization; n (%) Positive Negative |
6 (32%) 13 (68%) |
1 (7%) 13 (93%) |
0.0897P |
| Smoking; n (%) | 1 (5%) | 2 (10%) | 0.5483P |
| Alcohol use; n (%) | 0 (0%) | 0 (0%) | - |
| Antenatal corticosteroid prophylaxis; n (%) No Yes |
20 (100%) 0 (0%) |
13 (65%) 7 (35%) |
0.0036P |
| Gestational age [weeks]; Me (Q1-Q3) | 39 (38-40) | 30 (29-31) | <0.0001W |
| Birth weight [grams], Me (Q1-Q3) | 3020 (2843 - 3515) | 1400 (1225 - 1575) | <0.0001W |
| Mode of delivery; n (%) Caesarean section Vaginal Vaginal + caesarean section (twin pregnancy) |
11 (55%) 9 (45%) 0 (0%) |
18 (90%) 1 (5%) 1 (5%) |
0.0084F |
| Newborn’s gender; n (%) Female Female + Male (twin pregnancy) Male |
12 (60%) 0 (0%) 8 (40%) |
7 (35%) 1 (5%) 12 (60%) |
0.2049F |
Data are presented as mean ± SD or median (Q1–Q3), as appropriate, or as number (percentage). Between-group comparisons were performed using Student’s t-test or the Wilcoxon test for continuous variables and Pearson’s chi-squared test or Fisher’s exact test for categorical variables, as appropriate. Normality was assessed using the Shapiro–Wilk test.
Table 2.
Differentially abundant proteins in human milk during early lactation.
| UniProt Accession | Protein Name | Preterm vs Term (Fold Change) |
|---|---|---|
| Immunity | ||
| P02794 | Ferritin heavy chain | 1.5 |
| P04233 | HLA class II histocompatibility antigen gamma chain | 1.8 |
| P01717 | Immunoglobulin lambda variable 3-25 | -1.5 |
| P01861 | Immunoglobulin heavy constant gamma 4 | -1.6 |
| Signalling process | ||
| P0DTE7;P0DTE8;P0DUB6 | Alpha-amylase 1B;Alpha-amylase 1C;Alpha-amylase 1A | 1.5 |
| Q13444 |
Disintegrin and metalloproteinase domain-containing protein 15 | 1.8 |
| O00300 | Tumor necrosis factor receptor superfamily member 11B | -2.0 |
| Q6UX06 | Olfactomedin-4 | -2.3 |
| Other biological processes | ||
| P10451 | Osteopontin | 1.7 |
| Q10472 | Polypeptide N-acetylgalactosaminyltransferase 1 | -1.7 |
Differentially abundant proteins identified in milk samples collected within the first 10 days of lactation from mothers of preterm and term infants. Protein abundance differences are expressed as Preterm vs Term (fold change). Fold change values >1 indicate higher protein abundance in milk from mothers of preterm infants, whereas values <1 indicate lower abundance compared with term controls. Only proteins meeting the predefined statistical significance threshold (q-value < 0.05) are shown. Functional grouping was based on Gene Ontology biological process annotations using the DAVID bioinformatics resource; proteins may be involved in multiple biological pathways.
Table 3.
Differentially abundant proteins in human milk during later lactation.
| UniProt Accession | Protein Name | Preterm vs Term (Fold Change) |
|---|---|---|
| Immunity | ||
| P05109 | Protein S100-A8 | 5.6 |
| P06702 | Protein S100-A9 | 5.2 |
| P00738 | Haptoglobin | 1.6 |
| P06703 | Protein S100-A6 | 2.7 |
| P01861 | Immunoglobulin heavy constant gamma 4 | -1.6 |
| P01780 | Immunoglobulin heavy variable 3-7 | 1.8 |
| A0A075B6R9;A0A0C4DH68 | Probable non-functional immunoglobulin kappa variable 2D-24;Immunoglobulin kappa variable 2-24 | 1.7 |
| A0A075B6I0 | Immunoglobulin lambda variable 8-61 | 1.9 |
| Signalling process | ||
| O00592 | Podocalyxin | 2.1 |
| Q96DA0 | Zymogen granule protein 16 homolog B | 1.8 |
| P04114 | Apolipoprotein B-100 | -1.7 |
| P12273 | Prolactin-inducible protein | 1.7 |
| Q06481 | Amyloid-like protein 2 | 1.7 |
| P0DTE7;P0DTE8;P0DUB6 | Alpha-amylase 1B;Alpha-amylase 1C;Alpha-amylase 1A | 1.6 |
| Q9UNW1 | Multiple inositol polyphosphate phosphatase 1 | -1.5 |
| Q6WN34 | Chordin-like protein 2 | -1.6 |
| Other | ||
| P13796 | Plastin-2 | 2.5 |
| P07437 | Tubulin beta chain | 1.7 |
| P16671 | Platelet glycoprotein 4 | 1.7 |
Differentially abundant proteins identified in milk samples collected during the fifth week of lactation from mothers of preterm and term infants. Protein abundance differences are expressed as Preterm vs Term (fold change). Fold change values >1 indicate higher protein abundance in milk from mothers of preterm infants, whereas values <1 indicate lower abundance compared with term controls. Only proteins meeting the predefined statistical significance threshold (q-value < 0.05) are shown.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2026 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
