Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
J Hum Lact. Author manuscript; available in PMC 2017 April 26.
Published in final edited form as:
PMCID: PMC5405864

Effect of Promoter Polymorphisms on Cytokine Concentration in Preterm Breast Milk and Subsequent Infant Outcomes



Breast milk concentrations of immune components are variable between women and interleukin (IL) differences may be associated with infant outcomes. Molecular mechanisms for milk variability remain unknown.


The aims were to (1) examine the relationship between maternal IL genotypes and milk concentrations of IL4, IL6, and IL10, (2) describe the trajectories of milk IL change, (3) examine whether maternal IL genotypes predict IL trajectories and/or average weekly IL concentration, and (4) examine if weekly IL levels and/or IL trajectories are associated with infant outcomes.


Milk aliquots were collected from each feeding of mother's own milk and pooled weekly. DNA was extracted from 1 sample of each mother's breast milk whey (n = 64), and single nucleotide polymorphisms (SNPs) of IL genes were genotyped. Milk IL concentrations were measured and trajectory analysis examined IL milk change over time. Multivariate breast milk IL concentration analyses controlled for gestational age and prepregnancy body mass index. Multivariate infant outcome (n = 73) analyses controlled for gestational age and the ratio of human milk to total milk.


Trajectory analysis resulted in linear group shapes, with 2 distinct subgroups in IL6 and 3 subgroups in IL4 and IL10. Trajectory groups trended toward significance with calprotectin, intraventricular hemorrhage, and blood transfusions. Multivariate analyses resulted in trending associations between maternal SNPs and subsequent IL6 and IL10 milk levels. There was a trending relationship between IL milk levels and both fecal calprotectin and intraventricular hemorrhage.


Maternal IL SNPs may affect IL breast milk levels and IL milk levels may be associated with infant outcomes.

Keywords: breastfeeding, breast milk, calprotectin, interleukin, preterm infants, single nucleotide polymorphism, trajectory analysis

Well Established

Serum interleukin (IL) levels are influenced by single nucleotide polymorphisms (SNPs). Breast milk IL concentrations are highly variable between women and may be associated with disparate infant outcomes.

Newly Expressed

This study is the first to show a relationship between maternal IL SNPs and subsequent IL breast milk concentration. Variable milk IL levels were also associated with preterm infant outcomes.


Preterm birth complications are the leading cause of death among children younger than 5 years.1 Preterm infants face increased risks of pneumonia, retinopathy, necrotizing enterocolitis (NEC), and sepsis. Breast milk provides protection against these complications,2-7 and for this reason, the American Academy of Pediatrics recommends that all infants, particularly those weighing less than 1500 grams at birth, receive human milk.8 The robust immunological profile of breast milk may contribute to this protective influence on vulnerable infants.

Preterm infants no longer receive maternal immunological protection in utero; therefore, it is essential that they receive passive immunity through breast milk. Interleukins (ILs) are an integral part of the inflammatory response, and preterm infants are vulnerable to complications due, in part, to an underdeveloped immunological system. Interleukins are present in breast milk, although their concentration is highly variable between women.9 Breast milk from mothers with allergies is higher in IL4 concentration when compared to milk from mothers without allergies.10,11 IL4 is relevant to the preterm population because serum IL4 induces both antibody and IgE production,12 although it is unclear if exogenous milk IL4 has this effect on the breastfed infant. There is some evidence that varying IL4 in breast milk may contribute to the development of allergic dermatitis in healthy term infants.13 Serum IL4 levels have been associated with 2 single nucleotide polymorphisms (SNPs) at the gene's promoter region,14,15 although this relationship has never been explored in breast milk.

IL6, a pleiotropic cytokine, promotes the growth of monoclonal antibodies. IL6 is poorly regulated in preterm infants16; furthermore, higher milk IL6 levels are associated with decreased infant weight gain, percentage fat, and fat mass among healthy breastfed infants,17 although this work was exploratory and included a small sample size (n = 19). Two SNPs, also in the regulatory region, have been implicated in varying IL6 levels,18 although its effect on milk concentration has never been explored.

The role of IL10 in the prevention of inflammatory pathologies, as well as its function as a feedback regulator,19 illustrates its importance to preterm infants who are vulnerable to infection. IL10 concentration in breast milk is higher at 3 months postpartum among mothers who report allergies.20 Variable levels of IL10 in breast milk have been implicated in disparate neonatal outcomes, including immeasurable IL10 in milk fed to infants who developed NEC21 and higher milk levels associated with neonatal jaundice.22 IL10 is influenced by SNPs in the gene's promoter region23-26; however, the effect of SNPs on milk IL10 levels has never been examined.

Breast milk composition is influenced by gestational age at delivery and lactation stage, with higher cytokine milk levels that increase during lactation among preterm milk compared with term milk.27 This pattern of increasing cytokine concentration over time may explain why mother's own milk (MOM) is highly protective against infection, particularly NEC. Trajectory analysis of human milk ILs may illustrate the dynamic nature of breast milk, as this approach allows for estimation of multiple groups within a population, as opposed to a traditional regression that models only 1 mean within the population.28

The association between maternal SNPs and breast milk IL concentration has never been explored, despite evidence that milk IL concentration is variable between women. IL4, IL6, and IL10 are present in milk at highly variable amounts between women. Although limited, exploratory research suggests that milk IL variability may be associated with several infant outcomes. The purpose of this study was to, over the first 3 weeks postpartum and in a population who delivered preterm, (1) examine the relationship between maternal IL genotypes and cross-sectional (weekly) breast milk concentrations of IL4 (rs2243250 -589 T/C, rs2070874 -33 T/C), IL6 (rs1800795 -174 C/G, rs1800796 -572 G/C), and IL10 (rs1800871 -819 C/T, rs1800872 -592 A/C, rs1800896 -1082 G/A), (2) describe the trajectories of breast milk IL concentration change over time, (3) examine whether maternal IL genotypes are associated with breast milk IL trajectories, and (4) examine if weekly IL levels and/or IL trajectories predict infant outcomes.


Study Population

This ancillary study included women (n = 64) who delivered infants (n = 73, including multiples) with a birth weight < 1500 grams who were admitted to a level 3 neonatal intensive care unit (NICU) between 2011 and 2014. We were able to collect genotype data from DNA in prospectively collected breast milk samples over the first 3 weeks postpartum. Mothers with HIV, infants with major congenital anomalies, and moribund infants were excluded. Informed consent was obtained from all subjects and this study was approved by the University of Pittsburgh's Institutional Review Board.

The following variables were collected and available for analysis: maternal age, income, education, ethnicity, race, marital status, and pregnancy history. Medical record data provided information about the labor and delivery of the infant(s). Maternal prepregnancy weight was self-reported, and the height was obtained from the medical chart. The following equation, recommended by the Centers for Disease Control and Prevention, was used to obtain a body mass index (BMI) for each participant: weight (pounds) / [height (inches)]2 × 703. Infant data were obtained from the NICU medical record and included gender, ethnicity, gestational age at birth, birth weight, APGAR scores, ratio of MOM to total milk, and SNAP scores (illness severity). Continuous outcomes included length of stay, weight gain at 6 weeks, days on oxygen, and fecal calprotectin (see below). Categorical outcomes (yes/no) included sepsis, retinopathy of prematurity (ROP), NEC, intraventricular hemorrhage (IVH), blood transfusions, and feeding intolerance.

Breast Milk Collection and Whey Separation

Any volume and source of milk that the infant received (including MOM, donor breast milk, and formula) was recorded from the infant's medical record. Breast milk aliquots from each feeding were collected for up to 3 weeks and these aliquots were pooled weekly. All milk was collected and stored frozen until brought to the laboratory twice weekly. The pooled MOM was centrifuged, defatted, and filtered, and the whey was frozen at −80°C. Although infants received both donor milk and MOM, only MOM was examined and genotyped for this study.

Interleukin Measurement

Interleukin concentrations were measured using a bead-based assay on a MagPix instrument (Luminex) and Millipore kits (Emd Millipore) and are detailed in previous work.29 Interleukin concentrations were measured weekly in all pooled samples of MOM for the first 3 weeks postpartum. Each assay included a standard curve and quality controls, and all samples were done in duplicate. Mean intra-assay coefficients of variation for all cytokines ranged from 3.9% to 6.9%.

Fecal Calprotectin

Fecal calprotectin has been used as a biomarker of inflammation within the preterm population, as calprotectin is an indicator of neutrophil migration.30 Calprotectin levels were measured in weekly stool samples using the PhiCal™ Fecal Calprotectin Immunoassay (Geneva Diagnostics). Quality controls were done for all fecal calprotectin assays. Intra-assay coefficient of variation was 7.7%.


DNA was extracted from maternal breast milk whey using the Qiagen DNA Extraction Mini Kit. Genotype data were collected using TaqMan allele discrimination assays to genotype 7 functional promoter polymorphisms of IL4 (rs2070874, rs2243250), IL6 (rs1800795, rs1800796), and IL10 (rs1800871, rs1800872, rs1800896). We performed TaqMan allelic discrimination with the ABI Prism 7000 Sequence Detection System and SDS software v1.2.3 (Applied Biosystems Inc). Negative controls were included and a portion of the samples was repeated to confirm that they repeatedly discriminated into the same genotype. We also included duplicates and performed independent blinded double calls, and discrepancies were regenotyped. Blinded raw data were reexamined for the SNPs for which Hardy-Weinburg Equilibrium (HWE) was violated to rule out genotyping error.

Statistical Analysis

Preliminary analysis

Statistical analyses were performed using SAS (version 9.4). Univariate outliers were assessed using both frequency tables and graphical methods. Multivariate outliers were assessed using scatterplots. Missing data were assessed for both amount and pattern. Normality was assessed with the Shapiro-Wilks test and graphically at each time point. Data transformations were performed when a regression assumption was compromised. There were 8 sets of multiples, including 1 set of triplets. Milk-specific analyses excluded 1 twin, or 2 triplets, removed randomly, to ensure that each mother was represented only once (n = 64). For infant-specific aims, all infants were included in the analyses (n = 73).

Trajectory modeling

Trajectory modeling (PROC TRAJ) with the censored normal model was used to examine changes in breast milk. When determining the number of trajectory groups, all group orders were set to second order when fitting the maximum number of groups.31 After the number of groups was determined, trajectory shapes were evaluated in a stepwise manner, up to a second order polynomial. Participants were assigned to a trajectory group based on their highest posterior group probability. We compared Bayesian Information Criteria (BIC) between the models to determine the appropriate number of groups and trajectory weights.

Univariate analyses

Univariate analyses were performed for each association, and any relationship with a P value ≤ .20 was included in the multivariate model. Due to a small sample size and the need for ethnic subset analyses, which further decreased sample size, SNP-specific analyses included minor allele absence, with the exception of rs2243250. The minor alleles for rs2243250 are inconsistent between the represented ethnicities in this sample; therefore, genotypes were included in the rs2243250 analysis. The Fisher exact test was used to examine the relationships between (1) SNP and categorical infant outcomes, (2) SNP and IL trajectory group, and (3) IL trajectory group and categorical infant outcomes.

Multivariate analyses

We considered P values ≤ .05 significant, although we also considered trending relationships (P ≤ .10) in our discussion due to (1) the exploratory nature of a pilot study, (2) small sample size, and (3) race-specific subgroup analyses, which further decreased our sample size. Multivariate breast milk IL concentration analyses controlled for both gestational age and prepregnancy BMI. Multivariate infant outcomes analyses controlled for gestational age and ratio of MOM to total milk. Infant outcomes were examined using a multivariate approach only when there was a trending or significant relationship (P ≤ .10) between maternal SNP and IL milk concentration in the multivariate analysis. The association between maternal SNPs and IL trajectory grouping was examined without covariates because this relationship has never been explored. Effect sizes (semipartial omega squared [ωp2]) were calculated for continuous outcomes and odds ratios were calculated for categorical outcomes. Multivariate models that included the total population were examined using both minor allele absence and genotype. Continuous outcomes were examined using multiple linear regression and binary outcomes were examined with multiple logistical regression. The relationship between trajectory group and continuous infant outcomes was assessed by generating multiple contingency tables using IL trajectory group and outcomes. We also performed general linear regression models to measure these associations, since a sample size of at least 100 is ideal to perform trajectory analyses.31


Demographics and HWE

The average maternal age was 28.3 (± 6.8) years, the sample was mostly African American (39.68%) and multiparous (mean of 3 pregnancies ± 2.4), and the average prepregnancy BMI was 27.8 (± 7.3) (Tables 1 and and2).2). There were 3 SNPs not in HWE: rs2070874, rs2243250, and rs1800796 (data not shown). The infants were born at, on average, 28.3 (± 2.4) weeks gestation and weighed 1069.6 (± 216.8) grams. APGAR scores at 1 minute were, on average, 5.9 (± 1.9) and at 5 minutes, 7.4 (± 1.5). SNAP scores averaged 2.4 (± 0.5), they required oxygen for 15.2 (± 21.3) days, and their length of stay in the NICU was 70.5 (± 37.0) days.

Table 1
Maternal Demographics.a
Table 2
Infant Demographics.a

Interleukin Trajectory Modeling

Interleukin levels were natural log transformed to fulfill the normality assumption. The resulting model from the IL4 trajectory model included 3 groups, ordered (1) low linear (34.7%), (2) middle linear (46.5%), and (3) high linear (18.8%) (Figure 1). The IL6 trajectory resulted in 2 groups, ordered (1) low linear (49.7%) and (3) high linear (50.3%) (Figure 2). IL10 resulted in 3 groups, ordered (1) low linear (33.2%), (2) middle linear (46.5%), and (3) high linear (20.3%) (Figure 3).

Figure 1
lnIL4 versus Week Postpartum.
Figure 2
lnIL6 versus Week Postpartum.
Figure 3
lnIL10 versus Week Postpartum.

Maternal Interleukin SNPs and Breast Milk Interleukin Concentration/Interleukin Trajectory

Prepregnancy BMI was inversely associated with IL6 milk concentrations in the rs1800795 model for week 1 (P = .0411) and average IL6 (P = .0126), although this was seen only in African Americans. A similar association among African Americans was observed between BMI and IL6 milk concentration in the rs1800796 model for week 2 (P = .0288) (Table 3). Although there were no significant relationships between rs1800871 minor allele absence and IL10 milk concentration, BMI was inversely associated with milk IL10 levels at week 2 (P = .0350) and average IL10 (P = .0263). Body mass index was also related to IL10 milk concentration when examining minor allele absence of rs1800872 at week 1 (P = .0260), week 2 (P = .0250), and average IL10 (P = .015). Absence of minor allele rs1800896 was nearly significantly associated with milk IL10 levels at week 3 among African Americans (P = .0705, ωp2 = 0.1145) (Table 3; Figure 4). There was a trend toward significance between absence of minor allele for rs1800795 and IL6 milk concentration among Caucasians at week 3 (P = .0966, ωp2 = 0.1002) (Table 3; Figure 5). African Americans have a trend toward significance between rs1800796 minor allele absence and IL6 milk levels at week 2 (P = .0772, ωp2 = 0.0477) (Table 3; Figure 6). After controlling for gestational age and maternal BMI, there were no significant associations between IL4 genotypes and IL4 breast milk concentrations (data not shown).

Figure 4
IL10 Week 3 Milk Levels with and without Minor Allele rs1800896 (P = .0705), African Americans.
Figure 5
IL6 Week 3 Milk Levels with and without Minor Allele rs1800795 (P = .097), Caucasians.
Figure 6
IL6 Week 2 Milk Levels with and without Minor Allele rs1800796 (P = .0772), African Americans.
Table 3
Regression Model for Significant Interleukin Concentration with Genotype/Minor Allele Presence Adjusted for Gestational Age at Delivery and Prepregnancy BMI.a

Interleukin Trajectory, Interleukin Concentration, and Infant Outcomes

After controlling for gestational age and ratio of MOM to total milk received, there was a significant association between IL6 group 1 membership and IVH (odds ratio [OR] = 6.275; 95% confidence interval [CI], 1.076-36.584; P = .0412) (Table 4) and a trend toward significance with fecal calprotectin at week 3 (P = .0822, ωp2 = 0.0376) (Table 5). Infants who received breast milk from IL4 group 2 were 4.16 times more likely to receive a blood transfusion when compared with infants who received breast milk from trajectory group 3 (OR = 4.162; 95% CI, 0.778-22.277; P = .0712) (Table 4). There was a trending relationship between IL6 milk levels at week 2 and calprotectin at week 3 (P = .0978, ωp2 = 0.0832) (Table 6) among African Americans. There were no significant associations between (1) IL SNPs and interleukin trajectory group and (2) IL levels and categorical infant outcomes (data not shown).

Table 4
Regression Model for Categorical Infant Outcomes with Interleukin Trajectory Group Adjusted for Gestational Age and Ratio of Mother's Own Milk to Total Milk Administered (Total Population).
Table 5
Regression Model for Continuous Infant Outcomes with Interleukin Trajectory Group Adjusted for Gestational Age and Ratio of Mother's Own Milk to Total Milk Administered (Total Population).
Table 6
Regression Model for Continuous Infant Outcomes with Interleukin Concentration Adjusted for Gestational Age at Delivery and Ratio of Mother's Own Milk to Total Milk Administered.


This study suggests that maternal SNPs may influence IL milk concentration, and resulting IL milk levels affect continuous NICU outcomes. Our findings are consistent with previous work that demonstrates a functional effect of SNPs on IL concentration for rs1800795,32 rs1800796,18 and rs1800896.25,26 However, this is the first study to reveal this relationship in breast milk. Higher IL6 milk levels were associated with the IL6 -572C allele, and this is consistent with previous studies among acute coronary syndrome patients18 and postoperative coronary artery bypass patients.33 Our findings are shared with previous findings that show lower plasma IL6 levels among those with the IL6 -174G allele34; however, other studies have revealed opposite findings, with higher IL6 levels associated with the G allele among healthy donors,35 people with Eales' disease,36 and patients in septic shock.32 There are also shared findings between IL10 -1082G and subsequent lower IL10 levels, as demonstrated by Qaddourah et al23 in serum, and in vitro by liver transplant recipients.37 In another study, Capasso et al26 demonstrated a different association, with preterm infants who had the -1082G allele experiencing lower IL10 expression. The direction of the SNP's relationship on milk IL levels adds to the already conflicting body of literature. An explanation for these differing findings may be that epigenetic mechanisms also influence IL levels. It is unfortunate that neither this study nor others that examined these SNPs with relation to IL levels included epigenetic mechanisms in their study designs.

Trajectory analyses revealed that infants who received milk from the low-linear IL6 group 1 were more likely to develop IVH than infants who received milk from IL6 group 2. IL6 has been proposed as a strong candidate to modify the risk of perinatal brain injury.38 Preterm infants face deficient cerebral structural support and are vulnerable to brain injury. IL6 crosses the blood–brain barrier39 and has procoagulative properties.40,41 In addition, IL6 activation decreases vitamin K-dependent coagulation factors and subsequent IVH development.40 Neurodevelopmental complications among infants who experience IVH are lessened if they received breast milk.42 When adequate amounts of IL6 are present in breast milk, infants may receive exogenous coagulative protection against IVH, and neuroprotection to infants who had IVH may last into childhood. IL6 has been observed in high amounts in umbilical vein blood among infants who develop IVH,43 although other studies have found no such relationship.44,45 Conflicting infant IL6 levels and how they relate to IVH risk, coupled with our results, indicate that exogenous milk IL6 should be examined closely as potentially mediating IVH development.

Although this study was not powered to detect a relationship between IL milk levels and NEC, we were able to detect trending relationships between ILs and infant fecal calprotectin levels. The risk of NEC among preterm infants peaks between 13 and 21 days,46,47 and our findings suggest that this risk window is reflected in calprotectin levels. Calprotectin is derived mostly from neutrophils and monocytes,48 when white blood cell migration to the intestines activates neutrophils to release this protein.49 Fecal calprotectin has been directly associated with inflammation severity in the small intestine, including NEC.49-52 Fecal calprotectin is also involved in microbiota establishment in preterm infants,53 launching its potential influence on long-term outcomes. The trajectory analysis is consistent with these findings, as lower IL6 milk levels were associated with lower calprotectin levels. This is consistent with the pro-inflammatory role of IL6 in the gut, as illustrated by Hegazy and El-Bedewy,54 who ameliorated colitis in vitro by down-regulating IL6. The proinflammatory properties of IL6 include a role in neutrophil transition to monocyte infiltration during inflammation, suggesting a magnified effect when exogenous IL6 (eg, breast milk) is introduced. In vitro IL6 expression is increased in ileum mucosal tissue within the NEC population,55 and increased serum IL6 levels reflect the clinical severity of NEC.56,57 Higher IL6 levels may predispose the infant to a hyperinflammatory intestinal environment, thereby increasing calprotectin. Cury et al58 reported that increased levels of IL6 result in a loss of bowel homeostasis, which may lead to disease development.

Infants who received milk from IL4 group 2 were more likely to receive a blood transfusion when compared with infants who received milk from IL4 group 3. Although all groups are linear in shape, group 3 is higher, suggesting that those infants receive more IL4 via the breast milk. IL4 has a suppressive role on erythropoiesis.59,60 This is not consistent with our findings, which suggest that infants who receive less IL4 via breast milk are more likely to require a blood transfusion. This finding should be interpreted cautiously, since the overall effect of trajectory grouping on transfusions was not significant.

Among African Americans, the effect of prepregnancy BMI on milk IL10 levels remained consistent over time, and there were no significant differences in BMI between the 3 ethnic subgroups. Abdominal obesity has been associated with lowgrade inflammation and this increases plasma IL6.61 Maternal adiposity has been positively correlated with cord blood IL6 levels in a mixed-race sample,62 but a negative correlation was found in a sample of Mexican mothers,63 suggesting a race-specific influence of BMI on IL production. In the models for which BMI was significantly associated with milk IL levels, it was an inverse association, suggesting that higher BMIs result in lower IL6 milk concentration.

The effect of maternal BMI on milk IL levels suggests that, in addition to SNPs, environmental factors affect milk composition. Immunofactors present in milk are influenced by weight, including lower milk levels of TGF-β2 and sCD14 levels of overweight mothers when compared with normal weight mothers.64 It is estimated that between 50% and 75% of IL10 variability can be exlained by polymorphisms.65 Remaining contributors to IL10 production are unknown, although several studies suggest epigenetic influences.19 Epigenetic contributors to IL milk levels may explain why we did not uncover more SNP/IL relationships. For example, IL6 milk levels are influenced by mastitis,66 preeclampsia,67 cesarean section delivery,68 and maternal smoking.69 In addition, higher IL10 levels are associated with probiotics.70 These findings implore an epigenetic approach to understand other molecular mechanisms for breast milk variability.


There were several limitations to this study, most notably a small sample size. Despite trending relationships, the effect sizes ranged from small to medium, and this warrants future studies with larger sample sizes. To control for population substructure and different allele frequencies between ethnicities, we conducted subgroup analyses, which decreased power. The calprotectin analysis should be interpreted cautiously, since calprotectin levels follow their own trajectory patterns.50 This is relevant, as we examined cross-sectional calprotectin levels and not patterns over time. In addition, given the small sample size and exploratory nature of this pilot study, we did not correct for multiple comparisons, which may increase Type I error.

Much of this study is based on self-reported variables, including BMI, and study participants underreport weight.71 In addition, self-reported ethnicity does not adequately capture inherent biological differences, and ancestral markers are a more reliable way of obtaining biologically relevant information that accounts for admixture.72

Hardy-Weinburg Equilibrium was violated for 3 SNPs (rs2070874, rs2243250, rs1800796), although we report significant findings for only 1 of these (rs1800796). We were able to eliminate genotyping error; therefore, we believe HWE violation was due to a biased sample of women who delivered preterm infants, which enriched for the alleles under investigation. The SNPs included in this study have been implicated in a variety of obstetric complications.73-75

Infants in this study received donor breast milk, which was not included in the analysis. Of the 73 infants in the study, 31 (42.5%) received donor milk during the first 3 weeks of life (comprising, on average, 10.4%, 11.9%, and 8% of total milk received during weeks 1, 2, and 3, respectively). Infants who receive donor milk are exposed to variable amounts of protein and bioactive components.76 Groer et al77 measured ILs in breast milk using MAGPIX and found that donor milk retains some IL4, IL6, and IL10. These IL milk levels were much lower in donor milk compared with MOM, likely due to pasteurization and freeze-thaw cycles. This is relevant, as infants in our study were likely exposed to all three ILs from both MOM and donor milk. This could explain our finding that infants who received milk from IL4 trajectory group 2 were more likely to receive blood transfusions than infants who received milk from IL4 trajectory group 3. It is possible that the infants in group 2 received more donor breast milk, and therefore extra IL4, suppressing erythropoiesis; however, we did not measure donor IL concentration. Last, infant serum IL levels were not collected, and this may have better illustrated the potential relationship between exogenous milk IL and infant outcomes.

Future studies should include larger sample sizes to accommodate for additional confounders, including serum infant IL levels. A racially diverse sample size that is larger than the current study would allow for a more robust examination of the effect of minor alleles on subsequent milk IL levels.


Maternal SNPs may contribute to the immunological profile of breast milk, specifically ILs, which may also affect neonatal outcomes. African American IL milk levels may be influenced by prepregnancy adiposity. In addition, trajectory analyses revealed that IL milk concentration, although linear in shape, is variable in concentration. Future work should examine mechanisms for environmental influences on breast milk.


Funding: The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Targeted Research and Academic Training of Nurses in Genomics (NR00975909), the National Institute of Nursing Research (T32NR009759 and R21NR013094), the International Society of Nurses in Genetics Research Award, the Sigma Theta Tau Research Award (University of Pittsburgh), the Judith A. Erlen Research Award, and the Corrine M. Barnes Award.


Declaration of Conflicting Interests: The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.


1. World Health Organization. Preterm birth. [Accessed November 18, 2015]; Published 2014.
2. Barsam FJ, Borges GS, Severino AB, de Mello LM, da Silva AS, Nunes AA. Factors associated with community-acquired pneumonia in hospitalised children and adolescents aged 6 months to 13 years old. Eur J Pediatr. 2013;172(4):493–499. [PubMed]
3. Manzoni P, Stolfi I, Pedicino R, et al. Human milk feeding prevents retinopathy of prematurity (ROP) in preterm VLBW neonates. Early Hum Dev. 2013;89(suppl 1):S64–S68. [PubMed]
4. Maayan-Metzger A, Avivi S, Schushan-Eisen I, Kuint J. Human milk versus formula feeding among preterm infants: short-term outcomes. Am J Perinatol. 2012;29(2):121–126. [PubMed]
5. Corpeleijn WE, Kouwenhoven SM, Paap MC, et al. Intake of own mother's milk during the first days of life is associated with decreased morbidity and mortality in very low birth weight infants during the first 60 days of life. Neonatology. 2012;102(4):276–281. [PubMed]
6. Schanler RJ. Human milk supplementation for preterm infants. Acta Paediatr Suppl. 2005;94(449):64–67. [PubMed]
7. Cristofalo EA, Schanler RJ, Blanco CL, et al. Randomized trial of exclusive human milk versus preterm formula diets in extremely premature infants. J Pediatr. 2013;163(6):1592–1595. [PubMed]
8. American Academy of Pediatrics. Breastfeeding and the use of human milk. Pediatrics. 2012;115(2):496–506. [PubMed]
9. Lawrence RA, Lawrence RM. Breastfeeding: A Guide for the Medical Profession. Philadelphia, PA: Elsevier Mosby; 2005.
10. Bottcher MF, Jenmalm MC, Garofalo RP, Bjorksten B. Cytokines in breast milk from allergic and nonallergic mothers. Pediatr Res. 2000;47(1):157–162. [PubMed]
11. Marek A, Zagierski M, Liberek A, et al. TGF-beta(1), IL-10 and IL-4 in colostrum of allergic and nonallergic mothers. Acta Biochim Pol. 2009;56(3):411–414. [PubMed]
12. Pene J, Rousset F, Briere F, et al. IgE production by normal human lymphocytes is induced by interleukin 4 and suppressed by interferons gamma and alpha and prostaglandin E2. Proc Natl Acad Sci U S A. 1988;85(18):6880–6884. [PubMed]
13. Ochiai S, Shimojo N, Morita Y, et al. Cytokine biomarker candidates in breast milk associated with the development of atopic dermatitis in 6-month-old infants. Int Arch Allergy Immunol. 2013;160(4):401–408. [PubMed]
14. Nguyen DP, Genc M, Vardhana S, Babula O, Onderdonk A, Witkin SS. Ethnic differences of polymorphisms in cytokine and innate immune system genes in pregnant women. Obstet Gynecol. 2004;104(2):293–300. [PubMed]
15. Cabantous S, Poudiougou B, Oumar AA, et al. Genetic evidence for the aggravation of Plasmodium falciparum malaria by interleukin 4. J Infect Dis. 2009;200(10):1530–1539. [PubMed]
16. Currie AJ, Curtis S, Strunk T, et al. Preterm infants have deficient monocyte and lymphocyte cytokine responses to group B streptococcus. Infect Immun. 2011;79(4):1588–1596. [PMC free article] [PubMed]
17. Fields DA, Demerath EW. Relationship of insulin, glucose, leptin, IL-6 and TNF-alpha in human breast milk with infant growth and body composition. Pediatr Obes. 2012;7(4):304–312. [PMC free article] [PubMed]
18. Malarstig A, Lindahl B, Wallentin L, Siegbahn A. Soluble CD40L levels are regulated by the -3459 A>G polymorphism and predict myocardial infarction and the efficacy of antithrombotic treatment in non-ST elevation acute coronary syndrome. Arterioscler Thromb Vasc Biol. 2006;26(7):1667–1673. [PubMed]
19. Saraiva M, O'Garra A. The regulation of IL-10 production by immune cells. Nat Rev Immunol. 2010;10(3):170–181. [PubMed]
20. Prokesova L, Lodinova-Zadnikova R, Zizka J, et al. Cytokine levels in healthy and allergic mothers and their children during the first year of life. Pediatr Allergy Immunol. 2006;17(3):175–183. [PubMed]
21. Fituch CC, Palkowetz KH, Goldman AS, Schanler RJ. Concentrations of IL-10 in preterm human milk and in milk from mothers of infants with necrotizing enterocolitis. Acta Paediatr. 2004;93(11):1496–1500. [PubMed]
22. Zanardo V, Golin R, Amato M, et al. Cytokines in human colostrum and neonatal jaundice. Pediatr Res. 2007;62(2):191–194. [PubMed]
23. Qaddourah RH, Magdoud K, Saldanha FL, et al. IL-10 gene promoter and intron polymorphisms and changes in IL-10 secretion in women with idiopathic recurrent miscarriage. Hum Reprod. 2014;29(5):1025–1034. [PubMed]
24. Lowe PR, Galley HF, Abdel-Fattah A, Webster NR. Influence of interleukin-10 polymorphisms on interleukin-10 expression and survival in critically ill patients. Crit Care Med. 2003;31(1):34–38. [PubMed]
25. Yilmaz V, Yentur SP, Saruhan-Direskeneli G. IL-12 and IL-10 polymorphisms and their effects on cytokine production. Cytokine. 2005;30(4):188–194. [PubMed]
26. Capasso M, Avvisati RA, Piscopo C, et al. Cytokine gene polymorphisms in Italian preterm infants: association between interleukin-10 -1082 G/A polymorphism and respiratory distress syndrome. Pediatr Res. 2007;61(3):313–317. [PubMed]
27. Collado MC, Santaella M, Mira-Pascual L, et al. Longitudinal study of cytokine expression, lipid profile and neuronal growth factors in human breast milk from term and preterm deliveries. Nutrients. 2015;7(10):8577–8591. [PMC free article] [PubMed]
28. Arrandale V, Koehoorn M, MacNab Y, Kennedy SM. How to use SAS® Proc Traj and SAS® Proc Glimmix in respiratory epidemiology. Published 2006.
29. Groer M, Duffy A, Morse S, et al. Cytokines, chemokines, and growth factors in banked human donor milk for preterm infants. J Hum Lact. 2014;30(3):317–323. [PMC free article] [PubMed]
30. Kapel N, Roman C, Caldari D, et al. Fecal tumor necrosis factor-alpha and calprotectin as differential diagnostic markers for severe diarrhea of small infants. J Pediatr Gastroenterol Nutr. 2005;41(4):396–400. [PubMed]
31. Nagin DS, Tremblay RE. What has been learned from group-based trajectory modeling? Examples from physical aggression and other problem behaviors. Ann Am Acad Pol Soc Sci. 2005;602(1):82–117.
32. Tischendorf JJ, Yagmur E, Scholten D, et al. The interleukin-6 (IL6)-174 G/C promoter genotype is associated with the presence of septic shock and the ex vivo secretion of IL6. Int J Immunogenet. 2007;34(6):413–418. [PubMed]
33. Brull DJ, Montgomery HE, Sanders J, et al. Interleukin-6 gene -174g>c and -572g>c promoter polymorphisms are strong predictors of plasma interleukin-6 levels after coronary artery bypass surgery. Arterioscler Thromb Vasc Biol. 2001;21(9):1458–1463. [PubMed]
34. Jones KG, Brull DJ, Brown LC, et al. Interleukin-6 (IL-6) and the prognosis of abdominal aortic aneurysms. Circulation. 2001;103(18):2260–2265. [PubMed]
35. Fishman D, Faulds G, Jeffery R, et al. The effect of novel polymorphisms in the interleukin-6 (IL-6) gene on IL-6 transcription and plasma IL-6 levels, and an association with systemic-onset juvenile chronic arthritis. J Clin Invest. 1998;102(7):1369–1376. [PMC free article] [PubMed]
36. Sen A, Paine SK, Chowdhury IH, et al. Impact of interleu-kin-6 promoter polymorphism and serum interleukin-6 level on the acute inflammation and neovascularization stages of patients with Eales' disease. Mol Vis. 2011;17:2552–2563. [PMC free article] [PubMed]
37. Warle MC, Farhan A, Metselaar HJ, et al. Are cytokine gene polymorphisms related to in vitro cytokine production profiles? Liver Transpl. 2003;9(2):170–181. [PubMed]
38. Baier RJ. Genetics of perinatal brain injury in the preterm infant. Front Biosci. 2006;11:1371–1387. [PubMed]
39. Banks WA, Kastin AJ, Broadwell RD. Passage of cytokines across the blood-brain barrier. Neuroimmunomodulation. 1995;2(4):241–248. [PubMed]
40. Poralla C, Hertfelder HJ, Oldenburg J, Muller A, Bartmann P, Heep A. Elevated interleukin-6 concentration and alterations of the coagulation system are associated with the development of intraventricular hemorrhage in extremely preterm infants. Neonatology. 2012;102(4):270–275. [PubMed]
41. Poralla C, Traut C, Hertfelder HJ, Oldenburg J, Bartmann P, Heep A. The coagulation system of extremely preterm infants: influence of perinatal risk factors on coagulation. J Perinatol. 2012;32(11):869–873. [PubMed]
42. Gibertoni D, Corvaglia L, Vandini S, et al. Positive effect of human milk feeding during NICU hospitalization on 24 month neurodevelopment of very low birth weight infants: an Italian cohort study. PloS One. 2015;10(1):e0116552. [PMC free article] [PubMed]
43. Kassal R, Anwar M, Kashlan F, Smulian J, Hiatt M, Hegyi T. Umbilical vein interleukin-6 levels in very low birth weight infants developing intraventricular hemorrhage. Brain Dev. 2005;27(7):483–487. [PubMed]
44. Sorokin Y, Romero R, Mele L, et al. Umbilical cord serum interleukin-6, C-reactive protein, and myeloperoxidase concentrations at birth and association with neonatal morbidities and long-term neurodevelopmental outcomes. Am J Perinatol. 2014;31(8):717–726. [PMC free article] [PubMed]
45. Bhandari V, Buhimschi CS, Han CS, et al. Cord blood erythropoietin and interleukin-6 for prediction of intraventricular hemorrhage in the preterm neonate. J Matern Fetal Neonatal Med. 2011;24(5):673–679. [PMC free article] [PubMed]
46. Snyder CL, Gittes GK, Murphy JP, Sharp RJ, Ashcraft KW, Amoury RA. Survival after necrotizing enterocolitis in infants weighing less than 1,000 g: 25 years' experience at a single institution. J Pediatr Surg. 1997;32(3):434–437. [PubMed]
47. Llanos AR, Moss ME, Pinzon MC, Dye T, Sinkin RA, Kendig JW. Epidemiology of neonatal necrotising enterocolitis: a population-based study. Paediatr Perinat Epidemiol. 2002;16(4):342–349. [PubMed]
48. Selmi C, Cantarini L, Kivity S, et al. The 2014 ACR annual meeting: a bird's eye view of autoimmunity in 2015. Autoimmun Rev. 2015;14(7):622–632. [PubMed]
49. Yoon JM, Park JY, Ko KO, Lim JW, Cheon EJ, Kim HJ. Fecal calprotectin concentration in neonatal necrotizing enterocolitis. Korean J Pediatr. 2014;57(8):351–356. [PMC free article] [PubMed]
50. Josefsson S, Bunn SK, Domellof M. Fecal calprotectin in very low birth weight infants. J Pediatr Gastroenterol Nutr. 2007;44(4):407–413. [PubMed]
51. Dabritz J, Jenke A, Wirth S, Foell D. Fecal phagocyte-specific S100A12 for diagnosing necrotizing enterocolitis. J Pediatr. 2012;161(6):1059–1064. [PubMed]
52. Aydemir O, Aydemir C, Sarikabadayi YU, et al. Fecal calprotectin levels are increased in infants with necrotizing enterocolitis. J Matern Fetal Neonatal Med. 2012;25(11):2237–2241. [PubMed]
53. Rouge C, Butel MJ, Piloquet H, et al. Fecal calprotectin excretion in preterm infants during the neonatal period. PloS One. 2010;5(6):e11083. [PMC free article] [PubMed]
54. Hegazy SK, El-Bedewy MM. Effect of probiotics on pro-inflammatory cytokines and NF-kappaB activation in ulcer-ative colitis. World J Gastroenterol. 2010;16(33):4145–4151. [PMC free article] [PubMed]
55. Lu P, Struijs MC, Mei J, et al. Endoplasmic reticulum stress, unfolded protein response and altered T cell differentiation in necrotizing enterocolitis. PloS One. 2013;8(10):e78491. [PMC free article] [PubMed]
56. Goepfert AR, Andrews WW, Carlo W, et al. Umbilical cord plasma interleukin-6 concentrations in preterm infants and risk of neonatal morbidity. Am J Obstet Gynecol. 2004;191(4):1375–1381. [PubMed]
57. Morecroft JA, Spitz L, Hamilton PA, Holmes SJ. Plasma interleukin-6 and tumour necrosis factor levels as predictors of disease severity and outcome in necrotizing enterocolitis. J Pediatr Surg. 1994;29(6):798–800. [PubMed]
58. Cury DB, Mizsputen SJ, Versolato C, et al. Serum calprotectin levels correlate with biochemical and histological markers of disease activity in TNBS colitis. Cell Immunol. 2013;282(1):66–70. [PMC free article] [PubMed]
59. Thawani N, Tam M, Stevenson MM. STAT6-mediated suppression of erythropoiesis in an experimental model of malarial anemia. Haematologica. 2009;94(2):195–204. [PubMed]
60. Sawada K, Sato N, Koike T. Inhibition of GM-CSF production by recombinant human interleukin-4: negative regulator of hematopoiesis. Leuk Lymphoma. 1995;19(1-2):33–42. [PubMed]
61. Yudkin JS, Stehouwer CD, Emeis JJ, Coppack SW. C-reactive protein in healthy subjects: associations with obesity, insulin resistance, and endothelial dysfunction: a potential role for cytokines originating from adipose tissue? Arterioscler Thromb Vasc Biol. 1999;19(4):972–978. [PubMed]
62. Catalano PM, Farrell K, Thomas A, et al. Perinatal risk factors for childhood obesity and metabolic dysregulation. Am J Clin Nutr. 2009;90(5):1303–1313. [PubMed]
63. Vega-Sanchez R, Barajas-Vega HA, Rozada G, Espejel-Nunez A, Beltran-Montoya J, Vadillo-Ortega F. Association between adiposity and inflammatory markers in maternal and fetal blood in a group of Mexican pregnant women. Br J Nutr. 2010;104(12):1735–1739. [PubMed]
64. Collado MC, Laitinen K, Salminen S, Isolauri E. Maternal weight and excessive weight gain during pregnancy modify the immunomodulatory potential of breast milk. Pediatr Res. 2012;72(1):77–85. [PubMed]
65. Riiskjaer M, Nielsen K, Steffensen R, Erikstrup C, Forman A, Kruse C. Association of interleukin-10 promoter polymorphism and endometriosis. Am J Reprod Immunol. 2011;65(1):13–19. [PubMed]
66. Mizuno K, Hatsuno M, Aikawa K, et al. Mastitis is associated with IL-6 levels and milk fat globule size in breast milk. J Hum Lact. 2012;28(4):529–534. [PubMed]
67. Erbagci AB, Cekmen MB, Balat O, Balat A, Aksoy F, Tarakcioglu M. Persistency of high proinflammatory cytokine levels from colostrum to mature milk in preeclampsia. Clin Biochem. 2005;38(8):712–716. [PubMed]
68. Mehta R, Petrova A. Biologically active breast milk proteins in association with very preterm delivery and stage of lactation. J Perinatol. 2011;31(1):58–62. [PubMed]
69. Etem Pişkin I, Nur Karavar H, Araslı M, Ermiş B. Effect of maternal smoking on colostrum and breast milk cytokines. Eur Cytokine Netw. 2012;23(4):187–190. [PubMed]
70. Hoppu U, Isolauri E, Laakso P, Matomaki J, Laitinen K. Probiotics and dietary counselling targeting maternal dietary fat intake modifies breast milk fatty acids and cytokines. Eur J Nutr. 2012;51(2):211–219. [PubMed]
71. Connor Gorber S, Tremblay M, Moher D, Gorber B. A comparison of direct vs. self-report measures for assessing height, weight and body mass index: a systematic review. Obes Rev. 2007;8(4):307–326. [PubMed]
72. Yaeger R, Avila-Bront A, Abdul K, et al. Comparing genetic ancestry and self-described race in African Americans born in the United States and in Africa. Cancer Epidemiol Biomarkers Prev. 2008;17(6):1329–1338. [PMC free article] [PubMed]
73. Engel SA, Olshan AF, Savitz DA, Thorp J, Erichsen HC, Chanock SJ. Risk of small-for-gestational age is associated with common anti-inflammatory cytokine polymorphisms. Epidemiology. 2005;16(4):478–486. [PubMed]
74. Wu W, Clark EA, Stoddard GJ, et al. Effect of interleukin-6 polymorphism on risk of preterm birth within population strata: a meta-analysis. BMC Genet. 2013;14:30. [PMC free article] [PubMed]
75. Cochery-Nouvellon E, Nguyen P, Attaoua R, et al. Interleukin 10 gene promoter polymorphisms in women with pregnancy loss: preferential association with embryonic wastage. Biol Reprod. 2009;80(6):1115–1120. [PubMed]
76. Molinari CE, Casadio YS, Hartmann BT, Arthur PG, Hartmann PE. Longitudinal analysis of protein glycosylation and beta-casein phosphorylation in term and preterm human milk during the first 2 months of lactation. Br J Nutr. 2013;110(1):105–115. [PubMed]
77. Groer M, Duffy A, Morse S, et al. Cytokines, chemokines, and growth factors in banked human donor milk for preterm infants. J Hum Lact. 2011;30(3):317–323. [PMC free article] [PubMed]