Abstract
Background Observational studies have found associations between salt intake with obesity, body shape and composition; but the findings may be biased by residual confounding.
Objective To analyze the causal associations between salt intake and obesity measures in both sex-combined and sex-specific models.
Designs This was a two-sample Mendelian randomization study. Genome-wide association (GWA) studies of urinary sodium secretion (UNa, a surrogate for salt intake), body mass index (BMI), BMI-adjusted waist-to-hip ratio (WHR), body fat (BF) percentage and estimated glomerular filtration rate (eGFR) were identified. We initially extracted fifty SNPs associated with UNa at GWA significance level of 5 × 10−8, but further removed those SNPs with potential horizontal pleiotropy. Univariable MR and multivariable MR with adjustment for eGFR were performed. Inverse-variance weighted MR was performed as the primary analysis, with MR-Egger methods as sensitivity analysis. The potential bidirectional association between BMI and UNa was investigated. All exposure and outcomes were continuous, and the effect measure was linear regression coefficients (beta) and their 95% confidence intervals (95%CI).
Results UNa was causally associated with increased BMI in both men (eGFR-adjusted beta 0.443 (0.163 to 0.724)) and women (0.594 (0.333 to 0.855)). UNa caused BF percentage increase in men (0.622 (0.268 to 0.976)) and women (0.334 (0.007 to 0.662)). UNa significantly elevated BMI-adjusted WHR in men (0.321 (0.094 to 0.548)), but not in women (0.170 (−0.052 to 0.391)). Additionally, we found that BMI causally increased UNa (0.043 (0.023 to 0.063)).
Conclusions Salt intake increased BMI and BF percentage. Salt intake affects male body shape by increasing BMI-adjusted WHR, but showed no effects on female body shape. The bidirectional association between BMI and UNa suggested that salt reduction measures and weight reduction measures should be implemented simultaneously to break the vicious cycle and gain more health benefits.
Introduction
Obesity has been a global public health concern with increasingly substantial disease burdens (1–3). Obesity is believed to result from a complex relationship among psychosocial, biological and behavioral factors, such as genetic susceptibility, excessive calorie intake and physical inactivity (4). Other factors also contribute independently to development of obesity. (4,5) Obesity, commonly measured with body mass index (BMI, body weight in kilogram/squared height in meter), is a well-recognized risk factor for a range of health conditions and diseases, including cardiovascular diseases (6), diabetes (7,8) and cancers (9–11). BMI is easy to measure and calculate; nevertheless, body shape and body composition are sometimes stronger predictors for health outcomes than BMI (12,13).
Dietary salt intake is known to be associated with disease conditions, such as blood pressure (14), cardiovascular diseases (15), and diabetic complications (16). Salt intake can be measured with food frequency questionnaire or estimated from urinary sodium secretion (UNa), the latter of which is considered free of recall bias and more objective (17). Observational studies have suggested its effect on obesity, independent of total energy intake and fluid consumption. A systematic review published in 2015 (18) found that both dietary salt intake and UNa were positively associated with BMI, although with substantial heterogeneity. A recent cross-sectional study (19) showed that this positive association was directionally concordant across Asian and Western populations, with varying effect sizes. Another cross-sectional study (20) revealed an effect modified by sex, because the effect was observed only in women but not men. Salt intake was also demonstrated to be correlated with body shape and body composition measures, such as waist circumference (WC) (17,21), waist-to-hip ratio (WHR, waist circumference/hip circumference) (17), body fat (BF) mass and lean mass (20).
However, so far, the majority of the evidence has been generated from observational studies, mainly cross-sectional studies, in which residual confounding and reverse causation are likely. So far, few randomized trials on salt reduction have reported obesity-related outcomes (22,23). Mendelian randomization (MR), by using single nucleotide polymorphisms (SNPs) as instrumental variables to investigate the causal link between two phenotypes, is helpful in causal inference (24). The objective of this study is to examine the causal association between UNa with body mass, shape and composition using two-sample MR approach.
Materials and Methods
In this study, we used summary-level statistics from relevant GWA studies, where ethical approval and patient consent had been obtained. The exposure of interest was UNa, a commonly used surrogate marker for salt intake (25). The outcomes included BMI, body shape and body composition. For body shape, we used BMI-adjusted WHR as the primary outcome, while WC, hip circumference (HC), WHR, BMI-adjusted WC, BMI-adjusted HC were secondary outcomes. For body composition, we used BF percentage as the primary outcome, while whole-body lean mass (WLM) and appendicular lean mass (ALM) as secondary outcomes. Since sex has been suggested as a potential effect modifier in obesity-related outcomes (26–28), we investigated sex-specific effects in men and women separately for the outcomes, if possible.
SNP selection
A GWA study of 446,237 European-ancestry individuals from UK Biobank (29) identified 50 lead SNPs associated with UNa that reached significance level of 5 × 10−8 (Supplementary table 1). Spot urine sample collection and storage in UK Biobank have been described elsewhere (30). UNa concentration was measured by the ion-selective electrode method. The effect between genetic variants and log-transformed UNa concentration was adjusted for age, sex, array information and ancestral principal components. The SNPs are located in different gene regions and distributed independently (not in linkage disequilibrium), and explain 6.4% of total variance of the trait (29). If an SNP could not be matched in the GWA studies of an outcome, a proxy in linkage disequilibrium with the SNP (R2 > 0.80) would be identified; if no proper proxy was identified, the unmatched SNP would be removed from analysis. Relevant traits associated with the SNPs were also searched in the PhenoScanner v2 database at a significance level of 5 × 10−8 for potential pleiotropy.
Outcome data
We extracted relevant summary-level data of the outcomes from GWA studies of BMI (31), body shape (32) and body composition (33,34). The BMI GWA study (31) included 322,154 European-ancestry individuals (152,882 men and 171,963 women) from 125 cohorts. The effect was adjusted for age, squared age, ancestral principal components and study-specific covariates. The body shape GWA study (32) included 211,088 European-ancestry individuals (93,480 men and 116,742 women), and provided data for BMI-adjusted and unadjusted WC, HC, and WHR, additionally adjusted for age, squared age and ancestral principal components. The BF percentage GWA study (34) included 89,297 European-ancestry individuals (44,429 men and 45,525 women), and was adjusted for age, squared age and ancestral principal components. The GWA study for WLM and ALM (33) included 38,292 and 28,330 European-ancestry individuals, respectively, and was adjusted for age, sex, squared age, height, body fat percentage and ancestral principal components. Body composition (BF, WLM and ALM) was measured with either bioimpedance analysis or dual-energy X-ray absorptiometry (32,33). Sex-combined and sex-specific outcome data were available from the GWA studies of BMI, body shape and BF percentage. The GWA study for WLM and ALM provided sex-combined data only.
Covariate
Renal function is associated with both UNa secretion (35) and obesity (36,37), which is a potential confounder in the associations of interest. We identified a GWA study (38) of estimated glomerular filtration rate (eGFR), which included 567,460 European-ancestry individuals. Linear regression of log-transformed eGFR was adjusted for age, sex and ancestry principal components.
Statistical analysis
Previous studies (29) have suggested that some of the UNa-associated SNPs are related to adiposity phenotypes, which could introduce horizontal pleiotropy and violate one of the underlying assumptions of MR. To minimize potential pleiotropy, we further removed SNPs that: (1) were associated with the outcomes of interest at significance level of 5 × 10−8, or (2) were identified as outliers by the MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) outlier tests (39). Therefore, the included SNPs for each outcome may differ but remained a subset of the 50 SNPs identified by the UNa GWA study (29). Supplementary file 1 shows the characteristics of the SNPs included in the MR analysis for each outcome. F statistic was used to measure the strength of instrumental variables, and a F statistic > 10 suggests strong instrument variables. (40)
All GWA summary-level data were harmonized with the method of Hartwig et al. (41), so that the effect estimates of SNPs on exposure, outcome and covariate were shown for the same alleles. Univariable MR analyses were performed with inverse-variance weighted method as primary analysis, while weighted median method and MR-Egger regression method as sensitivity analysis. Multivariable MR analyses adjusted for eGFR with inverse-variance weighted method were performed as primary analysis, while MR-Egger method as secondary analysis. The intercept test for MR-Egger regression was used to examine potential residual horizontal pleiotropy, with a P value < 0.05 suggesting pleiotropy.
To investigate sex-specific effects, we extracted male-specific and female-specific outcome data from the GWA studies, where applicable (31,32,34), and performed MR analyses accordingly. Sex-specific effects on WLM and ALM were not estimated, because the data were unavailable. To examine the effect difference between men and women, we performed a heterogeneity test similar to that in meta-analysis, using Cochrane’s Q test. A P-value for Cochrane’s Q test < 0.10 indicates statistical significance.
To investigate potential bidirectional association between UNa and BMI, we performed MR analysis with BMI as the exposure and UNa as the outcome. Ninety-seven BMI-associated SNPs were extracted from the BMI GWA study (31), but 84 SNPs were finally included after removing those unmatched with UNa GWA data (29) and/or identified as outliers by the MR-PRESSO outlier tests. Similarly, univariable and multivariable MR analyses were performed.
A second way to deal with potential confounding effect of eGFR, alternative to the multivariable MR, was to remove the SNPs associated with eGFR. We found that rs1260326 (located in GCKR gene) was associated with eGFR at a GWA significance level of 5 × 10−8 in Phenoscanner v2 database; therefore, a sensitivity analysis by removing rs1260326 in univariable MR was performed as sensitivity analysis. Additional sensitivity analyses were performed by including outlier SNPs detected by the MR-PRESSO outlier tests. Since all outcomes were continuous, the effect measure was linear regression coefficient beta with its 95% confidence interval (CI), which should be interpreted as the change in units of standard deviation (SD) in the outcome when the log-UNa level increases 1 SD. Statistical significance was indicated by p value < 0.05. We did not apply Bonferroni correction, although we employed three primary outcomes, because their effects were estimated from three independent samples of individuals, instead of one. All statistical analysis was conducted with “MendelianRandomization” (version 0.4.1) and “MRPRESSO” (version 1.0) packages in R environment. Supplementary figure 1 shows the flowchart of data selection and analysis.
We planned to perform power calculations based on the method of Brion et al. (42), which required the beta estimates from our study and another beta estimates from observational analysis. However, the UNa GWA study analyzed log-UNa concentration (in log-mmol/L), while previous observational studies analyzed UNa concentration in a natural unit (mmol/L) or converted it to 24-hour UNa (mmol/d or g/d) or daily salt intake (g/d), thus making it impossible to make direct comparisons of the two beta estimates in the absence of individual participant data.
Results
The proportions of exposure variance explained by the instrument variables ranged from 4.036% to 5.565% (mean 4.780%), and the according F statistics ranged from 536.146 to 613.340 (mean 587.078), suggesting strong instrumental variables. (Table 1) Supplementary file 1 showed detailed information on the SNPs included in the MR analysis for each outcome.
Univariable MR
In sex-combined analyses, MR analyses with inverse-variance weighted method showed positive causal associations between UNa with BMI (beta (95%CI): 0.392 (0.149 to 0.634); Figure 1(A)), BMI-adjusted WHR (0.290 (0.101 to 0.478); Figure 2(A)) and BF percentage (0.475 (0.208 to 0.741); Figure 3(A)). UNa was not associated with BMI-adjusted WC, BMI-adjusted HC, WLM or ALM. MR-Egger intercept tests revealed no horizontal pleiotropy. (Table 1)
BMI was positively related to UNa with the effect size of 0.042 (0.022 to 0.062), consistent with the weighted median estimator (0.042 (0.019 to 0.066)). MR-Egger intercept test showed no pleiotropy (p = 0.119). This suggested a bidirectional association between UNa and BMI.
In men, UNa was positively associated with BMI (0.417 (0.113 to 0.720); Figure 1(B)), BMI-adjusted WHR (0.339 (0.114 to 0.564); Figure 2(B)) and BF percentage (0.617 (0.270 to 0.963); Figure 3(B)), but negatively with BMI-adjusted HC (−0.396 (−0.699 to −0.092)). No pleiotropy was observed according to MR-Egger intercept test, except WHR (p = 0.029), which was a secondary outcome. UNa was associated with WC but not BMI-adjusted WC, which suggested that BMI might mediate the association between UNa with WC. (Table 2)
In women, UNa increased BMI (0.571 (0.295 to 0.848); Figure 1(C)), BMI-adjusted WHR (0.249 (0.027 to 0.472); Figure 2(C)) and BF percentage (0.348 (0.027 to 0.669); Figure 3(C)). UNa was associated with HC and WC, but the associations became insignificant after adjustment for BMI. No pleiotropy was detected. (Table 1)
Multivariable MR
The results of multivariable MR analyses were largely consistent with the univariable analyses. Overall, UNa was associated with increased BMI and BF percentage in both men and women. Higher UNa affected male body shape by increasing BMI-adjusted WHR and decreasing BMI-adjusted HC. UNa seemed to have little effects on female body shape after adjusted for BMI and eGFR. Formal heterogeneity test showed significant difference between men and women on BMI-adjusted HC, but not other outcomes. (Table 2) The eGFR-adjusted effect of BMI on UNa was 0.043 (0.023 to 0.063) with no pleiotropy (p = 0.149), indicating a bidirectional association between BMI and UNa.
In sensitivity analysis including the SNPs that were identified as outliers by MR-PRESSO outlier tests, both the univariable and multivariable MR showed results consistent with primary results. (Supplement table 2 and supplementary table 3) Removing the eGFR-associated SNP rs1260326 also generated concordant results. (Supplementary table 4)
The scatter plots showed a potential outlier SNP rs4803378, located in gene CTC-490E21.10 and associated with urine creatinine according to PhenoScanner v2 database, the post hoc sensitivity analyses by removing it showed similar results to the primary results. (Supplementary table 5)
Discussion
In this MR study, salt intake, measured with UNa, was causally associated with increased BMI and BF percentage. Higher salt intake changed body shape towards central obesity in men, by reducing BMI-adjusted HC and increasing BMI-adjusted WHR, but probably not in women. Additionally, salt intake and BMI were bidirectionally positively related.
These results were consistent with previous findings. Moosavian et al.(18) found in a meta-analysis that both dietary salt intake and UNa were positively associated with BMI. This study computed the effect measure as the mean difference of BMI between highest and lowest categories of exposure, which introduced heterogeneity due to various cutoff values used and made it difficult to directly compare the effect size with other studies. Zhou et al.(19) found that for each 1gram per day (g/d) increase in salt intake, BMI increased by 0.42 and 0.52 kg/m2, overweight/obesity risk increased by 29% and 24% in UK and US population, respectively, after adjustment for total energy intake and physical activity..
Salt intake showed an effect on body composition by increasing BF percentage, which is consistent with previous findings. An increase of 1g/d salt intake was associated with an increase of 0.91kg BF mass in UK population (21), and of 0.79kg BF mass and 0.44% BF percentage in US population (20). Oh et al. (43) observed this positive association but only in people under 65 years old. Zhu et al. (44) found salt intake was associated with increased BF mass, BF percentage, and subcutaneous abdominal adipose tissue. Interestingly, a randomized trial (45) found that low-salt diet significantly decreased body weight and extracellular water, but reductions in BF mass and lean mass were insignificant. However, this trial was limited by small sample size (84 obese people) and short follow-up period (2 months), which might make it underpowered to detect significant reductions. Salt intake has an instant effect on water retention, but its effect on fat accumulation may take longer.
We found null effect of salt intake on lean mass, neither WLM nor ALM. This remains uncertain. Previously, Ma et al. (21) showed that 1 g/d increase in salt intake was associated with significant WLM increase by 0.32kg, but 1 g/2000kcal increase in salt intake density insignificantly reduced WLM by 0.007kg. Zhang et al. (20) found that salt intake significantly increased WLM by 0.53kg, but the effect became insignificant when using salt intake density (g/kcal) to adjust for energy intake. It seems that there are two possible effects of salt intake on body composition: (1) increase both fat and lean mass, but fat mass increases at a higher rate; (2) increase fat mass solely. The first explanation was supported by the findings that salt intake was associated with BF mass increase more than lean mass in both UK (0.91kg versus 0.32kg for fat and lean mass increase, respectively) (21) and US populations (0.79kg versus 0.53kg) (20). The second explanation was supported by the the null effect of salt intake on lean mass after adjustment for energy intake (by using salt intake density). (20,21) Our finding may also support the second explanation. Nevertheless, more research is warranted on this topic.
There are several potential explanations for how salt intake affects body weight and composition. It has been well recognized that excessive salt intake increases thirst and fluid consumption, thus increasing extracellular water while keeping urine volume almost unchanged. (46) Higher dietary salt intake also increases sweetened-beverage intake (47), food appetite and calorie intake (18,48), which all contribute to weight gain. In addition, high sodium intake has independent biological effects on fat accumulation. Lanaspa et al. (49) demonstrated that high salt intake activated the aldose reductase-fructokinase pathway in liver and hypothalamus, causing endogenous fructose production and the development of leptin resistance and hyperphagia, which finally cause obesity. Fonseca-Alaniz et al. (50) also observed that the high salt-induced adiposity in rats was characterized by high plasma leptin concentrations and adipocyte hypertrophy, which might be mediated by lipogenic capacity of white adipose mass. Lee et al. (51) revealed that high salt increased the expression of adipogenic and lipogenic genes (such as PPAR-γ, SREBP1c, ACC and C/EBPα), but decreased lipolysis gene expression (such as AMPK). However, its effect on lean mass requires more replications and mechanism research.
We found that salt intake had an impact on body shape. Body shapes more predictive of obesity-related mortality and health outcomes (52,53). Overall, salt intake increased both BMI-adjusted and unadjusted WHR, leading to central obesity, which is in line with previous findings (17,18,20,44). However, we observed an effect modification caused by sex. Salt intake increased WHR in both sexes, but a significant increase in BMI-adjusted WHR was seen in men only. Nam et al. (17) also found a similar sex difference. This may suggest that salt intake affects body shape independent of BMI in men, but its effect on female body shape is, at least partly, mediated by BMI.
This sexual dimorphism of body fat accumulation and distribution may be related to sex hormones (26,54), sex chromosome (55) and sex-specific autosomal genetic heritability (27,28,56,57). Estrogens favor fat accumulation in subcutaneous depot in women. Estrogens and their receptor regulation also enhance the expandability of adipose cells in subcutaneous depot and inhibits it in visceral depots. (26) Elevated free androgens have also been shown to be related to increased abdominal visceral fat accumulation and increased WC. (58) The Four Core Genotypes mouse model, has shown that X chromosome dosage is a risk factor for fat accumulation, but the presence of Y chromosome is not (56). Other studies have also shown that that adipocytes from different subcutaneous depots (abdominal versus gluteal) are developmentally distinct in a sex-dependent manner (59,60). The GWA study of BMI-adjusted WHR showed a higher heritability in women than in men (2.4% versus 0.8%) (32). This indicates that male body shape is more likely to be influenced by environmental factors (56), such as salt intake. However, mechanistic studies are needed to elucidate this issue.
In this MR study, we observed that high BMI was causally linked to higher UNa. High UNa could be resulted from two possible mechanisms: obesity-related glomerulopathy and high salt intake. Obesity is known to increase risk of glomerulopathy, characterized with altered renal hemodynamics, increased eGFR and tubular sodium reabsorption (61); however, the overall effect of glomerulopathy on UNa remains unclear. On the other hand, in our study, the association between BMI and UNa remained significant after adjustment for renal function, therefore, we believe that salt intake has independent mediation effect on the association between BMI and UNa, that is, BMI increases UNa by increasing salt intake, independent of renal function. Li et al. (62) found that overweight/obese people had lower salt sensitivity and higher salt preference, and consumed 2.3g/d more dietary salt than normal-weight people
The bidirectional association between obesity and salt intake indicates the formation of a vicious cycle, in which people with higher BMI are prone to consume more salt, which then further increases their BMI. Considering that high salt intake and obesity, both prevalent worldwide, are shared risk factors for many chronic diseases and that they also have their own adverse health consequences (63), urgent interventions are required to break this cycle in high-risk individuals.
Salt reduction and weight reduction programs have been widely implemented in various settings, but separately. Salt reduction programs do not consider weight change as an outcome, while weight reduction programs pay little emphasis on salt reduction measures. We are not advocating weight control by salt reduction alone. Instead, our findings provide evidence to combine salt reduction and weight reduction measures, to facilitate people to acquire a healthier lifestyle and achieve higher health benefits. Although the actual effect size of salt reduction on weight reduction has not been precisely estimated yet, which warrants future well-designed randomized trials, salt reduction per se costs less and brings additional health benefits. (22,43,64)
There are some limitations. First, we did not include all identified UNa-associated SNPs because unmatched data and pleiotropy issues; however, our results are consistent in primary analyses and sensitivity analyses, and with previous studies. Second, more accurate measurement of UNa is via 24-hour urine sample collection, but the UK Biobank uses spot urine sample. The 24-hour urine collection method is more expensive, time-consuming, labor-intensive and difficult for an epidemiological study as large as UK Biobank with half a million participants. Previous studies showed that spot UNa is an acceptable surrogate for 24-hour UNa and dietary salt intake, and has been commonly employed elsewhere (25,30). Third, there might still be residual pleiotropy, although we have removed SNPs that were significantly associated with outcomes and that were identified as outliers by the MR-PRESSO outlier test. MR-Egger intercept test showed pleiotropy for a few secondary outcomes (e.g., HC, WC) only. Fourth, the effect size estimated in this study was not directly comparable to other studies in which the effect was measured at natural units of outcome by 1 g/d salt intake increase. Although several equations have been proposed to estimate salt intake in g/d from 24h urine sodium secretion in g/d or spot urine sodium concentration in mmol/L (25), the UNa GWA study log-transformed the original spot UNa concentration, which makes back-transformation of beta estimate to natural unit impossible. Due to the same reason, formal power calculations were not performed as we planned, which was another limitation, but our study has employed strong instrumental variables and so far the largest GWA studies available. Lastly, we only included European-ancestry population, thus generalizing the results to other ethnicities requires cautions.
Conclusion
Salt intake increases BMI and BF percentage, and changes male body shape towards central obesity by increasing BMI-adjusted WHR. Given the bidirectional association between salt intake and obesity, salt reduction measures and weight reduction measures should be implemented together to achieve more health benefits.
Data Availability
Data are available upon request.
Conflict of interest
None.
Ethical approval
Not required.
Data sharing
Data described in the manuscript, code book and analytic code will be made available upon request.
- Abbreviations
- ALM
- appendicular lean mass
- BF
- body fat
- BMI
- body mass index
- eGFR
- estimated glomerular filtration rate
- GWA
- genome-wide association
- HC
- hip circumference
- MR
- Mendelian randomization
- PRESSO test
- Pleiotropy RESidual Sum and Outlier test
- SNP
- single nucleotide polymorphism
- UNa
- urinary sodium secretion
- WC
- waist circumference
- WHR
- waist-to-hip ratio
- WLM
- whole body lean mass