The objective of the present study was to test whether replacement of 24 ounces of sugar-sweetened soda per day (~80 g sugar/day) with an energy equivalent amount of reduced fat (2%) milk would increase LDL peak particle diameter in a cohort of overweight adolescent boys who were habitual soda drinkers. A two-period randomized crossover study was implemented with evaluation of lipids and lipoproteins, measures of adiposity, and blood pressure after each 3-week intervention period during which soda or milk was consumed on a daily basis in conjunction with the participants’ usual diets. Exploratory outcomes were insulin resistance as assessed by homeostatic model assessment of insulin resistance (HOMA-IR) and inflammation as assessed by C-reactive protein (CRP), as well as serum uric acid and serum transaminases. In addition, measurements were made of a panel of plasma sphingolipids (SPL), including ceramides (Cer), sphingoid bases (SB), and sphingomyelins (SM), that have been linked to a variety of metabolic diseases and inflammatory conditions [ 15 17 ].

Clinical trials have demonstrated deleterious effects of SSBs as well as fructose on lipid and lipoprotein profiles relative to milk and other beverages [ 8 10 ]. A study in free living adults showed that after 6 months, 32 ounces per day of SSBs (~100 g sucrose) resulted in increased fasting plasma triglyceride and total cholesterol levels compared to an equivalent amount of milk, water, or diet beverages [ 10 ]. Another study in adults randomly assigned to six different beverage interventions providing varying sources and amounts of sugars (fructose, glucose, sucrose) showed that daily consumption of SSB with 80 g sucrose or 80 g fructose for 3 weeks resulted in no significant differences in traditional lipid parameters, including plasma total, low-density lipoprotein (LDL), and high-density lipoprotein (HDL) cholesterol, and triglycerides [ 9 ], but there was a decrease in LDL peak particle size due to a shift in the distribution of LDL from larger to smaller species that are associated to a greater extent with cardiovascular disease risk [ 11 ]. In a cross-sectional study in children, total fructose was found to be the only dietary predictor of LDL particle size [ 12 ], and in children with obesity and metabolic syndrome, fructose restriction in exchange for complex carbohydrates resulted in increased LDL size [ 13 ]. Other cardiometabolic parameters, including blood pressure, fat volume, and inflammatory markers, have been shown to improve with replacement of SSBs with milk and other drinks, such as water and non-caloric diet drinks in adults [ 9 10 ]. These data are consistent with findings from the Harvard Pooling Project of Diet and Coronary Disease, which suggest a decrease in risk of coronary events when milk is substituted for SSBs [ 14 ]. It is not known, however, whether replacement of SSBs with milk improves markers of cardiometabolic risk in adolescents in a similar fashion.

Sugar-sweetened beverages (SSBs) are a major source of energy intake among adolescents, averaging approximately 9–10% of total daily kilocalories for children 12–19 years [ 1 ]. Historically, the increase in SSB consumption in the United States has paralleled the rise in rates of obesity [ 2 ] and related co-morbidities, including atherogenic dyslipidemia, hypertension, and insulin resistance in both adults [ 3 4 ] and children [ 5 6 ], suggesting that SSBs are detrimental to cardiometabolic health. Although SSB consumption in adolescents has declined slightly in recent years, absolute levels of intake remain high [ 1 ]. Furthermore, milk consumption among adolescents, which has declined in recent decades, remains below the recommended daily intake [ 1 ], an issue of particular concern given the nutrients that milk provides, including calcium, magnesium, potassium, and vitamin D. Notably, substitution models have indicated that replacing sugar-sweetened beverages with milk is inversely associated with weight gain in children and adolescents [ 7 ].

The primary outcome of the study was LDL peak particle diameter as an index of atherogenic dyslipidemia. Secondary outcomes included other lipid and lipoprotein measurements, measures of adiposity, and blood pressure. Exploratory outcomes included HOMA-IR, hsCRP, serum uric acid, transaminases, and sphingolipids. A sample size of 30 participants in a crossover design was determined to detect a minimum detectable difference of ~1% baseline values for LDL peak particle diameter, based on data from previous measurements in children, at 80% power and α = 0.05 [ 24 ]. This would have been sufficient to detect differences previously reported with high sucrose feeding [ 12 ]. Paired t-test or Wilcoxon signed-ranked test was used to evaluate comparisons between the soda and milk arms, with 2-tailed< 0.05 considered significant.

For sphingolipid analyses, plasma samples (100 μL) were spiked with 10 μL of internal standard mix and subsequently extracted based on protocols described in [ 23 ]. Ceramide/Sphingoid Internal Standard Mixture I (Avanti polar lipids; LM6002) was used as an internal standard mix. The detection system was composed of an Agilent 1290 binary gradient ultra-high pressure chromatography system coupled with an Agilent 6490 triple quadrupole mass spectrometer. Sphingolipid metabolites were resolved on a Zorbax RRHD Eclipse Plus C18 column (2.1 × 50 mm; 1.8 micron) fitted with a pre-column composed of identical matrix. Samples were eluted from the column using a binary gradient composed of mobile phase A (2 mM ammonium formate and 0.2% formic acid in 18 mL water) and mobile phase B (1 mM ammonium formate and 0.2% formic acid in MS-grade methanol) at a flow rate of 1 mL/min. Initial composition was 75% B, and this was increased to 80%, 85%, and 100% at 3, 3.1, and 10 min, respectively. The gradient was maintained at 100% B until 8.5 min and subsequently changed back to initial conditions until the end of the run at 10 min. Analysis was carried out using a multiple-reaction-monitoring (MRM) mode. For all compounds, the general source settings in the positive ionization modes were as follows: Gas temperature 200 °C; gas flow, 14 L min; nebulizer 20 psi; sheath gas temperature 250 °C; sheath gas flow 11 L min; capillary voltage 3000 V; and nozzle voltage; 0 V. The fragmentor voltage of 380 V and a dwell time of 15 ms were used for all mass transitions, and both Q1 and Q3 resolutions were set to nominal mass unit resolution.

Plasma triglycerides, total- and HDL-cholesterol, and glucose were measured by enzymatic endpoint analysis on a clinical chemistry analyzer (LIASYS 330) using methodology previously described [ 22 ]. LDL-cholesterol was calculated using the Friedewald equation. Lipid and lipoprotein cholesterol measurements were standardized through the CDC-NHLBI lipid standardization program. Apolipoproteins B and AI were analyzed by immunoturbidimetric assay using the ITA reagent kit. LDL peak particle diameter and plasma concentrations of very low-density lipoprotein (VLDL), intermediate-density lipoprotein (IDL), LDL, and HDL subclasses were analyzed using ion mobility, as previously described [ 11 ]. Inter-assay variation was reduced by inclusion of two in-house controls in each preparatory process and duplicate analysis. A coefficient of variation of less than 15% was maintained. Fasting insulin was measured by enzyme-linked immunosorbent assay, with two in-house quality control standards. HOMA-IR was calculated as [glucose (mg/dl) × insulin (mU/L)]/405. Serum uric acid, ALT, and AST were measured using a chemistry analyzer (LIASYS 330) and reagent kits from AMS Diagnostics (Brookfield, CT, USA). Quality control was maintained by the inclusion of two multi chemistry controls (AMS Diagnostics, Brookfield, CT, USA) with each run. hsCRP was measured at a clinical laboratory (Quest Diagnostics).

Using standard procedures, blood samples were collected from participants after fasting overnight for 12–14 h. Plasma and serum were separated immediately by centrifugation at 4 °C. Measurements included body weight, height, BMI, BMI z -score, and waist circumference measured twice midway between the iliac crest and bottom of the ribs. Body fat percentage was measured on a bioimpedance scale (Tanita TBF-551). Clinical blood pressure measurements were taken after participants remained in a sitting position for at least 5 min. Blood pressure was measured three times, with the last two measurements averaged, and expressed as absolute value and z -score.

Participants were asked to complete an online 24-h dietary recall (ASA24-Kids, NCI) [ 21 ] weekly to compare dietary intake at baseline (including washout periods), during the milk arm, and during the soda arm. Dietary intake data from completed recalls (ranging from one to three recalls per period) were averaged within each study period. Three participants did not complete any recalls for one or more of the study periods and were thus excluded from the dietary analyses. On average, participants completed 8 ± 1 recalls out of a possible 9 recalls for the entire study. Participants were advised to maintain a constant amount of physical activity throughout the study. They were provided with pedometers to estimate physical activity during the intervention periods. Pedometer counts were collected weekly.

Participants were recruited from local high schools, internet advertisements and social media postings, advertisements in Children’s Hospital publications, and referrals. Major eligibility criteria included being male, aged 13–18 years, Tanner stage ≥ 2, body mass index (BMI) between 85and 99percentile, self-reported intake of 24 oz per day or more of sugar-sweetened soda, no history of diabetes or hypertension, fasting triglyceride <3.39 mmol/L, no nicotine use, and no use of lipid or blood pressure lowering or diabetes medications. Participants were assigned their beverage sequence in randomly determined blocks of 2, 4, 6, and 8 individuals using a uniform random number generator by a statistician. Beverage sequences were kept in sealed envelopes and assigned to the participant by study staff 1–2 days prior to starting the first beverage intervention. Due to the nature of the intervention, participants and clinic staff were not blinded; however, laboratory staff and investigators were blinded until all data had been compiled. Figure 1 displays the details of participant recruitment and enrollment. Thirty participants were randomized and all thirty completed the study.

Participants visited the Cholesterol Research Center weekly to consult with staff, be weighed, pick up beverages, and fill out questionnaires assessing their beverage intake and diet. They were instructed to drink only their assigned beverages and water, and to maintain stable dietary habits throughout both arms of the study. Intake of non-milk dairy products was allowed in both arms of the study, with the requirement that it remained consistent for the duration of the study. Clinical staff texted participants daily to remind them to consume their assigned beverages. Participants and their parents were asked to consume similar foods for 3 days prior to each blood draw to minimize variation caused by other dietary factors.

This was a two-period randomized study in 30 overweight or obese adolescent males who were self-reported habitual consumers of SSBs conducted at the Cholesterol Research Center (Berkeley, CA, USA) from June 2014 to September 2016. As assigned beverages added to their habitual diet, participants consumed soda (24 oz per day; sweetened with high fructose corn syrup) or an energy equivalent amount of 2% milk (22 oz per day) in random order for 3 weeks each, separated by a ≥2-week washout period ( Figure 1 ). Previous studies have shown that a three-week dietary intervention is sufficient to induce changes in plasma lipids and lipoproteins [ 18 19 ] and changes in response to high fructose corn syrup supplementation were seen as early as two weeks in young adults [ 20 ]. Blood samples for laboratory measurements, including plasma lipids and lipoproteins, apolipoproteins B and A1, insulin, high sensitivity (hs) CRP, serum uric acid, alanine aminotransferase (ALT), and aspartate transaminase (AST) were obtained at the end of each 3-week period. During the washout period, participants were not provided any beverages and were allowed to return to their usual consumption of foods and beverages. Participants were required to remain weight-stable for the duration of the study.

Replacement of soda with milk resulted in decreased sugar, carbohydrate, and caffeine intake, and increased protein intake ( Table 2 ). Total fat intake did not differ between the soda and milk arms whereas saturated fat was significantly higher on the milk arm, attributable to the increased consumption of dairy fat. Cholesterol intake was not different between the soda and milk phases of the study. As would be expected with increased milk consumption, intake of calcium, magnesium, potassium and zinc was also higher during the milk versus soda arm. Although objective measures of compliance were not made, these self-reports of nutrient intake, based on weekly dietary recalls, suggest good compliance in consuming the test beverages. Dietary compliance scores were 4.8 ± 0.4 out of 5 for the soda intervention and 4.6 ± 0.6 for the milk intervention.

4. Discussion

Recent public health efforts have focused on reducing consumption of sugar-sweetened beverages in children and adolescents due to strong evidence of their deleterious effects on health. We sought to determine whether replacement of sugar-sweetened soda with reduced fat milk would improve markers of cardiometabolic health in adolescent boys. We chose to study adolescent males because, as a group, they consume the greatest amounts of sugar-sweetened beverages [ 26 ]. Furthermore, there are documented gender differences in insulin sensitivity and plasma lipid changes that occur during adolescence [ 27 ], with adolescent boys possessing a greater likelihood of expressing features of atherogenic dyslipidemia than girls [ 28 ]. We hypothesized that in habitual soda drinkers, the replacement of 24 oz of soda containing ~80 g high fructose corn syrup with an energy equivalent amount of reduced fat milk would increase LDL peak particle diameter. Smaller LDL particle diameter is a marker for atherogenic dyslipidemia [ 29 ] and is known to be modifiable by diet. However, in this short term, randomized dietary crossover trial, we did not observe a significant effect of replacing soda with milk on LDL particle diameter. Other biomarkers of cardiometabolic health, including triglycerides, total cholesterol, LDL-cholesterol, HDL-cholesterol, apolipoproteins B and AI, and lipoprotein subfraction concentrations, HOMA-IR as an indicator of glucose homeostasis, hsCRP as a marker for inflammation, and liver enzymes, were also unaffected by the intervention. However, systolic blood pressure and uric acid concentrations were significantly reduced with replacement of soda by milk.

The lack of change in LDL peak particle diameter with milk versus soda consumption in this study contrasts with previous data documenting decreases in LDL particle size or greater abundance of small, dense LDL in weight-stable individuals consuming low fat, high carbohydrate diets [ 18 24 ] and 80 g/day fructose or 80 g/day sucrose feeding [ 9 ]. The effects of milk or dairy consumption on LDL particle size are less clear, with studies showing neutral, positive, or adverse effects, and the type of dairy food and fat content modulating their effects [ 30 ]. As expected, LDL peak particle diameter in our study population was inversely and significantly correlated with plasma triglyceride and waist circumference as a measure of visceral adiposity (data not shown). However, it is notable that the mean LDL peak particle diameter after both the soda and milk intervention periods was relatively large, as measured by ion mobility, i.e., 22.4 nm, a size categorized as pattern A, which is characterized by a greater abundance of larger and more buoyant lipoprotein particles and associated with lower CVD risk than a predominance of smaller LDL (pattern B) [ 11 31 ]. There were only two individuals per treatment arm who exhibited pattern B, a prevalence of ~6%. While this is comparable to what we have previously reported in normal weight children and adolescents [ 32 ], even following consumption of a very high carbohydrate diet [ 24 ], this finding was unexpected given that study participants were overweight and obese. Two studies have reported a high prevalence (40% to 55%) of pattern B in overweight and obese children [ 33 34 ]. Thus, the teens in our study appeared to be protected from the adverse lipid profile that commonly accompanies overweight. The basis for the observed normolipidemia in our overweight and obese study cohort is not known, but it is possible that the teens in this study, who were not insulin resistant and had normal hsCRP levels, were metabolically fit in spite of being overweight [ 35 ]. Their average body fat percentage was in the range considered as ‘high adiposity’ (>75th percentile for age and sex) [ 36 37 ] and their physical activity, as measured by a pedometer, was not greater than in previous studies of normal and obese adolescents [ 38 39 ]. Their ability to maintain normal lipid levels may be related to intrinsic metabolic machinery that also allows them to efficiently burn fuel, including the sugar consumed from sodas. Previous studies in children and adolescents that showed significant lipid improvements with dietary fructose or sugar restriction focused on individuals with metabolic syndrome [ 40 ] and documented nonalcoholic fatty liver disease [ 41 42 ], suggesting that these effects may only be apparent in populations with higher cardiometabolic risk. Furthermore, a recent meta-analysis suggested that effects of high fructose consumption on lipid levels are only apparent in hypercaloric, and not isocaloric, conditions [ 43 ]. By design, the study participants maintained their weight during the study, which may have contributed to the lack of observed differences. We do not believe that our methods of recruitment were biased towards obese, but healthy, adolescents, although it is possible that there may have been unintentional selection for these individuals.

Despite the lack of significant effects on plasma lipoproteins of changing from SSB to milk, significant decreases in glycosphingolipids (LacCer and GluCer) were observed ( Figure 2 ). GluCer and LacCer are the two most abundant glycosphingolipids and are mostly (~60%) distributed in LDL particles [ 44 45 ]. Pre-clinical studies have shown beneficial effects of glycosphingolipid synthesis in improving insulin sensitivity [ 46 ] and atherosclerosis [ 47 ]. LacCer has also been found to be positively correlated with arterial stiffness among middle-aged men [ 48 ]. Increased GluCer has been shown to impair insulin sensitivity in adipocytes [ 49 ]. These results suggest that favorable shifts in LDL SPL composition may occur even in the absence of significant changes in levels of LDL particles.

51, As mentioned above, systolic blood pressure and uric acid were significantly reduced when milk was consumed in place of soda. These beneficial effects are consistent with results of previous studies that have shown similar responses in adults [ 10 50 ]. The metabolism of fructose has been reported to increase uric acid production, and serum uric acid levels have been shown to increase with sugar consumption [ 20 52 ]. The effects of increased fructose intake on blood pressure may be mediated at least in part by uric acid. Uric acid impairs endothelial function and vascular compliance [ 51 ], and epidemiological studies have reported an approximately 1.5- to 2-fold increased risk of hypertension with increased uric acid concentrations [ 51 ]. In children, a strong linear association between serum uric acid concentrations and systolic blood pressure has been observed [ 53 ]. Furthermore, a previous cohort study in adolescents reported associations between higher SSB intake and increases in serum uric acid levels and systolic blood pressure [ 54 ]. Studies in humans and animal models have demonstrated fructose-induced increases in uric acid levels and blood pressure [ 55 56 ]. Two recent intervention studies in children with metabolic syndrome [ 40 ] and nonalcoholic fatty liver disease [ 41 ] reported reduced diastolic and systolic blood pressure, respectively, with reduced fructose consumption. However, the effects on blood pressure in these studies may also have been mediated by reductions in body weight [ 40 ] and body fat [ 41 ] that occurred concomitantly.

Improvements in blood pressure with milk replacement of soda may also have been due to the beneficial effects of milk consumption. Epidemiological evidence has shown a favorable association of dairy product consumption with blood pressure [ 57 58 ]. A dose-response meta-analysis of prospective cohort studies found a 4% decrease in the incidence of hypertension for every 200 g/day of milk intake (our study provided 671 g/day), while high-fat dairy, total fermented dairy, yogurt, and cheese were not associated with hypertension incidence [ 57 ]. The results of clinical trials have been inconsistent, with a recent meta-analysis of randomized trials showing no effect of dairy on systolic blood pressure in interventions lasting one month to one year [ 59 ]; however, the analyses did not consider the source or type of dairy. In one recent trial, the known blood pressure lowering effects of the Dietary Approaches to Stop Hypertension (DASH) diet were shown to be independent of the saturated fatty acid content (provided as dairy) of the diets [ 22 ]. Beneficial effects of dairy on blood pressure may be due to milk proteins or other nutrients contained in dairy products, such as potassium, calcium, and magnesium, or bioactive peptides. Although there is limited evidence to support such effects, a meta-analysis of clinical trials showed a small but significant effect of milk proteins on reducing blood pressure [ 60 ].

There are several limitations to the present study including the short intervention period, restriction to males, and reliance on self-reports to evaluate dietary compliance. In addition, assessment of habitual physical activity by pedometer does not account for the intensity of activity or for calories burned for a given level of exertion. Given the apparent metabolic fitness of our study participants, it would be worthwhile in future studies to investigate the effects of milk vs. SSBs employing more comprehensive evaluations of energy metabolism, for example, through indirect calorimetry and exercise-challenged respiratory rates.