At age 2 years, the mean total milk intake was 2.6 (SD 1.2) servings per day, and mean total dairy intake was 4.3 (SD 1.5) servings per day. Of the 852 children, 452 (53.1%) predominantly drank whole milk, 226 (26.5%) drank 2% milk, and 174 (20.4%) drank 1%/skim milk. At age 3, 222 (26.1%) had a BMI ≥ 85th percentile and 76 children (8.9%) had a BMI ≥95th percentile. Among 113 children at age 2 with BMI ≥ 85th percentile, 81 (72%) had a BMI ≥ the 85th percentile at age 3. At age 3, weight-for-age z-score was 0.49 units (64th percentile, weight-for-age), a change from a mean weight-for-age z-score of 0.21 units (56th percentile, weight-for-age) at age 2. Mean height-forage z-score at age 3 was 0.27 units (58th percentile, height-for-age), a change from a mean height-for-age z-score of 0.68 units (68th percentile, height-for-age) at age 2.
Participant characteristics by type of milk intake are shown in . Among whole milk drinkers, 34% were of non-white race/ethnicity, compared with 20% and 13% of 2% and 1%/skim milk drinkers respectively. Mean birth weight was lowest among whole milk drinkers (3435 g, compared with 3552 g and 3590 g for 2% and skim/1% drinkers respectively). Mean BMI, BMI z-score, and % overweight at both ages 2 and 3 were lowest in the whole milk group, and highest in the 1%/skim milk group. Mean child 3-year height was slightly higher in the 1%/skim milk group (97.9 cm) compared with the 2% (97.3 cm) and whole milk (96.8 cm) groups. Mean energy intake was higher for whole milk drinkers than for the other milk groups. Mean calcium intake was 1069 g in the whole milk group, compared with 1171 g/day in the 1%/skim group. Mean fiber intake was slightly higher in the 1%/skim milk group (13.3 g per day) than in the 2% (12.6 g per day) and whole milk (12.2 g per day) groups. Other demographic characteristics were similar among the three groups.
| Table 1Characteristics of 852 children in Project Viva, by predominant type of milk intake at age 2 years |
shows adiposity outcomes at age 3 years by type and quantity of milk intake at age 2 years, unadjusted for covariates. In unadjusted analyses among whole milk drinkers, mean BMI z-score and the proportion of overweight children appeared to be lower for children who drank 2 or more servings per day, compared with those who drank fewer servings. Among 2% milk drinkers, there was no clear pattern of mean BMI z-score across categories of milk intake, but the proportion of overweight children was higher among children with greater milk intake, ranging from 21.7% for those drinking less than 1 serving per day, to 37.5% among children drinking 5 or more servings. Among 1%/skim milk drinkers, mean BMI z-score was higher across increasing categories of milk intake (0.42 vs. 0.74 units for <1 vs. >5 milk servings per day).
| Table 2Adiposity outcomes at age 3 years by category of milk intake at age 2 years among 852 participants from Project Viva |
Unadjusted and multivariable regression models using milk intake as a continuous variable are presented in and . In unadjusted analyses among whole and 2% milk drinkers, we did not detect an association between 2-year milk intake and 3-year BMI z-score (). After adjustment for covariates including BMI z-score at age 2, whole milk intake at age 2 was associated with a modest decrease in 3-year BMI z-score among whole milk drinkers (−0.09 units per serving, 95%CI −0.16, −0.01), with a similar trend among 2% milk drinkers (−0.08 units, 95%CI −0.17, 0.01). Among 1%/skim milk drinkers, the unadjusted effect estimate for each milk serving at age 2 was a 0.13 unit (95%CI −0.01, 0.26, p=0.06) increment in BMI z-score at age 3; adjustment for covariates including 2-year BMI z-score rendered this association nearly null (0.05 units, 95%CI −0.06, 0.16).
| Table 3Increment in BMI z-score at age 3 years associated with each additional serving per day of milk intake at age 2 years, using separate linear regression models for each milk type, among 852 participants from Project Viva |
| Table 4Odds of incident overweight (BMI ≥ 85th percentile) at age 3 years associated with each additional serving of milk intake at age 2 years, using separate logistic regression models for each milk type, among 645 participants from Project Viva |
To further examine possible reverse causality, we performed analyses restricting the study sample to 656 children who were not overweight at age 2 (BMI 5th to < 85th percentile for age and sex). In these analyses, we found no association between milk intake at age 2 and BMI z-score at age 3 (, Model 3). After adjustment for covariates, the increment in BMI z-score at age 3 years for each daily milk serving at age 2 was −0.05 units (95%CI −0.13, 0.02) for whole milk, −0.08 units (95%CI −0.17, 0.02) for 2% milk, and 0.00 units (95%CI −0.14, 0.14) for 1%/skim milk.
We found no association between milk intake at age 2 and incident overweight at age 3 (), regardless of type of milk intake. After adjustment for covariates (Model 3), the odds of incident overweight at age 3 for each daily milk serving at age 2 were 1.04 (95%CI 0.74, 1.44) for whole milk, 0.91 (95%CI 0.62, 1.34) for 2% milk, and 0.95 (95%CI 0.58, 1.55) for 1%/skim milk.
Neither total milk nor total dairy intake at age 2 was associated with BMI z-score or incident overweight at age 3. After adjustment for covariates (Model 3), the increment in BMI z-score at age 3 years for each daily milk serving at age 2 was −0.05 units (95%CI −0.10, 0.00). The odds of incident overweight at age 3 (Model 3) per daily total milk serving at age 2 was 1.01 (95%CI 0.76, 1.15). After adjustment for covariates (Model 3), the increment in BMI z-score at age 3 per daily serving of all dairy products at age 2 was −0.04 units (95%CI −0.08, 0.01). The odds of incident overweight at age 3 (Model 3) per daily total dairy serving at age 2 was 1.01 (95%CI 0.83, 1.23).