Hostname: page-component-586b7cd67f-t7fkt Total loading time: 0 Render date: 2024-11-25T16:03:19.597Z Has data issue: false hasContentIssue false

Low hydration status may be associated with insulin resistance and fat distribution: analysis of the Korea National Health and Nutrition Examination Survey (KNHANES) 2008–2010

Published online by Cambridge University Press:  19 March 2020

Hyang K. Min
Affiliation:
Department of Internal Medicine, Nowon Eulji Medical Center, Eulji University, Seoul, Republic of Korea
Hyun Y. Ko
Affiliation:
Department of Internal Medicine, Nowon Eulji Medical Center, Eulji University, Seoul, Republic of Korea
Jin T. Kim
Affiliation:
Department of Internal Medicine, Nowon Eulji Medical Center, Eulji University, Seoul, Republic of Korea
Lise Bankir
Affiliation:
Sorbonne Université, INSERM Unité 1138, Centre de Recherche des Cordeliers, Paris, France
Sung W. Lee*
Affiliation:
Department of Internal Medicine, Nowon Eulji Medical Center, Eulji University, Seoul, Republic of Korea
*
*Corresponding author: Sung W. Lee, fax +82-2-971-8212, email [email protected]
Rights & Permissions [Opens in a new window]

Abstract

We aimed to identify the association of hydration status with insulin resistance (IR) and body fat distribution. A total of 14 344 adults participated in the Korea National Health and Nutrition Examination Survey 2008–2010. We used urine specific gravity (USG) to indicate hydration status, and HOMA-IR (homoeostasis model assessment of IR) and trunk:leg fat ratio (TLR) as primary outcomes. In multivariate logistic regression, the OR per 0·01 increase in USG for high IR was 1·303 (95 % CI 1·185, 1·433; P < 0·001). In multivariate generalised additive model plots, increased USG showed a J-shaped association with logarithmic HOMA-IR, with the lowest Akaike’s information criterion score of USG 1·030. Moreover, increased USG was independently associated with increased trunk fat, decreased leg fat and increased TLR. In mediation analysis, the proportion of mediation effects of USG on TLR via IR was 0·193 (95 % CI 0·132, 0·285; P < 0·001), while the proportion of mediation effects of USG on IR via TLR was 0·130 (95 % CI 0·086, 0·188; P < 0·001). Increased USG, a sign of low hydration status and presumably high vasopressin, was associated with IR and poor fat distribution. Direct effect of low hydration status may be more dominant than indirect effect via IR or fat distribution. Further studies are necessary to confirm our findings.

Type
Full Papers
Copyright
© The Authors 2020

Vasopressin (AVP), also known as antidiuretic hormone, is well known for its role in conserving water via V2 receptors expressed in the renal collecting ducts(Reference Bankir1). However, there are several unfamiliar actions of AVP beyond its antidiuretic effect. First, via its original antidiuretic action, increased AVP was associated with glomerular hyperfiltration(Reference Bouby, Ahloulay and Nsegbe2,Reference Clark, Sontrop and Huang3) . We recently confirmed that ‘sub-morbid dehydration’, a condition in which AVP is presumably elevated, was associated with glomerular hyperfiltration in general Korean population(Reference Min, Sung and Lee4).

Second, AVP influences glucose metabolism in different ways. AVP stimulates glycogenolysis and gluconeogenesis via V1a receptors expressed in the liver(Reference Koshimizu, Nakamura and Egashira5,Reference Muscogiuri, Barrea and Annunziata6) . AVP stimulates the release of either glucagon or insulin, depending on extracellular glucose concentration, through the activation of V1b receptors in pancreatic islets(Reference Abu-Basha, Yibchok-Anun and Hsu7). AVP has also been shown to stimulate cortisol release through activation of the hypothalamus–pituitary–adrenal (H-P-A) axis, via V1b receptor expressed in the anterior pituitary, thus further influencing glucoregulatory mechanisms(Reference Koshimizu, Nakamura and Egashira5,Reference Muscogiuri, Barrea and Annunziata6) . This unfamiliar role of AVP in glucose metabolism has been recently supported in several clinical studies showing that increased copeptin, the C-terminal part of the AVP precursor peptide(Reference Christ-Crain and Fenske8), used as a surrogate marker of AVP, was associated with increased risk of insulin resistance (IR)(Reference Wannamethee, Welsh and Papacosta9Reference Jensen, Bjornstad and Johnson11), the metabolic syndrome(Reference Saleem, Khaleghi and Morgenthaler12Reference Vintilă, Gheorghiu and Caragheorgheopol14) and diabetes(Reference Muscogiuri, Barrea and Annunziata6,Reference Wannamethee, Welsh and Papacosta9,Reference Enhorning, Wang and Nilsson15Reference Roussel, El Boustany and Bouby18) .

Third, AVP may also affect fat metabolism. V1a receptor is expressed in both brown and white adipose tissues, while V1b receptor is expressed only in white adipose tissue(Reference Hiroyama, Aoyagi and Fujiwara19). It is suggested that activation of V1a receptor is associated with enhanced lipogenesis(Reference Hiroyama, Aoyagi and Fujiwara19), while activation of V1b receptor is associated with enhanced lipolysis(Reference Hiroyama, Fujiwara and Nakamura20) along with enhanced IR in adipose tissue(Reference Fujiwara, Hiroyama and Sanbe21). Unlike glucose metabolism, however, the potential effect of AVP on fat metabolism has rarely been evaluated in clinical studies(Reference Enhorning, Bankir and Bouby17).

Physiologically, a low hydration status is a major stimulus for the increase in AVP level(Reference Bankir1). Although hydration status can be evaluated in several ways(Reference Garrett, Rae and Fletcher22), urine specific gravity (USG) is frequently used(Reference Kavouras23), particularly in large-scale populations(Reference Baron, Courbebaisse and Lepicard24). Increased USG may indicate a low hydration condition and thus, presumably, elevated AVP(Reference Perrier, Vergne and Klein25). Therefore, we hypothesised that a low hydration status or presumably elevated AVP status may be associated with increased risk of IR and poor fat distribution. To investigate our study hypothesis, we analysed the relationship of USG with IR and fat distribution in community-dwelling Korean adults, using data from the Korea National Health and Nutrition Examination Survey (KNHANES).

Methods

Participants

The KNHANES has been performed periodically since 1998 to assess the health and nutritional status of the civilian, non-institutionalised Korean population. Participants are being selected using proportional allocation and systematic sampling with multistage stratification. This study used data from the KNHANES 2008–2010. Of 36 188 candidates, 29 235 agreed to participate in the survey (participation rate, 80·8 %). The protocol comprised a health questionnaire survey, health examination and nutrition survey. Among the 29 235 participants, 7664 people whose age was <20 years were excluded. Of 21 571 adult participants, 5210 people with missing body fat measurements, eighty-seven with missing anthropometric measurements, 668 with missing fasting insulin and glucose levels, 885 with missing USG values and 377 with fasting time <8 h were further excluded. Therefore, this study included 14 344 participants.

Ethics statement

The study protocol complied with the Declaration of Helsinki, and full approval of the study was obtained from the Institutional Review Board of the Korea Centers for Disease Control (Institutional Review Board numbers: 2008-04EXP-01-C, 2009-01CON-03-2C and 2010-02CON-21-C). All data were fully anonymised before we accessed them. KNHANES participants provided informed written consent to have their data used in research.

Urine specific gravity, insulin resistance and body fat measurements

All the blood and urine samples were collected in the morning after overnight fast. Samples were transported and analysed at the central laboratory (Neodin Medical Institute). USG and proteinuria were measured with a dipstick method using the UriSys 2400 analyser (Roche), and analysing range was 1·000 to 1·050 for USG and – to 4+ for proteinuria.

The homoeostasis model assessment (HOMA) was used to calculate IR using the following equation: HOMA-IR = fasting insulin (μIU/ml) × fasting glucose (FG) (mmol/l)/22·5(Reference Matthews, Hosker and Rudenski26). With the participant in the erect position, height was measured with a stadiometer (Seca 225; Seca) to the nearest 0·1 cm. Body weight was measured with the participant in a light gown with bare feet, with the use of digital scales (GL-6000-20; G-Tech) to the nearest 0·1 kg, and BMI was calculated by dividing the weight by the square of the height (kg/m2). Waist circumference (WC) was measured with a tape measure (Seca 200; Seca) to the nearest 0·1 cm at the midpoint between the lower border of the rib cage and the highest point of the iliac crest of the participant.

Regional fats were measured using dual-energy X-ray absorptiometry with the QDR 4500A analyzer (Hologic Inc.) and calculated using Hologic Discovery software (version 3.1). The subregions of the head, arms, legs and trunk were delineated during the whole-body scan for assessment of regional fat amount. The trunk was defined using a horizontal line below the chin and vertical lines passing through the femoral neck comprising the right and left ribs, thoracic and lumbar spines and pelvic area. The arm regions were separated from the trunk at the levels of the shoulder joint; arm fat was defined as the sum of fat mass in both arms, and leg fat as the sum of fat mass in both legs. Trunk:leg fat ratio (TLR) was used as a fat distribution marker(Reference Ito, Nakasuga and Ohshima27).

Other measurements

A standardised interview was conducted in the homes of the participants to collect information regarding demographic variables, medical history, medications used and other health-related variables. Blood pressure (BP) was measured three times in accordance with the standard protocol, and mean values of the three measurements were used as the representative BP. Total cholesterol, HDL-cholesterol, TAG, FG, and urinary Na levels were measured using the Hitachi 7600 Automatic Analyzer (Hitachi), and fasting insulin was measured using the 1470 WIZARD gamma counter analyzer (PerkinElmer). The estimated glomerular filtration rate was calculated using the Chronic Kidney Disease Epidemiology Collaboration equation(Reference Levey, Stevens and Schmid28). Plain water and energy intake were estimated using the single day 24 h recall method.

Definitions

The following five metabolic disorders were defined based on the recommendations of the International Diabetes Federation(Reference Alberti, Zimmet and Shaw29). Central obesity was defined as WC ≥ 90 cm for men and WC ≥ 80 cm for women. Raised BP was defined as systolic BP ≥ 130 mmHg, diastolic BP ≥ 85 mmHg, or treatment with antihypertensive drugs. Raised FG was defined as FG ≥ 5·6 mmol/l, treatment with insulin or oral anti-diabetic drugs, or a previous physician’s diagnosis of type 2 diabetes. Raised TAG was defined as TAG ≥ 1·7 mmol/l. Reduced HDL-cholesterol level was defined as HDL-cholesterol < 1·03 mmol/l for men and HDL-cholesterol < 1·29 mmol/l for women, or having a specific treatment for dyslipidaemia.

Previous CVD was defined as a physician’s diagnosis of cerebrovascular disease, angina or myocardial infarction. Alcohol drinking was defined as drinking alcoholic beverage more than twice a week. High-intensity exercise was defined as doing strenuous exercise (e.g. soccer, climbing) more than 3 d/week. Mid-intensity exercise was defined as doing moderate exercise (e.g. slow swimming, table tennis) more than 3 d/week. Proteinuria was defined as a protein level of 1+ or higher in the dipstick urinalysis. Furthermore, high IR was defined as a score in the highest quartile of the HOMA-IR (≥2·8).

Statistical analysis

For describing the baseline characteristics of the participants, we used sex-specific quartiles of USG because plasma vasopressin and urine concentration are known to be higher in men than women(Reference Perucca, Bouby and Valeix30). The first quartile of USG was <1·014 in women and <1·016 in men, the second quartile was 1·014–1·017 in women and 1·016–1·019 in men, the third quartile was 1·018–1·021 in women and 1·020–1·023 in men and the fourth quartile was ≥ 1·022 in women and ≥ 1·024 in men. The distributions of continuous variables were evaluated using histograms and Q–Q plots. The TAG, FG, fasting insulin levels and HOMA-IR score were not normally distributed. Normally distributed continuous variables were expressed as mean values and standard deviations, non-normally distributed continuous variables as medians (interquartile ranges (IQR)), and categorical variables as percentages. P trend was analysed for normally distributed continuous variables using a linear term of one-way ANOVA, for non-normally distributed continuous variables using the Jonckheere–Terpstra test, and for categorical variables using a linear-by-linear association. Differences were analysed using the Bonferroni post hoc analysis of one-way ANOVA for normally distributed continuous variables, the Mann–Whitney U test for non-normally distributed continuous variables, and the χ 2 test for categorical variables. OR and 95 % CI of USG for high IR were analysed using logistic regression analysis. The relationship between USG and regional fat depots was analysed using linear regression analysis. On multivariate analysis, covariates were chosen based on clinical and statistical relevance. For the potential non-linear association, a multivariate generalised additive model for Gaussian distributions was adapted to visualise the associations between USG and logarithmic HOMA-IR using ‘mgcv’ package in R Statistics (version 3.03). We used Akaike’s information criterion (AIC) as the primary measure of the model fit(Reference Hiroyama, Aoyagi and Fujiwara19). In the AIC analysis, lower scores within the data set indicated a better model fit, and the USG with the lowest AIC value for the association between USG and logarithmic HOMA-IR was assumed to be the threshold value. Mediation analysis was performed using ‘mediation’ package in R Statistics (version 3.03). Mediator and outcome were interchanged between logarithmic HOMA-IR and TLR, while treatment (USG) and pre-treatment covariates (age, sex, smoking and drinking status, suicide thought, previous CVD, proteinuria, estimated glomerular filtration rate, systolic and diastolic BP, WC, BMI, HDL-cholesterol, logarithmic TAG, lean body mass, arm and head fat, high- and mid-intensity exercise, plain water and energy intake, urine Na, and leucocytes) were fixed. We performed 2000 times simulation for the results. Average of all causal mediation effects and direct effects was presented. A P value of <0·05 was considered statistically significant. All analyses, unless otherwise specified, were performed using SPSS version 22 (release 2013; IBM Corp.).

Results

Among the 14 344 participants, 45·1 % were men and the mean values of age were 49·8 (sd 15·8) years. The mean values of BMI and WC were 23·7 (sd 3·3) kg/m2 and 81·4 (sd 9·8) cm, respectively. The prevalence rates of central obesity, raised TAG, reduced HDL-cholesterol, raised BP and raised FG were 37·7, 29·3, 47·6, 41·4 and 29·1 %, respectively. The mean values of body fat mass were 9·15 (sd 3·49) kg in the trunk region, 5·32 (sd 1·90) kg in the leg region, 1·99 (sd 0·77) kg in the arm region and 0·94 (sd 0·14) kg in the head region. The mean values of HDL-cholesterol were 1·24 (sd 0·28) mmol/l, and the medians of TAG, glucose and insulin levels were 1·22 (IQR 0·82, 1·84) mmol/l, 5·16 (IQR 4·83, 5·61) mmol/l and 61·4 (IQR 48·6, 79·2) pmol/l, respectively.

Mean values of USG were 1·017 (sd 0·006) in women and 1·019 (sd 0·005) in men (P < 0·001). The baseline characteristics of the participants according to sex-specific USG quartiles are depicted in Table 1. The age of the participants decreased, and the proportion of men increased as the sex-specific USG quartile increased. With the increase in sex-specific USG quartile, the rate of current smoking increased. People in the higher sex-specific USG quartile tended to drink less plain water, took more energy content and performed high-intensity exercise more. As the sex-specific USG quartile increased, the rates of previous CVD, raised BP, central obesity, reduced HDL-cholesterol, raised TAG and raised FG decreased. Moreover, as the sex-specific USG quartile increased, the estimated glomerular filtration rate increased.

Table 1. Clinical characteristics of the study population according to sex-specific urine specific gravity (USG) quartiles§

(Mean values and standard deviations for continuous variables; numbers and percentages for categorical variables; medians and interquartile ranges)

BP, blood pressure; FG, fasting glucose.

P trend was analysed for normally distributed continuous variables by a linear term of one-way ANOVA, for non-normally distributed continuous variables by Jonckheere–Terpstra tests, and for categorical variables by a linear-by-linear association. *, † and ‡ meant P < 0·008 when compared with Q1, Q2 and Q3, respectively, using Bonferroni post hoc analysis of one-way ANOVA for normally distributed continuous variables, Mann–Whitney U tests for non-normally distributed continuous variables and χ 2 tests for categorical variables.

§ Conversion factors from conventional units to SI units were ×6·945 for fasting insulin, ×0·0259 for HDL-cholesterol, ×0·0113 for TAG, ×0·0555 for FG.

To convert kcal to kJ, multiply by 4·184.

We analysed the relationship between USG and IR (Table 2). On univariate logistic regression analysis, per 0·01 increase in USG was associated with 1·130 times higher odds for high IR. This relationship between USG and high IR was independent from various confounders: the OR per 0·01 increase in USG for high IR was 1·303 (95 % CI 1·185, 1·432; P < 0·001). Furthermore, we identified a J-shaped association between USG and logarithmic HOMA-IR with the lowest AIC score of USG 1·030 (Fig. 1). We also explored the relationship between USG and four fat depots using linear regression analysis (Table 3). On multivariate analysis, increased USG was not associated with fat mass in the arm and head regions, but it was independently associated with increased trunk fat and decreased leg fat, resulting in an increased TLR.

Table 2. Association between urine specific gravity (USG) and high insulin resistance*

(Odds ratios and 95 % confidence intervals)

* The first quartile (Q1) of USG was <1·014 in women and <1·016 in men, the second quartile (Q2) was 1·014–1·017 in women and 1·016–1·019 in men, the third quartile (Q3) was 1·018–1·021 in women and 1·020–1·023 in men, and the fourth quartile (Q4) was ≥ 1·022 in women and ≥ 1·024 in men. High insulin resistance was defined as highest quartile of homeostasis model assessment of insulin resistance (≥2·8). OR and 95 % CI were calculated using logistic regression analysis. In model 1, cardio-metabolic confounders (systolic and diastolic blood pressure, waist circumference, BMI, HDL-cholesterol and TAG) were entered as covariates. In model 2, body composition indices (lean body mass, trunk fat, leg fat, arm fat and head fat) along with variables in model 1 were entered as covariates. In model 3, other variables (age, sex, smoking and drinking status, suicide thought, previous CVD, proteinuria, estimated glomerular filtration rate, high- and mid-intensity exercise, plain water and energy intake, urine Na, and leucocytes) were entered as covariates.

Fig. 1. Relationship between urine specific gravity and insulin resistance. The dashed line indicates 95 % CI for value of the smoothed logarithmic homoeostasis model assessment of insulin resistance (HOMA-IR) using multivariate generalised additive model analysis after adjusting for age, sex, smoking and drinking status, suicide thought, previous CVD, proteinuria, estimated glomerular filtration rate, systolic and diastolic blood pressure, waist circumference, BMI, HDL-cholesterol, TAG, lean body mass, trunk fat, leg fat, arm fat, head fat, high- and mid-intensity exercise, plain water and energy intake, and urine Na as covariates. AIC, Akaike’s information criterion.

Table 3. Relationship between urine specific gravity and four fat depots*

(β Values and 95 % confidence intervals)

* β Values and 95 % CI were analysed using linear regression analysis. In multivariate analysis, the covariates were logarithmic homeostasis model assessment of insulin resistance, systolic and diastolic blood pressure, waist circumference, BMI, logarithmic TAG, HDL-cholesterol, lean body mass, trunk fat, leg fat, arm fat, head fat, age, sex, smoking and drinking status, suicide thought, previous CVD, proteinuria, estimated glomerular filtration rate, high- and mid-intensity exercise, plain water and energy intake, urine Na, and leucocytes. When fat depots were chosen as an outcome, they were excluded from the model.

We performed mediation analyses to identify the mutual relationship among USG, IR and TLR (Figs. 2 and 3). The average direct effect per 0·01 increase in USG on TLR was 0·040 (95 % CI 0·026, 0·053; P < 0·001), while average causal mediation effect per 0·01 increase in USG via IR was 0·010 (95 % CI 0·007, 0·012; P < 0·001), resulting in the proportion of mediation effects of 0·193 (95 % CI 0·132, 0·285; P < 0·001) as shown in Fig. 2. The average direct effect per 0·01 increase in USG on logarithmic HOMA-IR was 0·052 (95 % CI 0·037, 0·067; P < 0·001), while average causal mediation effect per 0·01 increase in USG via TLR was 0·008 (95 % CI 0·005, 0·010; P < 0·001), resulting in the proportion of mediation effects of 0·130 (95 % CI 0·086, 0·188; P < 0·001) as shown in Fig. 3.

Fig. 2. Mediation analysis of insulin resistance (IR) for the relationship between urine specific gravity (USG) and fat distribution. IR was estimated by homoeostasis model assessment (HOMA). Outcome was trunk:leg fat ratio (TLR), treatment was USG, and mediator was logarithmic HOMA-IR. Overall effects of treatment on outcome were adjusted by mediator and pre-treatment covariates using multivariate linear regression analysis. Mediator was modelled by treatment and pre-treatment covariates using multivariate linear regression analysis.

Fig. 3. Mediation analysis of fat distribution for the relationship between urine specific gravity (USG) and insulin resistance. Outcome was logarithmic homoeostasis model assessment of insulin resistance (HOMA-IR), treatment was USG, and mediator was trunk:leg fat ratio (TLR). Overall effects of treatment on outcome were adjusted by mediator and pre-treatment covariates using multivariate linear regression analysis. Mediator was modelled by treatment and pre-treatment covariates using multivariate linear regression analysis.

Discussion

The most well-known function of AVP is antidiuresis via V2 receptor(Reference Bankir1). However, according to recent studies, AVP may have additional pleiotropic effects including AVP-induced glomerular hyperfiltration(Reference Bouby, Ahloulay and Nsegbe2,Reference Clark, Sontrop and Huang3) and AVP-induced hyperglycaemia(Reference Koshimizu, Nakamura and Egashira5,Reference Muscogiuri, Barrea and Annunziata6) . Theoretically, AVP may also suppress insulin sensitivity via V1b receptor(Reference Fujiwara, Hiroyama and Sanbe21), while increased AVP has antilipolytic effects via V1a receptor(Reference Hiroyama, Aoyagi and Fujiwara19), suggesting that increased AVP may participate in IR and obesity. AVP is difficult to measure(Reference Christ-Crain and Fenske8). The measurement of copeptin, used as a surrogate marker of AVP(Reference Christ-Crain and Fenske8), has not been incorporated in the ordinary practice. In this context, measuring the hydration status, based on USG, an index of urine concentration, may be a satisfactory phenotypic marker of AVP(Reference Perrier, Vergne and Klein25). Among several hydration markers(Reference Garrett, Rae and Fletcher22), USG is appropriate for large-scale populations(Reference Kavouras23,Reference Baron, Courbebaisse and Lepicard24) . Therefore, we hypothesised that high USG may be associated with increased risk of IR and poor fat distribution. In our large cohort, we identified that increased USG is independently associated with increased IR and poor fat distribution. We also found that the effect of USG on poor fat distribution is partially mediated by increased IR, while the hazard of USG on IR is partly mediated by poor fat distribution.

USG may be influenced by various conditions. If people consume more energy content, their waste production via urination may also increase, causing increased USG independent of water intake. Increased USG quartiles were associated with decreased cardiovascular risks, but this association was paralleled by a decrease in age. Because, in this study, mean USG was higher in men than women, in good agreement with the known sex difference in urine concentration(Reference Perucca, Bouby and Valeix30), the effect of sex was also considered. In this way, the genuine effect of USG on IR was appropriately adjusted for the potential confounders. In this study, increased USG was independently associated with increased odds for high IR, even after adjusting for age, sex, habitual patterns, various cardiovascular risks, and plain water and energy intakes. Specifically, per 0·01 increase in USG was associated with 30·3 % increased odds for high IR in the multivariate logistic regression analysis. Using a generalised additive model plot, we also identified that increased USG showed a J-shaped association with IR and the suggested threshold in the AIC analysis was USG ≥ 1·030. In our study population, the proportion of participants with USG ≥ 1·030 was only 2·2 %. Therefore, some may argue that although the hazard of low hydration status on IR may be statistically significant, the clinical significance of this low hydration may be limited. However, our results should be interpreted in the context of risk–benefit considerations. Several of the previous approaches to improve insulin sensitivity involved the use of various medications(Reference Canivell, Mohaupt and Ackermann10,Reference Hiroyama, Aoyagi and Fujiwara19,Reference Enhorning and Malan31) , which are not free from inadvertent side effects(Reference Wannamethee, Welsh and Papacosta9,Reference DeFronzo, Fleming and Chen32,Reference Ahmed33) . Although lifestyle modification with exercise is a must for people with diabetes, it is very difficult to achieve and maintain(Reference Wahl, Scalzo and Regensteiner34). In contrast, low hydration status can be easily improved simply by drinking more water(Reference Lemetais, Melander and Vecchio35,Reference Enhorning, Tasevska and Roussel36) . Therefore, the potential benefit of improved hydration status on glucose metabolism should be tested in future clinical trials. Actually, a recent study showed that, in adults with high copeptin levels, an increased hydration induced not only a significant decline in copeptin (as expected) but also a significant decline in glycaemia(Reference Enhorning, Brunkwall and Tasevska37).

In the present study, we found that increased USG was associated with increased trunk fat and decreased leg fat, resulting in increased TLR, after adjusting for age, sex, habitual patterns, various cardiovascular risks and plain water and energy intakes. The relationship between USG and TLR was largely attributed to the consistent relationship between USG and leg fat. Although Chang et al. have suggested that increased urine osmolality was associated with a higher BMI after analysing 9528 US citizens(Reference Chang, Ravi and Plegue38), they did not analyse the effect of hydration status on regional adiposity any further. Therefore, our study is the first to demonstrate that hydration status is related to body fat distribution. The differential effects of low hydration status on body fat reminded us of different results of AVP on body fat according to the effects mediated by the different types of AVP receptors; the activation of V1a receptor resulted in fat accumulation, whereas activation of V1b receptor resulted in lipolysis(Reference Koshimizu, Nakamura and Egashira5). Therefore, the influence of AVP on trunk fat may preferentially depend on the V1a receptor, whereas that on leg fat may preferentially depend on the V1b receptor. A possible differential expression of AVP receptors according to fat depots should be explored in the future experimental studies, and in clinical investigations.

AVP receptors are widely expressed in metabolic tissues, including adipose tissue, pancreas islet cells and the H-P-A axis(Reference Koshimizu, Nakamura and Egashira5). Poor fat distribution is tightly associated with IR(Reference Shah, Murthy and Abbasi39Reference Lee, Son and Kim42). Therefore, it is highly suspected that the relationship of USG with poor fat distribution may be mediated via IR, and the effect of USG on IR may be mediated via poor fat distribution. In our mediation analysis, however, only 19·3 % of the effect of USG on fat distribution was mediated via IR, whereas only 13·0 % of the effect of USG on IR was mediated via poor fat distribution. These results may imply that direct effects of AVP via V1a and V1b receptors on target tissues are a more dominant pathway than the indirect effects via IR or poor fat distribution. Therefore, improving hydration status may be an easy and powerful strategy to improve public health by restoring insulin-glucose metabolism and fat distribution.

This study has several limitations. First, it is a cross-sectional study; therefore, a causal relationship cannot be established. However, the study hypotheses were constructed based on the known causality between AVP and metabolism of glucose and fat demonstrated in animal studies. We aimed to revise the focus on AVP, from its mechanism to its clinical context, to expand clinical utility. Second, USG was measured using the first morning urine sample after fast. USG from a single urine sample may have a potential error because it may not be representative of 24 h hydration status(Reference Perrier, Demazieres and Girard43). In addition, because of this spottiness, USG may simply mean tubular concentrating capacity, not usual hydration status. However, the concentration of the first morning urine is well correlated with that of 24 h urine(Reference Youhanna, Bankir and Jungers44). Furthermore, in KNHANES, urine was always collected in morning after fast, suggesting urine may be standardised in sampling time. Therefore, although incomplete, USG in the first morning urine can be used as hydration marker, particularly in large-scale populations. Unfortunately, other hydration indices, including plasma and urine osmolality, bioimpedance analysis, and copeptin levels, were not measured because the KNHANES was not primarily designed to evaluate the hydration status. Finally, the generalisability of our results is limited because the study was conducted in a single country with a single ethnicity.

Besides these limitations, our study has several strengths. First, it is based on a well-defined very large population cohort (over 14 000 subjects) with a relatively uniform ethnic origin. Second, an extremely detailed evaluation of the metabolic phenotype was available, based on the measurement of a large number of variables. Third, same protocol to collect information and specimens was applied throughout study period. Furthermore, blood and urine samples were tested in the same central laboratory. Finally, this is the first study which suggested that hydration status, presumably AVP status, may affect body fat deposition differently depending on fat depots (trunk v. leg), which is a novel study subject to be confirmed by future animal and clinical studies.

In conclusion, this study revealed that a low hydration status, a presumably elevated AVP condition, was associated with increased IR and poor fat distribution. Direct effect of low hydration status may be more dominant than indirect effect via IR or fat distribution. Future studies are necessary to confirm our results.

Acknowledgements

There was no funding for this study.

Authors’ contributions were: H. K. M. and S. W. L. came up with the research idea and designed the study; H. Y. K. and J. T. K. anaysed and interpreted the data; S. W. L. performed the statistical analysis; H. K. M. and S. W. L. wrote the manuscript; L. B. provided supervision and mentorship. Each author contributed important intellectual content during manuscript drafting or revision and accepts accountability for the overall work by ensuring that questions pertaining to the accuracy or integrity of any portion of the work were appropriately investigated and resolved.

The authors declare that there are no conflicts of interest.

References

Bankir, L (2001) Antidiuretic action of vasopressin: quantitative aspects and interaction between V1a and V2 receptor-mediated effects. Cardiovasc Res 51, 372390.CrossRefGoogle ScholarPubMed
Bouby, N, Ahloulay, M, Nsegbe, E, et al. (1996) Vasopressin increases glomerular filtration rate in conscious rats through its antidiuretic action. J Am Soc Nephrol 7, 842851.Google ScholarPubMed
Clark, WF, Sontrop, JM, Huang, SH, et al. (2016) Hydration and chronic kidney disease progression: a critical review of the evidence. Am J Nephrol 43, 281292.CrossRefGoogle Scholar
Min, HK, Sung, SA, Lee, SY, et al. (2019) Sub-morbid dehydration-associated glomerular hyperfiltration: an emerging reality? Kidney Res Clin Pract 38, 196204.CrossRefGoogle Scholar
Koshimizu, TA, Nakamura, K, Egashira, N, et al. (2012) Vasopressin V1a and V1b receptors: from molecules to physiological systems. Physiol Rev 92, 18131864.CrossRefGoogle ScholarPubMed
Muscogiuri, G, Barrea, L, Annunziata, G, et al. (2018) Water intake keeps type 2 diabetes away? Focus on copeptin. Endocrine 62, 292298.CrossRefGoogle ScholarPubMed
Abu-Basha, EA, Yibchok-Anun, S & Hsu, WH (2002) Glucose dependency of arginine vasopressin-induced insulin and glucagon release from the perfused rat pancreas. Metabolism 51, 11841190.CrossRefGoogle ScholarPubMed
Christ-Crain, M & Fenske, W (2016) Copeptin in the diagnosis of vasopressin-dependent disorders of fluid homeostasis. Nat Rev Endocrinol 12, 168176.CrossRefGoogle Scholar
Wannamethee, SG, Welsh, P, Papacosta, O, et al. (2015) Copeptin, insulin resistance, and risk of incident diabetes in older men. J Clin Endocrinol Metab 100, 33323339.CrossRefGoogle ScholarPubMed
Canivell, S, Mohaupt, M, Ackermann, D, et al. (2018) Copeptin and insulin resistance: effect modification by age and 11 beta-HSD2 activity in a population-based study. J Endocrinol Invest 41, 799808.CrossRefGoogle ScholarPubMed
Jensen, T, Bjornstad, P, Johnson, RJ, et al. (2019) Copeptin and estimated insulin sensitivity in adults with and without type 1 diabetes: the CACTI study. Can J Diabetes 43, 3439.CrossRefGoogle ScholarPubMed
Saleem, U, Khaleghi, M, Morgenthaler, NG, et al. (2009) Plasma carboxy-terminal provasopressin (copeptin): a novel marker of insulin resistance and metabolic syndrome. J Clin Endocrinol Metab 94, 25582564.CrossRefGoogle ScholarPubMed
Enhorning, S, Struck, J, Wirfalt, E, et al. (2011) Plasma copeptin, a unifying factor behind the metabolic syndrome. J Clin Endocrinol Metab 96, E1065E1072.CrossRefGoogle ScholarPubMed
Vintilă, M, Gheorghiu, ML, Caragheorgheopol, A, et al. (2016) Increased copeptin levels in metabolic syndrome from a Romanian population. J Med Life 9, 353357.Google ScholarPubMed
Enhorning, S, Wang, TJ, Nilsson, PM, et al. (2010) Plasma copeptin and the risk of diabetes mellitus. Circulation 121, 21022108.CrossRefGoogle ScholarPubMed
Abbasi, A, Corpeleijn, E, Meijer, E, et al. (2012) Sex differences in the association between plasma copeptin and incident type 2 diabetes: the Prevention of Renal and Vascular Endstage Disease (PREVEND) study. Diabetologia 55, 19631970.CrossRefGoogle ScholarPubMed
Enhorning, S, Bankir, L, Bouby, N, et al. (2013) Copeptin, a marker of vasopressin, in abdominal obesity, diabetes and microalbuminuria: the prospective Malmo Diet and Cancer Study cardiovascular cohort. Int J Obes 37, 598603.CrossRefGoogle Scholar
Roussel, R, El Boustany, R, Bouby, N, et al. (2016) Plasma copeptin, AVP gene variants, and incidence of type 2 diabetes in a cohort from the community. J Clin Endocrinol Metab 101, 24322439.CrossRefGoogle Scholar
Hiroyama, M, Aoyagi, T, Fujiwara, Y, et al. (2007) Hypermetabolism of fat in V1a vasopressin receptor knockout mice. Mol Endocrinol 21, 247258.CrossRefGoogle ScholarPubMed
Hiroyama, M, Fujiwara, Y, Nakamura, K, et al. (2009) Altered lipid metabolism in vasopressin V1B receptor-deficient mice. Eur J Pharmacol 602, 455461.CrossRefGoogle ScholarPubMed
Fujiwara, Y, Hiroyama, M, Sanbe, A, et al. (2007) Insulin hypersensitivity in mice lacking the V1b vasopressin receptor. J Physiol 584, 235244.CrossRefGoogle ScholarPubMed
Garrett, DC, Rae, N, Fletcher, JR, et al. (2018) Engineering approaches to assessing hydration status. IEEE Rev Biomed Eng 11, 233248.CrossRefGoogle ScholarPubMed
Kavouras, SA (2002) Assessing hydration status. Curr Opin Clin Nutr Metab Care 5, 519524.CrossRefGoogle ScholarPubMed
Baron, S, Courbebaisse, M, Lepicard, EM, et al. (2015) Assessment of hydration status in a large population. Br J Nutr 113, 147158.CrossRefGoogle ScholarPubMed
Perrier, E, Vergne, S, Klein, A, et al. (2013) Hydration biomarkers in free-living adults with different levels of habitual fluid consumption. Br J Nutr 109, 16781687.CrossRefGoogle ScholarPubMed
Matthews, DR, Hosker, JP, Rudenski, AS, et al. (1985) Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 28, 412419.CrossRefGoogle ScholarPubMed
Ito, H, Nakasuga, K, Ohshima, A, et al. (2003) Detection of cardiovascular risk factors by indices of obesity obtained from anthropometry and dual-energy X-ray absorptiometry in Japanese individuals. Int J Obes Relat Metab Disord 27, 232237.CrossRefGoogle ScholarPubMed
Levey, AS, Stevens, LA, Schmid, CH, et al. (2009) A new equation to estimate glomerular filtration rate. Ann Intern Med 150, 604612.CrossRefGoogle ScholarPubMed
Alberti, KG, Zimmet, P & Shaw, J (2006) Metabolic syndrome – a new world-wide definition. A consensus statement from the International Diabetes Federation. Diabet Med 23, 469480.CrossRefGoogle ScholarPubMed
Perucca, J, Bouby, N, Valeix, P, et al. (2007) Sex difference in urine concentration across differing ages, sodium intake, and level of kidney disease. Am J Physiol Regul Integr Comp Physiol 292, R700R705.CrossRefGoogle ScholarPubMed
Enhorning, S & Malan, L (2019) Copeptin relates to a fatty liver and measures of obesity in a South African population with mixed ethnicities. Endocrine 65, 304311.CrossRefGoogle Scholar
DeFronzo, R, Fleming, GA, Chen, K, et al. (2016) Metformin-associated lactic acidosis: current perspectives on causes and risk. Metabolism 65, 2029.Google ScholarPubMed
Ahmed, MA (2016) Metformin and vitamin B12 deficiency: where do we stand? J Pharm Pharm Sci 19, 382398.CrossRefGoogle ScholarPubMed
Wahl, MP, Scalzo, RL, Regensteiner, JG, et al. (2018) Mechanisms of aerobic exercise impairment in diabetes: a narrative review. Front Endocrinol 9, 181.CrossRefGoogle ScholarPubMed
Lemetais, G, Melander, O, Vecchio, M, et al. (2018) Effect of increased water intake on plasma copeptin in healthy adults. Eur J Nutr 57, 18831890.CrossRefGoogle ScholarPubMed
Enhorning, S, Tasevska, I, Roussel, R, et al. (2019) Effects of hydration on plasma copeptin, glycemia and gluco-regulatory hormones: a water intervention in humans. Eur J Nutr 58, 315324.CrossRefGoogle ScholarPubMed
Enhorning, S, Brunkwall, L, Tasevska, I, et al. (2019) Water supplementation reduces copeptin and plasma glucose in adults with high copeptin: The H2O Metabolism Pilot Study. J Clin Endocrinol Metab 104, 19171925.Google ScholarPubMed
Chang, T, Ravi, N, Plegue, MA, et al. (2016) Inadequate hydration, BMI, and obesity among US adults: NHANES 2009–2012. Ann Fam Med 14, 320324.CrossRefGoogle ScholarPubMed
Shah, RV, Murthy, VL, Abbasi, SA, et al. (2014) Visceral adiposity and the risk of metabolic syndrome across body mass index: the MESA Study. JACC Cardiovasc Imaging 7, 12211235.CrossRefGoogle ScholarPubMed
Patel, P & Abate, N (2013) Body fat distribution and insulin resistance. Nutrients 5, 20192027.CrossRefGoogle ScholarPubMed
Neeland, IJ, Turer, AT, Ayers, CR, et al. (2012) Dysfunctional adiposity and the risk of prediabetes and type 2 diabetes in obese adults. JAMA 308, 11501159.CrossRefGoogle ScholarPubMed
Lee, SW, Son, JY, Kim, JM, et al. (2018) Body fat distribution is more predictive of all-cause mortality than overall adiposity. Diabetes Obes Metab 20, 141147.CrossRefGoogle ScholarPubMed
Perrier, E, Demazieres, A, Girard, N, et al. (2013) Circadian variation and responsiveness of hydration biomarkers to changes in daily water intake. Eur J Appl Physiol 113, 21432151.CrossRefGoogle ScholarPubMed
Youhanna, S, Bankir, L, Jungers, P, et al. (2017) Validation of surrogates of urine osmolality in population studies. Am J Nephrol 46, 2636.Google ScholarPubMed
Figure 0

Table 1. Clinical characteristics of the study population according to sex-specific urine specific gravity (USG) quartiles§(Mean values and standard deviations for continuous variables; numbers and percentages for categorical variables; medians and interquartile ranges)

Figure 1

Table 2. Association between urine specific gravity (USG) and high insulin resistance*(Odds ratios and 95 % confidence intervals)

Figure 2

Fig. 1. Relationship between urine specific gravity and insulin resistance. The dashed line indicates 95 % CI for value of the smoothed logarithmic homoeostasis model assessment of insulin resistance (HOMA-IR) using multivariate generalised additive model analysis after adjusting for age, sex, smoking and drinking status, suicide thought, previous CVD, proteinuria, estimated glomerular filtration rate, systolic and diastolic blood pressure, waist circumference, BMI, HDL-cholesterol, TAG, lean body mass, trunk fat, leg fat, arm fat, head fat, high- and mid-intensity exercise, plain water and energy intake, and urine Na as covariates. AIC, Akaike’s information criterion.

Figure 3

Table 3. Relationship between urine specific gravity and four fat depots*(β Values and 95 % confidence intervals)

Figure 4

Fig. 2. Mediation analysis of insulin resistance (IR) for the relationship between urine specific gravity (USG) and fat distribution. IR was estimated by homoeostasis model assessment (HOMA). Outcome was trunk:leg fat ratio (TLR), treatment was USG, and mediator was logarithmic HOMA-IR. Overall effects of treatment on outcome were adjusted by mediator and pre-treatment covariates using multivariate linear regression analysis. Mediator was modelled by treatment and pre-treatment covariates using multivariate linear regression analysis.

Figure 5

Fig. 3. Mediation analysis of fat distribution for the relationship between urine specific gravity (USG) and insulin resistance. Outcome was logarithmic homoeostasis model assessment of insulin resistance (HOMA-IR), treatment was USG, and mediator was trunk:leg fat ratio (TLR). Overall effects of treatment on outcome were adjusted by mediator and pre-treatment covariates using multivariate linear regression analysis. Mediator was modelled by treatment and pre-treatment covariates using multivariate linear regression analysis.