Literature DB >> 32696704

Effects of Sugar-Sweetened, Artificially Sweetened, and Unsweetened Beverages on Cardiometabolic Risk Factors, Body Composition, and Sweet Taste Preference: A Randomized Controlled Trial.

Cara B Ebbeling1, Henry A Feldman2, Sarah K Steltz1, Nicolle L Quinn2, Lisa M Robinson3, David S Ludwig1.   

Abstract

Background A 2018 American Heart Association science advisory indicated that, pending further research, artificially sweetened beverages (ASBs) may be an appropriate initial replacement for sugar-sweetened beverages (SSBs) during transition to unsweetened beverages (USBs). Methods and Results We randomly assigned 203 adults (121 males, 82 females; 91.6% retention), who habitually consumed SSBs, to 3 groups and delivered free SSBs, ASBs, or USBs to their homes for 12 months. Outcomes included serum triglyceride to high-density lipoprotein cholesterol ratio (primary), body weight, and sweet taste preference (experimental assessment, 0%-18% sucrose solutions). Change in serum triglyceride to high-density lipoprotein cholesterol ratio was not different between groups. Although overall change in weight also was not different between groups, we found effect modification (P=0.006) by central adiposity. Among participants in the highest tertile of baseline trunk fat but not other tertiles, weight gain was greater (P=0.002) for the SSB (4.4±1.0 kg, estimate±SE) compared with ASB (0.5±0.9 kg) or USB (-0.2±0.9 kg) group. Both sweetness threshold (-1.0±0.2% m/v; P=0.005) and favorite concentration (-2.3±0.4% m/v; P<0.0001) decreased in the USB group; neither changed in the SSB group. In the ASB group, sweetness threshold did not change, and favorite concentration decreased (-1.1±0.5% m/v; P=0.02). Pairwise comparison between the ASB and USB groups indicated a difference in sweetness threshold (P=0.015). Conclusions Replacing SSBs with noncaloric beverages for 12 months did not affect serum triglyceride to high-density lipoprotein cholesterol ratio. Among individuals with central adiposity, replacing SSBs with either ASBs or USBs lowered body weight. However, USBs may have the most favorable effect on sweet taste preference. Registration URL: https://www.clinicaltrials.gov; unique identifier: NCT01295671.

Entities:  

Keywords:  beverages; body weight; diet; dyslipidemia; sweet taste preference

Mesh:

Substances:

Year:  2020        PMID: 32696704      PMCID: PMC7792240          DOI: 10.1161/JAHA.119.015668

Source DB:  PubMed          Journal:  J Am Heart Assoc        ISSN: 2047-9980            Impact factor:   5.501


artificially sweetened beverage high‐density lipoprotein cholesterol randomized controlled trial sugar‐sweetened beverage triglyceride unsweetened beverage

Clinical Perspective

What Is New?

In this 12‐month intervention study, replacing sugar‐sweetened beverages (SSBs) with artificially sweetened beverages (ASBs) or unsweetened beverages (USBs) had no effect on cardiometabolic risk factors. In a subgroup analysis, individuals with fat around the midsection lost more weight with replacement of SSBs by either ASBs or USBs, and sweet taste preference decreased more when SSBs were replaced with USBs compared with ASBs. The null findings for cardiometabolic risk factors may relate to compensatory changes in other dietary sources of refined carbohydrate; in comparison with randomized controlled trials, the positive associations in observational studies typically include statistical control for carbohydrate source and other measures of dietary quality.

What Are the Clinical Implications?

Overall dietary quality should be considered in interventions involving SSBs to reduce risk for cardiometabolic disease. Consistent with American Heart Association recommendations, substitution of SSBs with ASBs or USBs may help to control body weight in susceptible individuals (eg, those with central adiposity). As replacements for SSBs, USBs decrease sweet taste preference more compared with ASBs, an effect that may have implications for promoting adherence to low‐sugar diets. With increasing awareness of the adverse health outcomes of consuming added sugars, , particularly sugar‐sweetened beverages (SSBs), much attention has been directed toward the effects of replacing SSBs with noncaloric options. The Dietary Guidelines for Americans 2015–2020 emphasizes the importance of drinking unsweetened beverages (USBs), most notably plain water. Pending further research, a science advisory from the American Heart Association indicated that artificially sweetened beverages (ASBs) may be an appropriate initial replacement for SSBs in adults who are habitual consumers, have a strong sweet taste preference, and consider USBs an undesirable alternative. Data are limited regarding the effects of consuming ASBs compared with USBs as replacements for SSBs. Randomized controlled trials (RCTs) of beverage interventions have compared consumption of SSBs versus a combination of ASBs and USBs, , only ASBs, , , and either ASBs or USBs in just 1 study, , with results showing a beneficial effect of various noncaloric beverages on body weight. Prospective observational studies of ASB consumption have yielded mixed results, with some indicating negative associations with weight gain, , , consistent with the cited RCTs, , , , , , and others indicating positive associations with weight gain, , dyslipidemia, , metabolic syndrome, , , type 2 diabetes mellitus, , , , and cardiovascular disease. , , , These positive associations, if causal, may relate to the metabolic effects of the synthetic sweeteners in ASBs (eg, acesulfame potassium, aspartame, and sucralose) that potently activate chemoreceptors for sweet taste at very low concentrations. Of special concern is the possibility that the high‐intensity sweetness of ASBs may increase the preference for sweet foods, in a manner similar to SSBs, with implications for long‐term food choices and dietary quality. , , The aim of this RCT was to compare effects of consuming SSBs, ASBs, and USBs in adults who habitually consumed SSBs. We hypothesized that, as replacements for SSBs, USBs versus ASBs would produce greater reductions in cardiometabolic risk factors, body weight and adiposity, and sweet taste preference. We also explored whether replacing SSBs with ASBs or USBs would be especially beneficial for individuals with high trunk fat mass, in view of the relationship between central adiposity and insulin secretion.

METHODS

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Study Design

We randomly assigned participants to 3 groups (SSB, ASB, USB) for 12 months. The randomization was stratified by sex, ethnicity/race (non‐Hispanic white versus other), and baseline body mass index (normal weight, 18.5–24.9; overweight, 25–29.9; obese, ≥30). Following randomization, participants received a 12‐month intervention to promote consumption of beverages consistent with group assignment. Study outcomes were assessed at baseline and 12 months. The institutional review board at Boston Children’s Hospital approved the study protocol. Participants provided written informed consent. The study was conducted at Boston Children’s Hospital between February 2011 and November 2015. The protocol history, documented in institutional review board amendments, with no change in the primary outcome, is presented in Table S1.

Participants

We screened adults aged 18 to 40 years who reported consuming at least 1 serving (12 fl oz) per day of SSBs and had a body mass index ranging from 18.5 to 40.0 kg/m2. Exclusion criteria included fasting blood glucose ≥110 mg/dL, physician diagnosis of a major medical or psychiatric illness, chronic use of any medication that could affect ≥1 study outcomes, and heavy smoking (>10 cigarettes/d). Additional exclusion criteria for females included: pregnancy (preceding 12 months) or plans to become pregnant during the study period, lactation (preceding 3 months), and change in hormonal contraceptives (preceding 3 months). Before baseline assessments, we obtained medical clearance from a primary care provider. Demographic information based on participant self‐report included: sex, date of birth, ethnicity (Hispanic or non‐Hispanic), race (white, black, Asian, multiple, or other), and total annual household income. We provided monetary incentive upon completion of assessments at baseline ($25) and 12 months ($275).

Interventions

The interventions comprised home deliveries of free beverages, specific messages pertaining to beverage consumption, and check‐in telephone calls. We asked each participant to select a desired beverage combination, consistent with group assignment, from a menu of options and ordered the desired combination from an online delivery service affiliated with a supermarket chain (Peapod, Stop & Shop, Quincy, MA). Unsweetened beverages included spring water, purified water, and sparkling water (with or without flavoring). Beverages were ordered at a frequency to ensure continuous availability, such that delivered beverages could be replacements for the SSBs consumed at baseline, as determined using a beverage frequency questionnaire adapted from a validated instrument. The intervention messages presented in Table 1 were conveyed to study participants in writing and during biweekly check‐in telephone calls. Core messages regarding consumption of delivered beverages were consistent across groups. Group‐specific messages provided instructions on which beverages to avoid. During the check‐in telephone calls, we reviewed beverage consumption, reiterated intervention messages, and encouraged compliance with intervention protocols. Discussions focused exclusively on beverage consumption, with no attention to other dietary behaviors. The target duration of each call was 10 minutes.
Table 1

Intervention Messages

Core Messages (For All Groups)Group‐Specific Messages

Drink the beverages delivered to your home.

The delivered beverages are a replacement for the sugar‐sweetened beverages that you usually drink. Please drink these beverages in place of the ones that you would have purchased prior to enrolling in the research study.*

Drink the delivered beverages at the same rate (number of servings per day) that you were consuming sugar‐sweetened beverages prior to enrolling in the research study.

Do not share the delivered beverages.

SSB Group

Do not drink artificially sweetened (“diet”) beverages.

ASB Group

Do not drink sugar‐sweetened beverages.

USB Group

Do not drink sugar‐sweetened or artificially sweetened (“diet”) beverages.

ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage.

To establish an individualized point of reference for the target number of delivered beverages to be consumed daily, we assessed baseline intake of SSBs using a beverage frequency questionnaire adapted from a validated instrument.

Intervention Messages Drink the beverages delivered to your home. The delivered beverages are a replacement for the sugar‐sweetened beverages that you usually drink. Please drink these beverages in place of the ones that you would have purchased prior to enrolling in the research study.* Drink the delivered beverages at the same rate (number of servings per day) that you were consuming sugar‐sweetened beverages prior to enrolling in the research study. Do not share the delivered beverages. SSB Group Do not drink artificially sweetened (“diet”) beverages. ASB Group Do not drink sugar‐sweetened beverages. USB Group Do not drink sugar‐sweetened or artificially sweetened (“diet”) beverages. ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage. To establish an individualized point of reference for the target number of delivered beverages to be consumed daily, we assessed baseline intake of SSBs using a beverage frequency questionnaire adapted from a validated instrument.

Intervention Fidelity

Procedures to maintain intervention fidelity included group‐specific scripts for reiterating intervention messages, obtaining beverage orders, and negotiating delivery times; systems for cross‐checking beverage orders with supermarket receipts; guides to provide structure for check‐in calls, with adequate flexibility for addressing situations unique to each participant; protocols for documenting each check‐in call; and regular study team meetings to discuss strategies for promoting adherence without compromising differentiation among groups. We also digitally recorded check‐in calls for periodic review.

Process Evaluation

Implementation of the beverage interventions was evaluated based on numbers of completed beverage deliveries and check‐in calls. Three unannounced telephone interviews (2 weekdays, 1 weekend day) were conducted at baseline and again at 12 months to assess dietary intake and physical activity for the day preceding each call. The interviewer was masked to group assignment. Dietary data were collected using the Nutrition Data System for Research software versions 2010–2014, and final data were generated with version 2015 (Nutrition Coordinating Center, University of Minnesota, Minneapolis). Healthy Eating Index 2015 scores (total and subcomponent) were calculated according to a standardized protocol using Nutrition Data System for Research output and code from the Nutrition Coordinating Center. Physical activity data were collected, and daily metabolic equivalent was calculated as previously described.

Outcomes

Study outcomes were measured during 2 visits at each time point (baseline, 12 months), following a 12‐hour overnight fast. Personnel who assessed study outcomes were masked to group assignment. A secure web‐based application (Research Electronic Data Capture), hosted at Boston Children’s Hospital, was used for data management. The prespecified primary outcome was ratio of serum triglyceride (TG) to high‐density lipoprotein cholesterol (HDL‐C) concentration (TG:HDL‐C). Reducing refined sources of dietary carbohydrate, including sugar, has been shown to decrease TG and increase HDL‐C, , , possibly mediated by attenuated postprandial insulinemia. A lipid profile characterized by hypertriglyceridemia and low HDL‐C correlates with insulin resistance and is a component of the metabolic syndrome. Previous studies indicate that the TG:HDL‐C is an independent predictor of coronary atherosclerosis, risk of myocardial infarction, and cardiovascular events and all‐cause mortality. Other prespecified outcomes included blood levels of low‐density lipoprotein cholesterol (direct determination by enzymatic spectrophotometric assay), high‐sensitivity C‐reactive protein, fibrinogen, uric acid, alanine aminotransferase, glucose, and insulin. Insulin sensitivity and β‐cell function, expressed as percentages of values for a normal reference population, were evaluated by homeostasis model assessment using glucose and insulin data. All biochemical analyses were done in facilities certified by the Clinical Laboratory Improvement Amendments. Resting blood pressure was measured by auscultation. Body weight and height were measured using an electronic scale and a wall‐mounted stadiometer. Whole‐body fat mass (henceforth referred to as fat mass) and trunk fat mass (henceforth referred to as trunk fat) were assessed by dual‐energy X‐ray absorptiometry (Discovery A, Hologic, Inc., Bedford, MA). According to our original proposal, insulin secretion was specified as an effect modifier; however, because of scheduling challenges and the need to streamline study visits, we did not measure insulin secretion. In a previous study, we measured insulin concentration 30 minutes after a 75‐g dose of oral glucose (insulin‐30), as a proxy measure of insulin secretion, and noted a significant association between trunk fat and log insulin‐30 (C.B. Ebbeling, H.A. Feldman, and D.S. Ludwig, unpublished data, 2007: N=97; Pearson r=0.46; P<0.0001). Based on this relationship, we conducted exploratory analyses to evaluate effect modification by trunk fat when comparing effects of consuming SSBs, ASBs, or USBs on body weight and fat mass. These analyses are consistent with the American Heart Association science advisory, which notes that replacing SSBs with noncaloric options may be particularly beneficial for weight control among individuals with central adiposity due to high levels of visceral fat. We evaluated sweet taste preference using 10 samples of solutions ranging in sucrose concentration (%m/v) from 0% to 18%. The samples were made using distilled water, flavored with lemon juice, and served at 4°C. On the morning of the taste testing protocol, personnel in the clinical research center at Boston Children’s Hospital prepared and served a standard breakfast (Table S2). Ninety minutes after breakfast, we instructed participants to sip, taste, swallow, and rate the sweetness of each sample in order of increasing sucrose concentration. Sweetness ratings were obtained using a 10‐cm visual analog scale, with verbal anchors at 0 (not at all sweet) and 10 (extremely sweet). Participants consumed a small cracker and sip of water after rating each solution. At the end of the protocol, participants tasted all of the samples again and selected 1 as an overall favorite concentration.

Adverse Events

An adverse event was defined as any symptom or safety concern requiring medical attention, as reported by a participant.

Statistical Analysis

Baseline characteristics and process data were compared across the three study groups by 1‐way ANOVA for continuous measures and Fisher’s exact test for discrete variables. Longitudinal analysis of anthropometric, dietary, and biochemical measures was conducted by repeated‐measures ANOVA. The independent variables were study group and visit (baseline or 12 months), with an interaction term to test for differential change across groups. All models were adjusted for sex, ethnicity, race, and baseline age. A compound‐symmetric covariance structure was employed for the repeated‐measures model, chosen over alternatives that either failed to account for within‐participant correlation (independence), were equivalent to compound symmetry in this 2‐point design (autoregressive, Toeplitz, spatial power), or required additional parameters unjustified by evidence of nonuniform variance (heterogeneous, unstructured). From parameters of the fitted model, we calculated covariate‐adjusted mean changes for each group, with standard error. The mean changes were compared within group to 0 by t test (mean÷standard error), and across groups by the F‐test for interaction. Variables displaying a sharply skewed distribution were log‐transformed for analysis. For descriptive purposes, in place of the unadjusted mean and standard deviation, we retransformed the mean and quartile boundaries of the log‐transformed distribution and reported them as median and interquartile range. The covariate‐adjusted mean changes and standard errors from analysis of log‐transformed variables (Δ±SE) were retransformed for reporting as percentages (100%×exp(Δ)±100%×exp(Δ)×(exp(SE)−1)). To investigate effect modification, we employed analysis of variance with 12‐month change as the dependent variable. The independent variables were study group and the value of the effect modifier at baseline, expressed as tertiles (irrespective of study group), with an interaction term to test whether the group differences varied by tertile. The model was adjusted for sex, ethnicity, race, and baseline age. To characterize sweet taste preference at each time point, we used nonlinear regression analysis to fit a sigmoid curve to each participant's ratings of the 10 sucrose solutions, and calculated sweetness threshold as the concentration corresponding to 5 cm on the 10‐cm visual analog scale (Figure S1).We calculated each participant’s 12‐month change for this parameter and the designated favorite concentration, and then compared the mean changes across study groups. To minimize the influence of extreme values, we applied robust (outlier‐resistant) regression, an iterative procedure by which data points with greater deviation from the fitted mean are given lower weight on subsequent iterations. Among the available variants of this procedure, we employed M‐estimation, most suitable for a discrete independent variable; the quasi‐triangular bisquare weighting function; and a modified Wald test for comparing study groups. We repeated the analysis with adjustment for sex, ethnicity, race, and baseline age, and corroborated the results with the nonparametric Kruskal–Wallis test. We used SAS software (version 9.4, SAS Institute, Cary, NC) for all computations and specified P<0.05 as indicating statistical significance. In cases where the hypothesis of 3 equal groups was rejected, the principle of closed testing allowed us to conduct pairwise comparisons between groups using P<0.05 as critical value while preserving the family‐wise type I error rate of 5% per outcome variable. We applied the 5% family‐wise type I error rate to each study outcome independently, without adjustment for multiple comparisons, following the rationale of Glantz and Slinker. When designing the RCT, we specified a sample size of 270 participants to detect, with 80% power, a mean difference for change in TG:HDL‐C between groups as small as 24%, based on associations observed in the Framingham Heart Study (and a magnitude of effect deemed to have clear clinical implications). Drs Ebbeling, Feldman, and Ludwig had full access to all the data and take responsibility for its integrity and data analysis.

RESULTS

We randomly assigned 203 participants (121 males and 82 females) to an intervention group (Figure 1). The overall retention rate at 12 months was 91.6%. The 17 dropouts were uniformly distributed among study groups (7/67 SSB group, 7/67 ASB group, 3/69 USB group; P=0.32). Dropouts did not differ significantly from the 186 completers with respect to the baseline characteristics shown in Table 2, with minor exceptions. Compared with completers, dropouts more often declined to report household income (11/17 versus 39/186; P<0.001) or education (9/17 versus 7/186; P<0.001), but these discrepancies did not vary significantly across study groups. Dropouts were on average slightly younger than completers (24.2±4.7 versus 27.1±5.7 years, mean±SD; P=0.04), but this difference likewise did not vary significantly across study groups. All parametric analyses were adjusted for age, further obviating any selection bias from this source.
Figure 1

Participant flow. ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage.

Table 2

Baseline Characteristics of Study Participants by Beverage Group

Beverage Group
SSBASBUSB
N676769
Categorical variables* N (%)
Sex
Male40 (59.7)40 (59.7)41 (59.4)
Female27 (40.3)27 (40.3)28 (40.6)
Ethnicity
Hispanic10 (14.9)7 (10.4)8 (11.6)
Race
White38 (56.7)33 (49.3)32 (46.4)
Black14 (20.9)15 (22.4)11 (15.9)
Asian6 (9.0)6 (9.0)14 (20.3)
Multiple/unknown/other9 (13.4)13 (19.4)12 (17.4)
Annual household income
<$30 00013 (19.4)19 (28.4)20 (29.0)
$30 000–$59 99918 (26.9)16 (23.9)12 (17.4)
$60 000–$89 9998 (11.9)9 (13.4)7 (10.1)
≥$90 00010 (14.9)5 (7.5)16 (23.2)
Not reported/unknown18 (26.9)18 (26.9)14 (20.3)
Education
Some high school1 (1.5)0 (0.0)3 (4.4)
High school or GED5 (7.5)5 (7.5)4 (5.8)
Some college or vocational school17 (25.4)20 (29.9)14 (20.3)
Associate’s degree4 (6.0)4 (6.0)3 (4.4)
Bachelor’s degree19 (28.4)19 (28.4)20 (29.0)
Some graduate school or degree16 (23.9)14 (20.9)19 (27.5)
Not reported/unknown5 (7.5)5 (7.5)6 (8.7)
Continuous variablesMean±SD
Age, y25.9±5.126.7±5.727.9±6.0
Weight, kg75.5±15.676.8±16.777.5±16.1
Height, cm170.8±9.2171.3±9.6170.3±10.5
Body mass index, kg/m2 25.8±4.726.1±5.226.6±4.6

ASB indicates artificially sweetened beverage; GED, General Educational Development; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage.

Percentages may not sum to 100 because of rounding.

Ethnicity and race reported by participants.

Participant flow. ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage. Baseline Characteristics of Study Participants by Beverage Group ASB indicates artificially sweetened beverage; GED, General Educational Development; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage. Percentages may not sum to 100 because of rounding. Ethnicity and race reported by participants.

Process Data

Intervention intensity did not differ among groups (Table S3). Participants retained in the study (N=186) received 7.4±4.1 beverage deliveries and completed 19.2±4.5 of 24 planned check‐in calls (mean±SD). Changes in dietary intake are presented in Table 3. Daily consumption of SSBs increased by 1 serving in the SSB group (P<0.001) and declined to almost 0 in the ASB and USB groups (P<0.001). Reflecting these changes, consumption of added sugars and energy from SSBs increased in the SSB group (P=0.002 for added sugars; P<0.001 for energy from SSB) and declined in the ASB and USB groups (P<0.001). Consumption of ASBs increased by 1.5 servings (12 fl oz/serving) in the ASB group (P<0.001) and did not change in the SSB and USB groups. At 12 months, participants in the ASB group consumed (mean±SD) 185.2±202.1 mg/day of aspartame, 268.35±559.89 mg/day of sucralose, and negligible amounts of other artificial sweeteners (including acesulfame potassium and saccharine). Consumption of USBs increased for the USB group (P<0.001) and remained constant for the SSB and ASB groups.
Table 3

Self‐Reported Dietary Intake and Physical Activity

VariableStudy GroupUnadjusted DataAdjusted Data*
Baseline12 mo Mean Change±SE P Within Group P Between Groups
NMean±SDN Mean±SD
Beverages (12 fl oz/serving)
SSBs, servings/dSSB671.4±1.1602.5±1.81.0±0.2<0.001<0.001
ASB671.5±1.4600.1±0.3−1.4±0.2<0.001
USB691.8±1.3650.2±0.5−1.6±0.2<0.001
ASBs, servings/dSSB670.2±0.4600.1±0.2−0.2±0.10.17<0.001
ASB670.1±0.3601.6±1.41.5±0.1<0.001
USB690.2±0.7650.0±0.3−0.2±0.10.14
USBs, servings/dSSB671.6±1.5601.6±1.5−0.0±0.30.92<0.001
ASB671.8±1.6601.9±1.90.1±0.30.70
USB692.1±1.6654.3±2.12.2±0.3<0.001
Sugar
Total, g/dSSB67115.8±50.760139.0±64.223.7±7.40.002<0.001
ASB67122.4±57.96063.6±37.0−57.6±7.4<0.001
USB69123.2±57.06562.8±34.6−60.3±7.1<0.001
Added, g/dSSB6775.7±44.06098.5±61.723.2±7.40.002<0.001
ASB6780.9±51.56034.1±28.4−46.6±7.4<0.001
USB6986.2±47.36530.0±23.1−56.1±7.2<0.001
Energy
Total, kcal/dSSB672053±638602207±577130±790.10<0.001
ASB672225±676601894±576−312±79<0.001
USB692113±689651874±591−247±760.002
From SSBs, kcal/dSSB67202±16760338±256138±25<0.001<0.001
ASB67207±1906021±43−186±25<0.001
USB69242±164659±31−232±24<0.001
Physical activity
Total physical activity (MET)SSB671.58±0.19601.59±0.210.01±0.030.680.90
ASB671.64±0.21601.65±0.260.01±0.030.80
USB691.56±0.19651.58±0.240.03±0.030.37

ASB indicates artificially sweetened beverage; MET, metabolic equivalent; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage.

Data analyzed using repeated measures analysis of variance, adjusted for covariates (sex, ethnicity, race, age).

Dietary recall data missing for one participant in USB group, who could not be reached for telephone interviews at 12 months.

Self‐Reported Dietary Intake and Physical Activity ASB indicates artificially sweetened beverage; MET, metabolic equivalent; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage. Data analyzed using repeated measures analysis of variance, adjusted for covariates (sex, ethnicity, race, age). Dietary recall data missing for one participant in USB group, who could not be reached for telephone interviews at 12 months. Changes in Healthy Eating Index total score and subcomponent scores are presented in Table S4. Change in Healthy Eating Index total score did not differ between groups (P=0.14) but improved in the USB group (P=0.03) on the basis of exploratory within‐group analyses. Changes in the added sugars score differed between groups (P<0.001), consistent with data presented in Table 3. The between‐group effect for the refined grains score was of borderline significance (P=0.06), attributable to decreased consumption in the SSB group (P=0.05) and no change in the ASB and USB groups.

Study Outcomes

Blood lipids, homeostasis model assessment variables, and body weight and fat mass are presented in Table 4. Change in TG:HDL‐C from baseline to 12 months, the prespecified primary outcome, was not significantly different between groups (P=0.65). Likewise, changes in homeostasis model assessment variables (insulin sensitivity, P=0.38; β‐cell function, P=0.49), body weight (P=0.66), fat mass (P=0.27), and other outcomes (Table S5) did not differ significantly between groups. We observed mean (±SE) differences for change in TG:HDL‐C of –5.4%±7.5 for the ASB versus SSB group and –6.2%±7.2 for the USB versus SSB group. We had only 8% power to detect a group effect of this magnitude with the attained sample size, which was smaller than proposed because of recruitment challenges. Nevertheless, enrollment of additional participants to attain the proposed sample size would not have substantially enhanced power to detect such small differences.
Table 4

Blood Lipids, Homeostasis Model Assessment Variables, and Body Weight and Fat Mass

Variable* Study GroupUnadjusted DataAdjusted Data §
Baseline12 mo Mean Change±SE P Within Group P Between Groups
N Mean±SD Median (IQR) N Mean±SD or Median (IQR)
Blood lipids
TG:HDL‐C (primary outcome)* SSB671.37 (0.87–2.16)601.46 (0.98–2.19)3.2±5.70.560.65
ASB671.30 (0.78–1.98)601.19 (0.86–1.87)−2.4±5.40.65
USB691.51 (1.08–2.51)661.45 (0.99–2.31)−3.2±5.10.52
LDL‐C, mg/dLSSB6799.6±26.860102.5±30.51.2±2.60.650.67
ASB67101.5±28.060104.0±30.21.7±2.60.51
USB69109.9±33.366109.1±29.4−1.3±2.50.61
Homeostasis Model Assessment
Insulin sensitivity, %* SSB6789.1 (69.9–161.5)6084.3 (59.3–143.0)−7.5±6.10.220.38
ASB6681.8 (58.7–125.0)6087.5 (60.4–148.8)4.9±6.90.46
USB6985.2 (55.8–122.9)6587.4 (57.5–126.5)−0.6±6.30.93
β‐cell function, %* SSB67120.6 (89.7–141.9)60127.6 (88.4–153.8)0.7±4.20.860.49
ASB66126.3 (98.6–155.8)60122.1 (90.9–147.5)−5.7±4.00.15
USB69122.3 (99.2–163.9)65117.7 (96.4–154.1)−4.3±3.90.27
Body weight and fat mass
Weight, kgSSB6775.5±15.66078.0±17.41.2±0.60.030.66
ASB6776.8±16.76077.0±17.10.6±0.60.32
USB6977.5±16.16678.3±17.10.7±0.50.22
Whole‐body fat mass, kgSSB6422.0±9.64424.2±11.11.0±0.50.030.27
ASB6522.5±9.84422.8±10.20.1±0.50.81
USB6624.0±8.85223.5±8.20.1±0.40.81

ASB indicates artificially sweetened beverage; LDL‐C, low‐density lipoprotein cholesterol; SSB, sugar‐sweetened beverage; TG:HDL‐C, triglyceride to high‐density lipoprotein cholesterol ratio; and USB, unsweetened beverage.

TG:HDL‐C and homeostasis model assessment log‐transformed for analysis, results retransformed for reporting as described below.

Missing data: Homeostasis model assessment missing for 1 participant in ASB group at baseline (missing glucose) and 1 participant in USB group at 12 months (missing insulin); dual‐energy x‐ray absorptiometry whole body fat mass missing for 8 participants at baseline (including 1 dropout in SSB group) and 44 participants at 12 months (scanned on a replacement scanner that was not adequately calibrated with the original scanner); 2 additional participants did not complete the dual‐energy x‐ray absorptiometry scan at 12 months (1 in SSB group, 1 in USB group).

Unadjusted mean±standard deviation; in the case of variables log‐transformed for analysis, median and interquartile range (25th, 75th percentiles).

Mean change±standard error and P‐values, from repeated‐measures analysis of variance, adjusted for prespecified covariates (sex, ethnicity, race, age). For variables log‐transformed for analysis, adjusted change in mean log value and its standard error (Δ±SE) expressed as percentages: 100%×exp(Δ)±100%×exp(Δ)×(exp(SE)−1).

Blood Lipids, Homeostasis Model Assessment Variables, and Body Weight and Fat Mass ASB indicates artificially sweetened beverage; LDL‐C, low‐density lipoprotein cholesterol; SSB, sugar‐sweetened beverage; TG:HDL‐C, triglyceride to high‐density lipoprotein cholesterol ratio; and USB, unsweetened beverage. TG:HDL‐C and homeostasis model assessment log‐transformed for analysis, results retransformed for reporting as described below. Missing data: Homeostasis model assessment missing for 1 participant in ASB group at baseline (missing glucose) and 1 participant in USB group at 12 months (missing insulin); dual‐energy x‐ray absorptiometry whole body fat mass missing for 8 participants at baseline (including 1 dropout in SSB group) and 44 participants at 12 months (scanned on a replacement scanner that was not adequately calibrated with the original scanner); 2 additional participants did not complete the dual‐energy x‐ray absorptiometry scan at 12 months (1 in SSB group, 1 in USB group). Unadjusted mean±standard deviation; in the case of variables log‐transformed for analysis, median and interquartile range (25th, 75th percentiles). Mean change±standard error and P‐values, from repeated‐measures analysis of variance, adjusted for prespecified covariates (sex, ethnicity, race, age). For variables log‐transformed for analysis, adjusted change in mean log value and its standard error (Δ±SE) expressed as percentages: 100%×exp(Δ)±100%×exp(Δ)×(exp(SE)−1). To explore individual differences in response according central adiposity, we divided the cohort into tertiles of baseline trunk fat. As shown in Figure 2 (and Table S6), trunk fat was a significant effect modifier (group×time×trunk fat tertile) for change in body weight (P=0.006) and fat mass (P=0.004). Among participants with the most trunk fat (tertile 3), body weight and fat mass increased significantly more in the SSB group compared with the other 2 groups (P=0.002). There were no discernible group effects for these outcomes among participants in the lower tertiles of trunk fat.
Figure 2

Effect modification by baseline trunk fat for changes in body weight and fat mass. Each bar indicates 12‐month mean change±standard error, from repeated‐measures analysis of variance, adjusted for prespecified covariates (sex, ethnicity, race, age). Within each tertile of trunk fat, P (bottom) tests for difference in mean change across beverage groups. Interaction P (top) tests for difference in beverage effect across tertiles. ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage.

Effect modification by baseline trunk fat for changes in body weight and fat mass. Each bar indicates 12‐month mean change±standard error, from repeated‐measures analysis of variance, adjusted for prespecified covariates (sex, ethnicity, race, age). Within each tertile of trunk fat, P (bottom) tests for difference in mean change across beverage groups. Interaction P (top) tests for difference in beverage effect across tertiles. ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage. Changes in sweet taste preference are shown in Figure 3. Comparison across the 3 beverage groups by robust regression analysis showed a significant difference in the mean change for both sweetness threshold (P=0.008) and favorite concentration (P=0.006). Sweetness threshold did not change discernibly in the SSB or ASB group but decreased significantly in the USB group (P=0.001), with pairwise comparison indicating a significant difference between the ASB and USB groups (P=0.015). Similarly, favorite concentration did not change discernibly in the SSB group but decreased significantly in the USB group (P<0.0001) and ASB group (P=0.02), with the difference between the ASB and USB groups of borderline statistical significance (P=0.06). Adjustment for covariates left the between‐group difference significant for both measures (P=0.01 and P=0.003, respectively). The nonparametric Kruskal–Wallis test confirmed the group differences in distribution (P=0.03 for sweetness threshold; P=0.009 for favorite concentration; Table S7).
Figure 3

Changes in sweet taste preference. Robust (outlier‐resistant) regression analysis provided 12‐month mean change±standard error and P values testing for difference in mean change between groups (top) and within groups (bottom). ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage.

Changes in sweet taste preference. Robust (outlier‐resistant) regression analysis provided 12‐month mean change±standard error and P values testing for difference in mean change between groups (top) and within groups (bottom). ASB indicates artificially sweetened beverage; SSB, sugar‐sweetened beverage; and USB, unsweetened beverage. Three adverse events were documented during the study (hospitalization for body aches; hospitalization for asthma; headache and mood swings attributable to a toothache). Each of these events was deemed unrelated to participation in the RCT.

DISCUSSION

Our study aimed to address a major controversy of relevance to dietary guidelines for the public: Are artificially sweetened beverages equivalent to unsweetened beverages as replacements for sugar‐sweetened beverages? To address this controversy, we conducted an RCT of well‐differentiated interventions, controlling for intervention intensity. We found no overall group differences for changes in TG:HDL‐C and other prespecified cardiometabolic risk factors. Although body weight and fat mass also did not differ among groups, baseline trunk fat was a significant effect modifier for these outcomes. As such, among individuals with central adiposity, replacing SSBs with either ASBs or USBs had a favorable effect on body weight and fat mass. Overall, USBs were a better replacement than ASBs for decreasing sweet taste preference. Length of the intervention period, level of beverage exposure, and individual susceptibility warrant careful consideration when comparing results to data from prospective observational studies , , , and previous RCTs. , , , , , , , The benefits of replacing consumption of SSBs with ASBs or USBs on cardiometabolic risk factors may require longer periods of study for the general population. In prospective observational studies, significant associations between SSB or ASB consumption and mortality have been observed in several cohorts over long‐term follow‐up periods. , , , For example, SSB consumption was positively associated with all‐cause and cardiovascular disease mortality in the Health Professionals’ Follow‐up Study (28 years of follow‐up) and Nurses’ Health Study (34 years of follow‐up), with mortality increasing by 10% with each additional daily serving of SSB. In addition, ASB consumption was associated with all‐cause and cardiovascular disease mortality in the Nurses’ Health study among women who were consuming at least 4 servings per day. Statistical models of beverage substitutions in the Women’s Health Initiative (mean follow‐up of 8.4 years) and European Prospective Investigation into Cancer and Nutrition (mean follow‐up of 10.8 years) indicate reductions in risk for type 2 diabetes mellitus when replacing SSBs with USBs but no benefit when replacing SSBs with ASBs. The intervention messages in the present RCT focused on replacing SSBs consumed at baseline with provided beverages. While relevant from a public health perspective, these messages may have resulted in more variability in consumption and thus less consistent exposure compared with interventions specifying an absolute daily intake. In a Danish RCT of healthy adults, participants were instructed to consume provided beverages at a rate of 1 L (≈36 fl oz) per day. , , At this level of daily consumption for just 6 months, regular cola compared with aspartame‐sweetened cola or water caused greater increases in triglycerides, total cholesterol, and uric acid, but not HDL‐cholesterol and insulin sensitivity. On average, participants in the present RCT did not consume the daily volumes specified in the Danish RCT, even with the unintended increase in consumption of SSBs (in the SSB group) leading to a total volume of ≈30 fl oz per day on average (≈2.5 servings, rather than ≈1.5 servings reflecting baseline consumption). Certain individuals may be particularly susceptible to the adverse effects of dietary carbohydrate on deposition of fat tissue and weight gain and thus more likely to benefit from replacing consumption of SSBs with noncaloric options. , , , , The present RCT indicates that consumption of SSBs had an adverse effect on body weight and fat mass among individuals with high baseline trunk fat, likely attributable in part to increased consumption in the SSB group. Insulin secretion may be one key component of complex mechanisms underlying susceptibility in that individuals with higher central adiposity are more likely to have high initial insulin secretion in response to sugar consumption. As summarized previously, consumption of high‐glycemic‐load sources of carbohydrate, such as SSBs, may promote weight gain by raising the postprandial ratio of serum insulin to glucagon, resulting in increased hunger and decreased energy expenditure. Dietary changes to reduce glycemic load may have the most pronounced effect among individuals with high trunk fat, in whom the postprandial insulin response to oral glucose may be greatest. , In the present RCT, similar changes in body weight among participants with high trunk fat who consumed ASBs or USBs are consistent with studies indicating that mixed meals containing sucralose or aspartame do not raise postprandial blood glucose or insulin levels to the same extent as those containing sucrose. , The observed decrease in sweetness threshold among participants who consumed USBs is consistent with the findings of Wise et al. In their RCT, a dietary intervention to reduce consumption of sugar for 3 months altered perceived sweet taste intensity, such that puddings and beverages containing specified sucrose concentrations were perceived as more sweet in the intervention compared with control group. In the present RCT, favorite concentration also decreased with USBs, and to a lesser extent ASBs. Change in sweet taste preference, achieved by replacing SSBs with USBs, may provide a mechanism for promoting adherence to prescribed low‐sugar diets in the context of comprehensive behavioral intervention programs. Improved Healthy Eating Index total score for the USB group is consistent with the change in overall dietary quality observed by Hedrick et al with an intervention aimed solely at reducing consumption of SSBs. Group effects must be interpreted in the context of the unintended mean increase in consumption of SSBs. With delivery of free SSBs to their homes, participants in the SSB group exhibited propensity to increase, rather than maintain, baseline levels of consumption. Neither sweetness threshold nor favorite concentration changed for the SSB group, suggesting that increased exposure to SSBs did not have an effect on sweet taste preference. While contributing to differentiation in beverage intake among groups (and thus confidence in the null outcomes), increased intake in the SSB group arguably would not threaten external validity for the positive outcomes, as consumption remained within the prevailing range for 60% of adults aged 20 to 39 years who consume SSBs in the United States. In light of the effects of dietary carbohydrate on TG and HDL‐C observed in previous studies, , the spontaneous decrease in consumption of refined grains for the SSB group may have contributed to the null finding. In a recent epidemiologic study, consumption of SSBs was directly associated with adverse changes in TG and HDL‐C in analyses adjusted for dietary quality. To detect the potential effects of beverage consumption on these variables in RCTs, more attention may be needed toward behavioral strategies for controlling intake of other foods (most notably, sources of refined carbohydrate) that could attenuate the independent effects of beverage consumption in intervention studies. Strengths of this RCT include an intervention targeting a single dietary behavior (beverage consumption), home delivery of beverages to promote differentiation in consumption across study groups, examination of several biomarkers of cardiometabolic risk, evaluation of effect modification by baseline central adiposity, a novel protocol for assessing sweet taste preference, a diverse sample (≈50% nonwhite), and high retention rates across groups. Limitations include limited power to observe small effects because of study design (single site rather than multisite study) and recruitment challenges, inability to mask participants to study group assignments, lack of biomarkers of compliance, and reliance on self‐report to assess dietary intake and physical activity (with high likelihood of inaccurate reporting and possibly differential misreporting among intervention groups related to factors such as social desirability bias). To speculate, participants with a strong desire to be viewed favorably by others (high social desirability bias) may have demonstrated more underreporting of actual energy intake in response to interventions focusing on consumption of noncaloric beverages (ASBs and USBs) compared with SSBs. Also, the study was not designed to compare the effects of different artificial sweeteners on study outcomes. Emerging data indicate that, while some metabolic effects are consistent, others vary depending on choice of artificial sweetener. In conclusion, replacing consumption of SSBs with either ASBs or USBs for 12 months had no effect on cardiometabolic risk factors. Among individuals with central adiposity, replacing SSBs with either ASBs or USBs had a favorable effect on body weight and adiposity, consistent with prior findings. , , , , As stated in the American Heart Association science advisory, replacing SSBs with ASBs may be an appropriate initial recommendation for susceptible adults who habitually consume SSBs and consider USBs an undesirable alternative because of a strong sweet taste preference. However, USBs were a better replacement than ASBs for decreasing sweet taste preference, particularly sweetness threshold, a finding with plausible implications for promoting adherence to prescribed low‐sugar diets. In light of epidemiologic data, , , the benefits of eliminating consumption of SSBs and the differential effects of ASBs and USBs on cardiometabolic risk factors may require longer periods of study for the general population.

Sources of Funding

The RCT was funded by grants from the National Heart, Lung, and Blood Institute (R01HL104215) and National Institute of Diabetes and Digestive and Kidney Diseases (K24DK082730 awarded to Dr Ludwig); the National Center for Research Resources (M01RR02172); the Harvard Catalyst Clinical and Translational Science Center (UL1RR025758), and the New Balance Foundation. The views expressed in this article are those of the authors and do not necessarily represent the official views of the sponsors.

Disclosures

Dr Ludwig received royalties for books on obesity and nutrition that recommend a low‐glycemic‐load diet. The remaining authors have no disclosures to report. Tables S1–S7 Figure S1 Click here for additional data file.
  61 in total

1.  Research electronic data capture (REDCap)--a metadata-driven methodology and workflow process for providing translational research informatics support.

Authors:  Paul A Harris; Robert Taylor; Robert Thielke; Jonathon Payne; Nathaniel Gonzalez; Jose G Conde
Journal:  J Biomed Inform       Date:  2008-09-30       Impact factor: 6.317

2.  Reproducibility and validity of dietary patterns assessed with a food-frequency questionnaire.

Authors:  F B Hu; E Rimm; S A Smith-Warner; D Feskanich; M J Stampfer; A Ascherio; L Sampson; W C Willett
Journal:  Am J Clin Nutr       Date:  1999-02       Impact factor: 7.045

3.  Glycemic effect of a single high oral dose of the novel sweetener sucralose in patients with diabetes.

Authors:  N H Mezitis; C A Maggio; P Koch; A Quddoos; D B Allison; F X Pi-Sunyer
Journal:  Diabetes Care       Date:  1996-09       Impact factor: 19.112

4.  Long-Term Consumption of Sugar-Sweetened and Artificially Sweetened Beverages and Risk of Mortality in US Adults.

Authors:  Vasanti S Malik; Yanping Li; An Pan; Lawrence De Koning; Eva Schernhammer; Walter C Willett; Frank B Hu
Journal:  Circulation       Date:  2019-04-30       Impact factor: 29.690

5.  Low-calorie sweetener consumption is increasing in the United States.

Authors:  Allison C Sylvetsky; Jean A Welsh; Rebecca J Brown; Miriam B Vos
Journal:  Am J Clin Nutr       Date:  2012-08-01       Impact factor: 7.045

Review 6.  The Carbohydrate-Insulin Model of Obesity: Beyond "Calories In, Calories Out".

Authors:  David S Ludwig; Cara B Ebbeling
Journal:  JAMA Intern Med       Date:  2018-08-01       Impact factor: 21.873

Review 7.  Measurement of Blood Pressure in Humans: A Scientific Statement From the American Heart Association.

Authors:  Paul Muntner; Daichi Shimbo; Robert M Carey; Jeanne B Charleston; Trudy Gaillard; Sanjay Misra; Martin G Myers; Gbenga Ogedegbe; Joseph E Schwartz; Raymond R Townsend; Elaine M Urbina; Anthony J Viera; William B White; Jackson T Wright
Journal:  Hypertension       Date:  2019-05       Impact factor: 10.190

8.  Reduced dietary intake of simple sugars alters perceived sweet taste intensity but not perceived pleasantness.

Authors:  Paul M Wise; Laura Nattress; Linda J Flammer; Gary K Beauchamp
Journal:  Am J Clin Nutr       Date:  2015-11-25       Impact factor: 7.045

9.  Effects of a low-glycemic load vs low-fat diet in obese young adults: a randomized trial.

Authors:  Cara B Ebbeling; Michael M Leidig; Henry A Feldman; Margaret M Lovesky; David S Ludwig
Journal:  JAMA       Date:  2007-05-16       Impact factor: 56.272

10.  Effects of Sugar-Sweetened, Artificially Sweetened, and Unsweetened Beverages on Cardiometabolic Risk Factors, Body Composition, and Sweet Taste Preference: A Randomized Controlled Trial.

Authors:  Cara B Ebbeling; Henry A Feldman; Sarah K Steltz; Nicolle L Quinn; Lisa M Robinson; David S Ludwig
Journal:  J Am Heart Assoc       Date:  2020-07-22       Impact factor: 5.501

View more
  10 in total

Review 1.  Molecular aspects of fructose metabolism and metabolic disease.

Authors:  Mark A Herman; Morris J Birnbaum
Journal:  Cell Metab       Date:  2021-10-06       Impact factor: 27.287

Review 2.  High versus low-added sugar consumption for the primary prevention of cardiovascular disease.

Authors:  Sara Bergwall; Anna Johansson; Emily Sonestedt; Stefan Acosta
Journal:  Cochrane Database Syst Rev       Date:  2022-01-05

3.  Different Food Sources of Fructose-Containing Sugars and Fasting Blood Uric Acid Levels: A Systematic Review and Meta-Analysis of Controlled Feeding Trials.

Authors:  Sabrina Ayoub-Charette; Laura Chiavaroli; Qi Liu; Tauseef Ahmad Khan; Andreea Zurbau; Fei Au-Yeung; Annette Cheung; Amna Ahmed; Danielle Lee; Vivian L Choo; Sonia Blanco Mejia; Russell J de Souza; Thomas Ms Wolever; Lawrence A Leiter; Cyril Wc Kendall; David Ja Jenkins; John L Sievenpiper
Journal:  J Nutr       Date:  2021-08-07       Impact factor: 4.687

Review 4.  The effects of low-calorie sweeteners on energy intake and body weight: a systematic review and meta-analyses of sustained intervention studies.

Authors:  Peter J Rogers; Katherine M Appleton
Journal:  Int J Obes (Lond)       Date:  2020-11-09       Impact factor: 5.551

5.  Interactions of Non-Nutritive Artificial Sweeteners with the Microbiome in Metabolic Syndrome.

Authors:  Valerie Harrington; Lilian Lau; Alexander Crits-Christoph; Jotham Suez
Journal:  Immunometabolism       Date:  2022-04-18

6.  Metabolic Syndrome and the Risk of Preclinical Heart Failure: Insights after 17 Years of Follow-Up from the STANISLAS Cohort.

Authors:  Abhinav Sharma; Amir Razaghizad; João Pedro Ferreira; Jean-Loup Machu; Erwan Bozec; Nicolas Girerd; Patrick Rossignol; Faiez Zannad
Journal:  Cardiology       Date:  2022-02-23       Impact factor: 2.342

7.  Sugar Is the Key Cause of Overweight/Obesity in Sugar-Sweetened Beverages (SSB).

Authors:  Lianlong Yu; Han Zhou; Fengjia Zheng; Jian Song; Yutong Lu; Xiao Yu; Changsheng Zhao
Journal:  Front Nutr       Date:  2022-06-28

Review 8.  Effect of Important Food Sources of Fructose-Containing Sugars on Inflammatory Biomarkers: A Systematic Review and Meta-Analysis of Controlled Feeding Trials.

Authors:  XinYe Qi; Laura Chiavaroli; Danielle Lee; Sabrina Ayoub-Charette; Tauseef A Khan; Fei Au-Yeung; Amna Ahmed; Annette Cheung; Qi Liu; Sonia Blanco Mejia; Vivian L Choo; Russell J de Souza; Thomas M S Wolever; Lawrence A Leiter; Cyril W C Kendall; David J A Jenkins; John L Sievenpiper
Journal:  Nutrients       Date:  2022-09-26       Impact factor: 6.706

9.  Effects of Sugar-Sweetened, Artificially Sweetened, and Unsweetened Beverages on Cardiometabolic Risk Factors, Body Composition, and Sweet Taste Preference: A Randomized Controlled Trial.

Authors:  Cara B Ebbeling; Henry A Feldman; Sarah K Steltz; Nicolle L Quinn; Lisa M Robinson; David S Ludwig
Journal:  J Am Heart Assoc       Date:  2020-07-22       Impact factor: 5.501

Review 10.  The role of sugar-sweetened beverages in the global epidemics of obesity and chronic diseases.

Authors:  Vasanti S Malik; Frank B Hu
Journal:  Nat Rev Endocrinol       Date:  2022-01-21       Impact factor: 47.564

  10 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.