Novel Strategies for Assessing Associations Between Selenium Biomarkers and Cardiometabolic Risk Factors: Concentration, Visit-to-Visit Variability, or Individual Mean? Evidence From a Repeated-Measures Study of Older Adults With High Selenium

Introduction

Cardiovascular disease (CVD) is a leading cause of death in East Asia and the world. In 2019, an estimated 18.6 million [95% confidence interval (CI): 17.1–19.7 million] CVD-related deaths occurred globally (1). Cardiometabolic risk factors, including dyslipidaemia, hypertension, diabetes, and obesity, are the primary causes of CVD (2, 3). Because of the antioxidant effects of selenium, selenium supplements and high-selenium foods are widely consumed to prevent the development of cardiometabolic risk factors (4). However, evidence on the influence of selenium on cardiometabolic risk factors remains inconclusive, and selenium may even increase cardiometabolic risk (5–7). Thus, whether selenium should be promoted for its cardiometabolic protective effects remains uncertain.

A cross-sectional study of healthy adults in two Chinese counties with different selenium intake habits indicated that serum selenium was positively correlated with serum glucose in those with selenium deficiency (median: 58 μg/L). However, no significant relationship between serum selenium and glucose was observed in those without such deficiency (median: 103 μg/L) (8). Additionally, a cross-sectional analysis using data from the National Health and Nutrition Examination Survey (NHANES) reported a U-shaped association between plasma selenium (mean: 137.1 μg/L) and the likelihood of diabetes, with the lowest risk noted for a concentration of ~122 μg/L (9). An experimental study discovered that the dose-dependent effects of selenium range from antioxidation and anti-inflammation to the promotion of oxidative stress and insulin resistance (10). Therefore, two crucial factors must be emphasized regarding the effect of selenium on cardiometabolic risk: the U-shaped associations between selenium levels and cardiometabolic risk factors and the broad individual variation in selenium concentrations. According to the study of Reference Man by the International Commission on Radiological Protection, a mean whole blood selenium exceeding 130 μg/L is defined as high selenium (11). However, most relevant epidemiological studies have been conducted on populations with moderate or deficient selenium levels, with means of blood selenium concentrations ranging from 47 to 127.5 μg/L (12–14). Few epidemiological studies have been conducted in populations with high selenium (>130 μg/L), limiting the understanding of the role of high selenium in the development of cardiometabolic risk factors (15, 16). Selenium-rich soil with concentrations exceeding 0.4 mg/kg has been discovered in areas across Beijing (17), providing an excellent opportunity to investigate the association between high selenium levels and cardiometabolic risk in a real-world setting. Compared with young adults, fewer older adults (aged 50 years or older) smoke or consume alcohol (18, 19). In addition, the metabolic and physiological functions, including excretion function, of older adults are gradually impaired during the aging process (20, 21). Therefore, compared with other age groups, older adults have higher selenium levels and more cases of high selenium (22–24). Moreover, aging adults experience multisystem functional impairment and increasing susceptibility to multiple chronic diseases (25). Therefore, examining older adults in a study of the cardiometabolic health effects of selenium is warranted.

Selenium assessment is mainly based on selenium concentrations in blood and urine. Increasing evidence suggests inconsistent associations among urinary selenium, circulating selenium (whole blood selenium, plasma selenium, or serum selenium), and cardiometabolic risk factors. A case–control study conducted in Wuhan, China, indicated no association between urinary selenium (mean: 20.47 μg/g) and blood pressure (BP) (26). However, another case–control study conducted in Wuhan demonstrated a U-shaped association between plasma selenium (median: 92.66 μg/L) and hypertension (12). Therefore, the associations of multiple selenium measures (e.g., whole blood and urinary selenium) with cardiometabolic risk factors should be examined simultaneously. Three cross-sectional analyses using data from the NHANES reported diverse associations of serum selenium with blood pressure and glucose and lipid profiles. U-shaped associations of serum selenium (mean: 137.1 μg/L) with systolic BP (SBP) and pulse pressure were observed (27). In addition, high serum selenium (mean: 137.1 μg/L) was associated with high glucose levels (9). Another study discovered that total cholesterol and triglyceride levels increased with serum selenium (median: 192.99 μg/L); furthermore, low-density lipoprotein cholesterol (LDL-C) was non-linearly associated with serum selenium, but high-density lipoprotein cholesterol (HDL-C) did not vary with serum selenium (5). This evidence suggests that multiple cardiometabolic risk factors should be measured to comprehensively characterize the associations of selenium levels with cardiometabolic health and evaluate the consistency of results to draw robust conclusions.

Changes in diet, lifestyle, or daily activities dynamically affect biomarkers and physiological parameters and these variations influence health endpoints (28). Visit-to-visit variability (VVV) in BP and glucose and lipid profiles has been associated with CVD and mortality (29–31). Selenium homeostasis is essential for a wide range of cellular functions, such as modulation of the cell cycle and apoptosis, redox balancing, and protein and DNA synthesis (32). Therefore, investigating the VVV or individual mean (IM) of selenium levels may be beneficial, while exploring the associations between selenium and cardiometabolic risk factors. Such considerations may elucidate a crucial means of reducing misclassifications in exposure assessment and further understanding the associations between selenium and cardiometabolic health. Because VVV and IM data cannot be calculated in a traditional cross-sectional study or a cohort study measuring only baseline exposure, a repeated-measures study is required.

Using a three-wave repeated-measures study of 201 older adults with high selenium residing in Beijing, China, from 2016 through 2018, we explored the associations between selenium biomarkers (i.e., whole blood and urinary selenium concentration) and cardiometabolic risk factors. On the basis of the selenium concentrations measured during three clinic visits, selenium homeostasis factors, namely VVV and IM, were considered. Furthermore, we investigated the associations of the VVV and IM of selenium levels with cardiometabolic risk factors.

Methods

Study Setting and Population

The study was conducted in Beijing, which is in the northern North China Plain. As previously mentioned, selenium-rich soil has been discovered across the Beijing area from north to south. We selected five communities (Qian Nantai, Liu Hegou, Dongcheng, Chaoyang, and Fangshan) from four regions of Beijing from north to south to serve as the study settings. Among five communities, Qian Nantai, Liu Hegou are located in rural area. Dongcheng, Chaoyang and Fangshan are located in urban area. Further details regarding the study setting and design are presented in Figure 1.

FIGURE 1

Figure 1. (A) Five study locations within Beijing. (B) Overview of study design and data collection.

The inclusion and exclusion criteria were the same as those described in a previous study (33). We attempted to sample a population of older adults with stable selenium statuses and the ability to complete three clinic visits by recruiting individuals aged 50 years or older who had lived in the community for more than 5 years and were unlikely to leave during the study period. Those who were unable to complete a questionnaire survey, had received a malignant tumor diagnosis, or had CVD or liver disease were excluded. In addition, a questionnaire inquired into the participants' dietary supplement intake, and the participants did not consume any dietary selenium supplements. All three measurements were conducted in the winter to control for confounding seasonal factors. After the first visit in November or December 2016, 201 participants were enrolled. Two follow-up surveys were conducted in November 2017 and January 2018. Ultimately, 83% (n = 167) of the study participants completed all three visits, and 17% (n = 34) completed only two visits. In total, 569 observations were analyzed.

The study protocol was reviewed and approved by the Institutional Ethics Committee of the Institute of Basic Medicine at the Chinese Academy of Medical Sciences. Each participant provided written consent before participating.

Measurement of Whole Blood Selenium and Urinary Selenium

Species Sampling

The participants fasted overnight (≥8 h) before each examination. A qualified nurse obtained peripheral blood samples from all participants at the same time of day (8:00–9:00 am) to control for variations in circadian rhythm. Whole blood and serum were collected for whole blood selenium concentration measurements and serum biochemical analyses. First morning urine samples (10 mL) were collected from the participants in trace element—free containers. All urine and blood samples were stored at −80°C for subsequent analysis.

Measurements of Blood and Urine Samples

Before analysis, the whole blood samples were diluted (1:10 ratio) in a digestive solution composed of 0.5% nitric acid and 0.01% Triton X-100. The urine samples were diluted (1:10) in a digestive solution composed of 0.5% nitric acid and 0.02% Triton X-100. Both the whole blood selenium and urinary selenium concentrations were measured through inductively coupled plasma mass spectrometry (Nexion 300D, PerkinElmer, USA). In addition, we applied quality control protocols to ensure the accuracy of our analyses. First, all the samples were measured three times. Second, standard reference materials were measured every 20 samples to ensure that the whole blood and urinary selenium measurements were consistent with certified concentrations before subsequent samples were assayed. The limits of detection (LODs) for whole blood and urinary selenium were determined to be 2.15 and 0.14 μg/L, respectively. The concentrations in all samples were higher than the LODs.

Adjustment of Urinary Selenium Concentration

In epidemiological studies of environmental contaminants measured in urine, adjustment for creatinine to control for the effect of urine dilution is common, but the optimal approach is debated (34). Traditionally, urinary biomarker concentrations have been standardized through division by the urinary creatinine concentration, in accordance with the assumption that creatinine excretion is approximately constant across individuals and time (35, 36). However, increasing evidence suggests that creatinine levels can vary with time and individual characteristics, including sex, race, age, and body mass index (BMI); thus, conventional creatinine adjustment may yield misleading conclusions. Therefore, we employed a three-step approach, enabling systematic individual differences in long-term creatinine levels to be used to estimate covariate-adjusted, standardized urinary selenium concentrations (37). Urinary creatinine concentrations were measured using the sarcosine oxidase method. First, natural log–transformed creatinine was regressed on factors that can affect urine dilution. We included age, sex, race, weight, and height as covariates in the regression model (38). The measured creatinine concentration was divided by the fitted creatinine value obtained from the regression model to obtain a ratio representing the covariate-independent residual effect of hydration on creatinine. We finally standardized the urinary selenium concentration by dividing it by this ratio. This method specifically controls the covariate-independent, short-term multiplicative effect of hydration on urine dilution, avoids potential collider bias resulting from the use of urinary creatinine or specific gravity in adjusting for urine dilution, and increases the accuracy of exposure assessment (38, 39).

Calculation of VVV and IM of Selenium Biomarkers

The VVV and IM of the selenium biomarkers were calculated to estimate the fluctuation and homeostasis of the biomarkers across three visits. The IM of the selenium biomarkers for each participant was the mean selenium concentrations from the three visits. No gold standard regarding the appropriate indicators of the VVV of selenium biomarkers exists (40). According to the concentrations of the selenium biomarkers across the three visits, we calculated four indices of VVV, namely the standard deviation (SD), coefficient of variation (CV), average real variability (ARV), and variability independent of the mean (VIM), to comprehensively assess variability. SD and CV are simple and informative indices of VVV (41). SD was used to indicate the general variation between the selenium biomarker levels and the mean values of the participants. However, the order of measurement is not considered in SD calculation (42). ARV was used to overcome the limitations of SD by accounting for the order of measurements when quantifying VVV (29). However, ARV was correlated with the selenium concentrations. The CVs of the selenium biomarkers across the three visits were computed from the mean and SD. Therefore, the CVs might be correlated with the mean selenium concentrations (43). Hence, the VIM was calculated through logarithmic curve fitting to eliminate the correlation with the biomarkers, ensuring that we could evaluate the impact of the VVV separately from the effect of the selenium biomarker levels themselves (29). The formulae for the four VVV indices are as follows:

$\begin{array}{ll}\begin{array}{ccc}SD& =& \frac{\sqrt{{\displaystyle {\sum }_{k=1}^n{(S{e}_k-IM)}^2}}}{n-1}\end{array}& (1)\end{array}$

$\begin{array}{ll}CV=\frac{SD}{IM}& (2)\end{array}$

$\begin{array}{ll}ARV=\frac{1}{n-1}{\displaystyle \sum _{k=1}^{n-1}|S{e}_{k+1}-S{e}_k|}& (3)\end{array}$

$\begin{array}{ll}\ln (S{D}_i)={\beta }_0+{\beta }_1\ln (IM)& (\mathrm{4-1})\end{array}$

$\begin{array}{ll}VI{M}_i= {\overline{M}}^{{\beta }_1}\frac{S{D}_i}{I{M}^{{\beta }_1}}& (\mathrm{4-2})\end{array}$

where Se denotes the selenium concentration; IM denotes the individual mean of the concentration; n is the total number of measurements (n = 2 or 3); k indicates the visit number (k = 1, 2, or 3); Se_K is the selenium concentration recorded during visit k; $\overline{M}$ is the overall mean of the selenium concentration for the study population; SD_i is the SD of the selenium measurements for individual i; and VIM_i is the VIM of the selenium measurements for individual i. In total, we calculated eight selenium VVV indicators: the SD, CV, ARV, and VIM of whole blood and urinary selenium (denoted as SD-WBSe, CV-WBSe, ARV-WBSe, VIM-WBSe, SD-USe, CV-USe, ARV-USe, and VIM-USe, respectively).

Outcome Assessment

Physical Measurement

The participants' upper arm BP was measured using calibrated mercury sphygmomanometer at least three times after the participants had rested for 15 min in a room. SBP and diastolic BP (DBP) were obtained at the appearance and disappearance of Korotkoff sounds, respectively. The interval between measurements was required to exceed 2 min. Additional measurements were taken if the differences among measurements exceeded 5 mmHg. The mean of the final two measurements was calculated for subsequent analysis. Mean arterial pressure (MAP) was calculated because MAP is associated with the risks of adverse cardiovascular events (44). The following formula was used to estimate MAP:

$\begin{array}{ll}MAP=(2\times DBP+SBP)/3& (5)\end{array}$

Height was measured using a verified stadiometer, with each participant in a standing position with shoes removed, shoulders relaxed, head facing forwards, and back facing the wall. Weight was measured with the participants wearing little clothing, by using a certified body composition analyser. The participants' weights in kilograms were divided by their heights in meters squared to obtain their BMIs. For waistline and hipline measurements, each participant was asked to wear minimal clothing and stand with their feet close together, arms at their side, and body weight evenly distributed. Waistline was measured in centimeters at the midpoint between the lower margin of the last palpable rib and the top of the iliac crest using stretch-resistant tape. Hipline measured in centimeters around the widest portion of the buttocks, with the tape parallel to the floor (45). The waistline measurements were divided by the hipline measurements to obtain waist-hip ratios. All the measurements were performed twice and averaged for reliability and accuracy.

Clinical Laboratory Examination

The serum samples were sent to the Clinical Laboratory Center of Peking Union Medical College Hospital for biochemical testing. Lipid (total cholesterol, triglycerides, HDL-C, and LDL-C) and glucose, uric acid, and high-sensitivity C-reactive protein (hs-CRP) levels were measured using a Beckman Coulter analyser (AU2700, Beckman Coulter, Brea, CA, USA).

Metabolic Syndrome Score

The metabolic syndrome score (MSS), calculated from various factors related to metabolic syndrome, is an accurate indicator of cardiometabolic risk (46). The MSS is correlated with long-term risk for diabetes and CVD and is widely employed in clinical practice to identify patients with high such risks (47, 48). To account for potential sex-based differences in the MSSs, we computed the sex-specific MSS (MSS-sex) for each participant to comprehensively evaluate their cardiometabolic health (49). MSS-sex was calculated by summary of sex-specific standardized Z-scores for each of five components including MAP, HDL-C, triglyceride, waistline, and glucose, generated from the population being studied (50).

Ascertainment of Covariates

Questionnaires were used to obtain the demographic characteristics, including age (years), sex (male or female), race (Han or others), and educational attainment (illiteracy or primary school, or junior school, or senior high school and above); lifestyle characteristics, such as cigarette smoking (current, former and never) and alcohol consumption (current, former and never) habits; and health status. The participants' communities were coded as categorical variables. Additionally, because this was a longitudinal study, controlling for temporal confounders was necessary. Therefore, sampling months were recorded and adjusted for as potential confounders.

In accordance with China's medical guidelines for the prevention and treatment of dyslipidaemia in adults (51), we considered hypercholesterolaemia, hypertriglyceridemia, hyperalphalipoproteinaemia, hyperbetalipoproteinaemia, or the use of any antidyslipidaemic agents (e.g., statins) , or both to indicate dyslipidaemia. Hypercholesterolaemia was defined as total cholesterol ≥6.2 mmol/L, hypertriglyceridemia was defined as triglycerides ≥2.3 mmol/L, hyperalphalipoproteinaemia was defined as HDL-C ≤ 1.0 mmol/L, and hyperbetalipoproteinaemia was defined as LDL-C ≥ 4.1 mmol/L. Hypertension was indicated by any antihypertensive medication prescription or BP ≥140/90 mmHg measured under standardized clinical conditions (52). Diabetes was defined as fasting glucose ≥7.0 mmol/L, 2-h glucose ≥11.1 mmol/L, or current antidiabetic medication use (both insulin or oral antidiabetic drugs) (53).

Dietary intake information was collected using a food frequency questionnaire (FFQ), which inquired into the frequency and amount of food consumed over the previous 30 days. The FFQ comprised 81 food-related items, and its reproducibility and validity were verified in a prior study (54). The daily total energy intake (kcal/day) was determined from the FFQ data and Chinese Tables of Food Composition (55).

Statistical Analysis

Statistical Description

The participants' general characteristics, demographic features, health status, VVV and IM of selenium biomarkers, and cardiometabolic risk factors were documented. The categorical variables are expressed as numbers (percentages), and the continuous variables are expressed as means ± SDs. To comprehensively profile the distribution of the selenium biomarkers, geometric means and percentiles for the selenium biomarkers were calculated. In addition, distribution characteristics of selenium biomarkers by urban and rural areas were described. Unadjusted linear mixed-effects regression models of the selenium biomarkers by area (urban area, rural area) were developed to assess difference of selenium biomarkers concentrations between two areas.

The variability in the selenium biomarkers and cardiometabolic risk factors during the study period was also calculated. We calculated intraclass correlation coefficients (ICCs) of the selenium biomarkers, defined as the ratio of the interindividual variance to the total variance, and corresponding 95% CIs. The ICCs were used to evaluate the reproducibility of the selenium biomarkers across the three visits as follows: 0.75 ≤ ICC ≤ 1.00 suggested excellent reproducibility, 0.40 ≤ ICC <0.75 suggested fair to good reproducibility, and 0.00 ≤ ICC <0.40 suggested poor reproducibility (56, 57). The distributions of the selenium biomarker values across the three visits were recorded and represented using violin plots. Spearman correlation coefficients between each pair of biomarkers measured at each visit were calculated. Unadjusted linear mixed-effects regression models of the selenium biomarkers and cardiometabolic risk factors by visit (modeled as a nominal categorical variable) were developed to assess global significance for each visit. The variability of the selenium biomarkers and cardiometabolic risk factors was helpful for both constructing models and explaining results in the statistical analyses.

Associations of Selenium Biomarkers With Cardiometabolic Risk Factors

Covariate Selection

We applied a directed acyclic graph to determine covariates adjusted in the models a prior, which was the minimum adjustment sets selected by graphical criteria such as the so-called “back-door” (58, 59). The covariates included were age, sex, race, educational attainment, cigarette smoking habits, alcohol consumption habits, BMI, sampling month, community, and daily total energy intake.

Exposure–Response Curves

The generalized additive mixed models (GAMMs) with embedded restricted cubic spline smoothers with three degrees of freedom were conducted to plotted exposure–response curves for the selenium biomarkers and cardiometabolic risk factors. The results of GAMMs were summarized to provide more numeric information about associations between selenium biomarkers and cardiometabolic risk factors. The linear and non-linear trends illustrated by the curves also helped guide the development and interpretation of the main analysis.

Main Analysis

To account for potentially skewed distributions, the concentrations of whole blood and urinary selenium were natural log–transformed before statistical analysis. Regarding outcomes, the levels of the four lipid indicators and of glucose and uric acid were natural log–transformed to reduce skewness in further analysis. Associations between the selenium biomarkers and cardiometabolic risk factors were assessed using linear mixed-effects regression models with participant-specific random intercepts to account for within-participant correlation in this repeated-measures study. For each model, one selenium biomarker was regressed on one cardiometabolic risk factor at the same study visit. According to the distributions of the selenium biomarkers and the non-linear and linear associations between them and the cardiometabolic risk factors, the selenium concentrations were included in the models as natural log–transformed continuous variables or categorical variables indicating tertiles, respectively. Linear trend p-values were derived by modeling the median of each selenium biomarker tertile as a continuous variable in the adjusted models. To address the multiple testing problem, in which the more inferences are made, the more likely erroneous inferences become, and reduce the probability of type I error, we adjusted the raw p-values on the basis of the false discovery rate (60).

Stratified Analysis

To account for the residual confounding and effect modification of the main analysis, we conducted analyses stratified by area (urban area or rural area), comorbidities (dyslipidaemia, diabetes, and hypertension), demographic characteristics [age (<65 or ≥65 years) and sex (male or female)], and BMI (<24 or ≥24 kg/m²), while simultaneously controlling for the same covariates as in the main analysis. The differences between strata were examined by estimating values and 95% CIs as follows (61).

$\begin{array}{ll}(\beta 1-{\beta }_2)\pm 1.96\times \sqrt{{(S{E}_1)}^2+{(S{E}_2)}^2}& (6)\end{array}$

where β₁ and β₂ are the effect estimates attributed to each subgroup or stratum (e.g., the effects for participants with and without dyslipidaemia, respectively) and SE₁ and SE₂ are the corresponding standard errors. In the stratified analysis, the selenium concentrations were included in the models as natural log–transformed continuous variables.

Sensitivity Analysis

We conducted several sensitivity analyses to assess the robustness of the findings in the main analysis. First, to address any residual confounding attributed to chronic disease status, we additionally adjusted for history of hypertension, diabetes, or dyslipidaemia in the same model of main analysis. Second, blood selenium concentrations reportedly differ significantly with BMI (62), and a higher BMI increases the risk of cardiometabolic disease (63). Therefore, we repeated our analysis while excluding BMI from the covariates to elucidate the impact of BMI on the associations. Third, hs-CRP, an indicator of systemic inflammation, may also influence both selenium levels and cardiometabolic risk (64, 65). Therefore, we additionally adjusted for hs-CRP in the sensitivity analysis to explore its influence on the associations. In the sensitivity analysis, the selenium concentrations were incorporated into the models as natural log–transformed continuous variables.

Associations of VVV and IM of Selenium Biomarkers With Cardiometabolic Risk Factors

The linear mixed-effects regression models employed in the main analysis were fitted to assess the associations of the VVV and IM of the selenium biomarkers with the cardiometabolic risk factors. A single-factor models (Model 1) respectively including individual indicator of VVV or IM in selenium biomarkers were used to assess the associations of the VVV or IM of the selenium biomarkers with the cardiovascular risk factors. Furthermore, mutually adjusted models (Model 2) considering both the VVV and IM of the selenium biomarkers were used to evaluate whether the effects of VVV in selenium biomarkers on cardiovascular risk factors are independent of IM. Both models were adjusted for the same covariates as described for the main analysis.

All statistical tests were two sided, and significance was set at p < 0.05. All analyses were conducted using the lmerTest, gamm4, splines, corrplot, dagitty, and ggplot2 packages in R (version 4.1.0).

Quality Assurance and Control

Standard operating procedures were followed throughout the research project. The interviewers, nurses, and physicians were all trained prior to the fieldwork. To ensure that the questionnaires collected personal information accurately, the principal investigators surveyed the interviewers with the same questionnaires to measure consistency. When each visit commenced, all the related data were recorded by qualified and trained staff (i.e., physicians, registered nurses, or trained interviewers). Data cleaning, quality assessment, and processing were performed independently by two qualified statistical analysts after completion of the surveys.

Results

Participant Characteristics

A summary of the participants' demographic data is presented in Table 1. Among the 201 participants included in the study, 125 (62.2 %) were female. The mean (SD) age and BMI of the participants were 64.60 (9.10) years and 25.10 (3.50) kg/m², respectively. The distributions of whole blood and urinary selenium concentration across the study period are detailed in Table 2. The geometric means (geometric SDs) of urinary and whole blood selenium were 18.00 (1.87) and 134.30 (1.19) μg/L, respectively. The ICCs (95% CIs) for whole blood and urinary selenium were 0.54 (0.44, 0.61) and 0.12 (0.02, 0.22), respectively, indicating good reproducibility of the whole blood measurements and poor reproducibility of the urinary measurements (Table 2). The characteristics of the cardiometabolic risk factors and the VVV and IM of the selenium biomarkers of the participants are presented in Table 3. The means (SDs) of the SD-WBSe, ARV-WBSe, and VIM-WBSe in population were 14.00 (8.70), 17.10 (12.90), 14.00 (8.10) μg/L, respectively. By comparison, the means (SD) of the SD-USe, ARV-USe, and VIM-USe in population were lower at 11.90 (12.90), 14.30 (15.80), and 10.80 (5.00) μg/L, respectively. By contrast, the mean (SD) of the CV-USe was higher than CV-WBSe in population [0.48 (0.26) vs. 0.10 (0.06)]. Furthermore, the mean (SD) of IM-WBSe and IM-USe in population were 136.00 (20.20) and 22.30 (11.60) μg/L, respectively.

TABLE 1

Table 1. Demographic characteristics and health status of participants throughout the study period.

TABLE 2

Table 2. Limits of detection, distribution characteristics, and reproducibility of the selenium biomarkers.

TABLE 3

Table 3. Summary of cardiometabolic risk factors and selenium homeostasis of participants.

The variations in the selenium biomarkers and cardiometabolic risk factors across the study period were also assessed, revealing three key findings: First, the global significance test revealed that whole blood selenium and most cardiometabolic risk factors (except SBP) significantly differed across the three visits (Supplementary Tables 1, 2). The difference of selenium biomarkers concentrations between urban and rural area was significant (Supplementary Table 3). Second, the overall correlation between the urinary selenium concentrations across three visits was negligible to moderate (r_s: −0.09 to 0.56). With respect to whole blood selenium, the moderate positive correlations were found between whole blood selenium levels across three visits (r_s: 0.61 to 0.65) which was consistent with better reproducibility of whole blood selenium compared to urinary selenium. Furthermore, there were negligible to low positive correlations between whole blood selenium concentrations and urinary selenium concentrations across three visits (r_s: 0.07 to 0.43; Supplementary Figure 2). Final, greater variance in the distribution of whole blood selenium than that of urinary selenium was also observed, as indicated in Supplementary Figure 3.

Associations of Selenium Biomarker Concentrations With Cardiometabolic Risk Factors

Exposure–Response Curves

The exposure–response curves of the selenium biomarkers and cardiometabolic risk factors obtained from the GAMMs are presented in Supplementary Figures 4, 5. Whole blood selenium was linearly associated with total cholesterol, triglycerides, HDL-C, uric acid, SBP, waistline, hipline and MSS-sex but non-linearly associated with LDL-C, glucose, DBP, MAP, and waist-hip ratio (Supplementary Figure 4). Consist with the exposure–response curves, estimated degrees of freedom (EDF) of LDL-C, glucose, DBP, MAP, and waist-hip ratio were 2.16, 2.14, 1.55, 1.46, and 2.22 which were significantly higher than 1 (Supplementary Table 4). Urinary selenium was linearly associated with cardiometabolic risk factors, except for waistline (Supplementary Figure 5). Similarly, EDF of waistline were 1.62 (Supplementary Table 4). Because of the diverse linear and non-linear relations of the selenium biomarkers with the various cardiometabolic risk factors, the selenium biomarkers were respectively considered continuous variables and categorized into tertiles in the linear mixed-effects regression models.

Main Analysis

In this three-wave repeated-measures study of 201 older adults with high selenium levels, we discovered that selenium levels were significantly associated with cardiometabolic risk factors after correction for multiple comparisons and adjustment for age, sex, race, educational attainment, cigarette smoking habits, alcohol consumption habits, BMI, sampling month, community, and total daily energy intake (Table 4). Primarily positive associations were observed with whole blood selenium: whole blood selenium was positively associated with total cholesterol (0.22, 95% CI: 0.12, 0.33), LDL-C (0.28, 95% CI: 0.13, 0.42), glucose (0.22, 95% CI: 0.10, 0.34), and uric acid (0.16, 95% CI: 0.04, 0.28) levels. These positive associations remained significant even with whole blood selenium concentration treated as a categorical variable (p trend < 0.01 for total cholesterol, LDL-C, and glucose; p trend = 0.02 for uric acid). Similarly, the results of GAMMs showed that associations between whole blood selenium and total cholesterol, LDL-C, glucose and uric acid were significant.

TABLE 4

Table 4. Associations of selenium biomarkers with cardiometabolic risk factors, analyzed using linear mixed-effects regression models.

For urinary selenium, only negative associations were observed and only with BP-related factors. Urinary selenium was negatively associated with SBP (−3.37, 95% CI: −5.68, −1.01), DBP (−1.69, 95% CI: −3.00, −0.39), and MAP (−2.32, 95% CI: −3.76, −0.88). However, these negative associations were non-significant when urinary selenium concentration was treated as a categorical variable (p trend = 0.08, 0.14, and 0.06 for SBP, DBP, and MAP, respectively). The results of GAMMs showed that associations between urinary selenium and BP indicators were significant.

Stratified Analyses

We conducted analyses stratified by area, comorbidities (dyslipidaemia, diabetes, and hypertension status) as well as by age, sex, and BMI (Supplementary Figures 6, 7). The associations of whole blood selenium with cardiometabolic risk factors did not differ significantly with area, BMI or dyslipidaemia status but differed significantly with sex, age, diabetes status, and hypertension status (Supplementary Figure 6). The associations of urinary selenium with cardiometabolic risk factors did not differ significantly with area, age, sex, or BMI but differed significantly with dyslipidaemia, diabetes, and hypertension status (Supplementary Figure 7).

Sensitivity Analyses

The sensitivity analyses controlling for a history of hypertension, diabetes, or dyslipidaemia (Supplementary Figure 8) and for hs-CRP yielded results similar to those of our main analysis (Supplementary Figure 10). A significant positive association between whole blood selenium and triglyceride levels was also discovered after BMI was excluded from the covariates (Supplementary Figure 9).

Associations Between VVV and IM of Selenium Biomarkers and Cardiometabolic Risk Factors

We observed that SD-WBSe was positively correlated with triglycerides (0.009, 95% CI: 0.000, 0.017) and MSS-sex (0.047, 95% CI: 0.005, 0.088). In addition, ARV-WBSe was positively associated with glucose (0.003, 95% CI: 0.000, 0.005) and MSS-sex (0.034, 95% CI: 0.006, 0.061) in Model 1. However, these associations were non-significant after additional adjustments for the impact of IM in Model 2. Only a marginal association of ARV-WBSe with total cholesterol (−0.003, 95% CI: −0.005, −0.000) was observed in Model 2 (Table 5).

TABLE 5

Table 5. Changes and 95% confidence intervals of cardiometabolic risk factors associated with one-unit increments in VVV and IM of whole blood selenium.

Consistent with results of the main analysis, the IM of whole blood selenium was positively correlated with total cholesterol (0.002, 95% CI: 0.001, 0.004), triglycerides (0.007, 95% CI: 0.004, 0.011), and LDL-C (0.002, 95% CI: 0.000, 0.004) in Model 1, and these associations remained significant after adjustment for the impact of VVV in Model 2. The positive association between the IM of whole blood selenium and glucose (0.002, 95% CI: 0.000, 0.004) was also observed in Model 1, and this association remained significant after adjustment for SD-WBSe, CV-WBSe, and VIM-WBSe in Model 2. Furthermore, we discovered a positive association between the IM of whole blood selenium and MSS-sex (0.022, 95% CI: 0.002, 0.041) in Model 1, and this association was significant after adjustment for CV-WBSe and VIM-WBSe in Model 2 (Table 5).

We did not observe any robust associations of the VVV of urinary selenium with cardiometabolic risk factors, discovering only a marginal association of SD-USe with total cholesterol (−0.004, 95% CI: −0.007, −0.000) after adjusting for the IM of urinary selenium in Model 2. The IM of urinary selenium was positively related to triglycerides (0.014, 95% CI: 0.002, 0.025) and glucose (0.006, 95% CI: 0.000, 0.011) levels after adjustment for SD-USe in Model 2. Although urinary selenium concentration was negatively associated with SBP, DBP, and MAP in the main analysis, the IM of urinary selenium was not associated with SBP, DBP, or MAP. After additional adjustment for VIM-USe in Model 2, the estimated strengths of the associations of the IM of urinary selenium with SBP, DBP, and MAP were −0.029 (95% CI: −0.267, 0.208), −0.050 (95% CI: −0.175, 0.076), and −0.040 (95% CI: −0.188, 0.109), respectively. A negative association between the IM of urinary selenium and waistline (−0.105, 95% CI: −0.195, −0.010) was also discovered, and the association remained significant after adjustment for CV-USe, ARV-Use, and VIM-USe in Model 2. Moreover, a higher IM of urinary selenium was associated with a smaller hipline (−0.085, 95% CI: −0.161, −0.005) after additional adjustment for VIM-USe in Model 2 (Table 6).

TABLE 6

Table 6. Changes and 95% confidence intervals of cardiometabolic risk factors associated with one-unit increments in VVV and IM of urinary selenium.

Discussion

Interpretation of Results

Associations of Selenium Biomarkers Concentrations With Cardiometabolic Risk Factors

In this three-wave repeated-measures study of older adults with high selenium levels, we discovered that whole blood selenium was positively associated with total cholesterol, LDL-C, glucose, and uric acid. Urinary selenium was not associated with any cardiometabolic risk factors except for BP among older adults with high selenium levels.

According to epidemiological studies conducted in China and other countries, the participants in our study had high selenium, with a geometric mean (geometric SD) whole blood selenium level of 134.30 (1.19) μg/L. Although this figure is lower than that recorded in the Enshi area, which is well-known for selenium poisoning in China (66), it exceeds the median level obtained in Wuhan (92.66 μg/L) (12) and mean level observed in Shandong (120 μg/L) (24). In addition, our study population had a lower selenium level than the mean of NHANES respondents in the United States (2015–2016: 191 μg/L) (16) but a higher level than that obtained in the United Kingdom (mean: 86.856 μg/L) (67). The high selenium of our participants may have two explanations: (1) As mentioned, the study settings have selenium-rich soil. (2) Our participants were older adults and they were reported to preferred cereals rich in selenium (68). The low rates of alcohol and cigarette consumption and metabolic and physiological function (e.g., excretion) impairment that occurs during the aging process may also have contributed to the participants' high selenium levels (69, 70).

Consistent with our study, a cross-sectional study of 8,198 Chinese participants from rural areas discovered that serum selenium (mean: 120 μg/L) was positively correlated with triglycerides (0.24, 95% CI: 0.17, 0.31), total cholesterol (0.57, 95% CI: 0.49, 0.65), LDL-C (0.37, 95% CI: 0.32, 0.42), and HDL-C (0.12, 95% CI: 0.08, 0.16) (24). A cross-sectional analysis using NHANES data (2011–2012, N = 2,287) revealed that serum selenium (mean: 192.99 μg/L) was positively associated with LDL-C (0.063, 95% CI: 0.016, 0.110) (5). However, a cross-sectional analysis using NHANES data (2007–2014, N = 7,597) revealed that dietary selenium intake was negatively associated with total cholesterol and HDL-C (71). Furthermore, few longitudinal studies have been conducted, and those that have been conducted have reported controversial results. A 21-year follow-up analysis of the Young Finns Study and an 8-year follow-up analysis of the Olivetti Heart Study discovered positive cross-sectional associations between selenium and lipid concentrations but no such longitudinal relations (72, 73). The inconsistency between the results of our study and those of the two longitudinal studies might be partly attributable to the different population characteristics, substantial disparity in selenium levels, and potential non-linearity of the associations of selenium with lipid profiles. First, the participants of the Young Finns Study and the Olivetti Heart Study were children (aged 3–18 years) and adult males, whose physiological characteristics differ from those of a population of older adults. Second, the mean (SD) selenium concentrations reported in those two studies were 74.3 (14.0) and 77.5 (18.4) μg/L, respectively, which differed from the 134.30 (1.19) μg/L observed in our study. Therefore, large repeated-measures studies of participants with different selenium statuses and demographic characteristics are warranted to further investigate the associations of selenium and lipid profiles and the exposure–response curves.

Whole blood selenium and glucose were positively associated in older adults with high selenium. A systematic review of 15 observational studies enrolling 32,728 participants observed a positive association between selenium concentration and the odds ratio (OR) for diabetes, with a summary OR of 2.03 (95% CI: 1.51, 2.72) (74). Specifically, a cross-sectional analysis of 8,142 middle-aged adults (mean serum selenium: 121.5 μg/L) residing in Linyi, China, revealed that, compared with a low-selenium group (<124.9 μg/L), the ORs for elevated fasting serum glucose for two high-selenium groups (124.9–143.9 and >143.9 μg/L) were 2.31 (95% CI: 1.37, 3.90) and 2.67 (95% CI: 1.59, 4.48) (75). Moreover, a cross-sectional study of NHANES data (1999–2006, N = 41,474) reported that serum selenium concentration (mean: 129 μg/L) was positively associated with plasma glucose (12.454, 95 % CI: 4.122, 20.786) (76). Another NHANES-based (2003–2004, N = 917) cross-sectional analysis of adults aged over 40 years revealed a higher OR of serum selenium–related diabetes (mean: 137.1 mg/L) in women (OR: 5.99) than in men (OR: 2.30) (9). In combination with these results, our study suggests that an individual's sex may affect the selenium–glucose association.

The associations between selenium and other cardiovascular risk factors such as uric acid have not been widely studied. A positive association between whole blood selenium and uric acid among older adults with high selenium levels was observed in the present study. Consistent with our findings, a cross-sectional study of 1,406 Han Chinese adults revealed a positive association between serum selenium and the odds of hyperuricaemia, with an OR of 1.50 (95% CI: 1.01, 2.23) (77). In addition, selenium intake exceeding the recommended amount was positively correlated with the uric acid levels in the first (mean: 53.99 μg/day) and second (mean: 58.93 μg/day) trimesters of pregnancy for 95 Polish women (78). Previous findings suggest that the relationships of selenium biomarkers with uric acid may be confounded by diet and physical activity (79). All three visits in our repeated-measures study were conducted in winter to minimize the effect of seasonally relevant confounding factors, including diet and physical activity which could vary by seasons. High uric acid has been identified as an independent risk factor for long-term CVD events, CVD-related death, and all-cause mortality (80, 81). Furthermore, the positive relationship between selenium and uric acid observed in this study should be interpreted with caution, and additional comprehensive studies are required.

Few studies have explored the relationship between urinary selenium and BP. However, in the present study, we observed that urinary selenium as a continuous variable was negatively associated with SBP, DBP, and MAP among older adults with high selenium levels; however, these negative associations were non-significant when urinary selenium was treated as a categorical variable. A cross-sectional study conducted in Wuhan, China (N = 823), revealed a positive association between increased ORs for hypertension and urinary selenium quartiles (geometric mean: 19.8 μg/g creatinine) by using a single-metal regression model (p trend < 0.05) (26). However, another case–control study conducted in Wuhan (N = 1,004) reported that urinary selenium (geometric mean: 20.47 μg/g creatinine) was not associated with hypertension (82). This disparity can be partly attributed to the different study design and the fact that urinary selenium accounts for the majority of absorbed selenium that is not retained and is easily influenced by diet and physical activity (83).

Associations of VVV and IM of Selenium Biomarkers With Cardiometabolic Risk Factors

In this three-wave repeated-measures study of older adults with high selenium status, the reproducibility of the whole blood selenium measurements was high, and that of urinary selenium was poor. Because the collection of blood samples is invasive, the reproducibility of whole blood selenium measurement has rarely been studied. A longitudinal analysis of urine samples collected from 11 Chinese men at days 0, 1, 2, 3, 4, 30, 60, and 90 suggested poor reproducibility of urinary selenium (ICC = 0.03) (84). The reproducibility of whole blood and urinary selenium measurements may be affected by the kinetics of selenium absorption, distribution, metabolism, and elimination (85). Selenium is absorbed and accumulates in the human body, mainly in the liver, kidneys, and blood. The elimination of selenium in urine is also influenced by diet, lifestyle, and activity (83). Therefore, the reproducibility of whole blood selenium measurement is superior to that of urinary selenium. The reproducibility of the selenium biomarkers observed in the present study suggests that single measurements may not accurately reflect individual exposure to selenium over time, at least among older adults with high selenium, and repeated-measures studies with even more waves should be employed to minimize exposure misclassification.

We explored the relations of the VVV and IM of selenium biomarkers with cardiometabolic risk factors in the present study. Among the older adults with high selenium levels, the VVV of the selenium biomarkers was not significantly associated with cardiometabolic risk factors after adjustments for IM. The VVV of whole blood selenium, including SD-WBSe and ARV-WBSe, was associated with cardiometabolic factors, including glucose, triglycerides, and MSS-sex, but these associations were non-significant after adjustment for IM. VIM-WBSe was also not associated with any cardiometabolic risk factors. Four indicators were calculated to comprehensively evaluate the VVV of the selenium biomarkers. However, in whole blood selenium, these four indicators of VVV were not all associated with any particular cardiometabolic risk factors. Moreover, no stable association between the VVV of urinary selenium and any cardiometabolic risk factor was discovered. Although we did not observe associations between the VVV of the selenium biomarkers and cardiometabolic risk factors among older adults with high selenium levels, future research is suggested. This three-wave repeated-measures study was conducted from November 2016 through January 2018 and captured the long-term trends of the selenium biomarkers. However, repeated-measures studies conducted in different seasons are warranted to characterize seasonal biomarker variability. The VVV of the selenium biomarkers remains an informative measure for selenium biomarkers, and the potential relations of this measure with cardiometabolic risk factors should not be ignored.

We also observed that the IM of whole blood selenium was positively correlated with cardiometabolic factors, including total cholesterol, triglycerides, LDL-C, and MSS-sex, and these associations remained significant after adjustment for VVV. The IM of urinary selenium was also significantly associated with cardiometabolic risk factors after adjustment for VVV. Specifically, the IM of urinary selenium was positively related with triglycerides and glucose and negatively associated with waistline and hipline. In short, the IM of selenium biomarkers can be used in conjunction with the absolute levels for prediction of cardiometabolic risk factors among older adults with high selenium levels. The IM of the selenium biomarkers reflected the centralized levels across the study period. In sum, cardiometabolic risk is associated with exposure to selenium over time among older adults with high selenium levels. Furthermore, a pooled analysis of two prospective population-based cohort studies (i.e., the Health Retirement Study and the English Longitudinal Study of Aging, n = 6,237) revealed that a 1% increase in mean glycosylated hemoglobin A1c was associated with a faster rate of memory function decline (−0.041, 95% CI: −0.071, −0.012) (40). These and our findings suggest that IM in biomarkers and physiological parameters may have an impact on health end points and other studies exploring the relationship between nutrition and disease should consider the impact of IM.

In conclusion, the associations between whole blood and urinary selenium and cardiometabolic risk factors are inconsistent. Except for VVV, concentrations and IM of whole blood selenium were associated with cardiometabolic risk factors. In contrast, only VVV of urinary selenium was significantly associated with cardiometabolic risk factors. These differences could be explained for several reasons. First, selenium enters into the human body mainly by three routes (i.e., digestion, inhalation and skin absorption) (86). Whole body retention studies following oral administration of sodium selenite have indicated that human bodies absorb and accumulate selenium, primarily in the liver, kidneys, and blood (half-life: almost 9 months). The urinary excretion of selenium lasted ~1 week, with a half-life of 8–9 days (87). Hence, whole blood selenium reflected internal burden of selenium and urinary selenium reflected absorbed but unutilized selenium. This selenium metabolism may partly explain the differences in the associations observed in urinary and whole blood selenium. Second, the elimination of urinary selenium is also influenced by diet, lifestyle, and activity (83). Our data also found that the reproducibility of urinary selenium measurement was lower to that of whole blood selenium. Tough we observed negative associations of urinary selenium and SBP, DBP, and MAP, these associations were non-significant when urinary selenium concentration was treated as a categorical variable. Finally, single measurements may not accurately reflect selenium status. IM can capture the centralizing trend of multiple repeated measurements. We did not observe negative association between IM of urinary selenium and BP-related factors. However, both whole blood selenium concentration and IM of whole blood selenium were positively associated with total cholesterol, LDL-C, glucose, and uric acid. Hence, association of whole blood selenium are more stable than that of urinary selenium. Even though the available evidence can partially explain the inconsistent results in the associations observed in urinary and whole blood selenium, experimental animal studies exploring the underlying mechanisms are warrant.

As an essential trace element, selenium is incorporated into selenoproteins that have a wide range of dose-dependent effects. The U-shaped associations between selenium levels and cardiometabolic risk factors must be emphasized (10). A case–control study conducted in Wuhan demonstrated a U-shaped association between plasma selenium (median: 92.66 μg/L) with central obesity and high blood pressure (12). As above mentioned, another cross-sectional analysis using data from NHANES reported a U-shaped association between plasma selenium and the likelihood of diabetes, with the lowest risk noted for a concentration of ~122 μg/L (9). Due to the high blood selenium status (GM:134.30 μg/L), the positive association between whole blood selenium and blood glucose was observed in our study. Additionally, the U-shaped association was observed between selenium status and all-cause and cancer mortality among 13,887 adult participants in the NHANES (1988-1994), with the lowest risk noted for a concentration of ~135 μg/L (88). In order to analyse the two-way effect of selenium on health, especially cardiometabolic risk factors, more epidemiological studies of participants with wide selenium concentration statuses and demographic characteristics were warranted. The exposure response curve of selenium biomarkers and cardiometabolic risk factors among both deficient and excessive selenium levels were needed to be profiled.

Potential Mechanisms

The biological mechanisms explaining selenium's involvement in the pathogenesis of cardiometabolic diseases are not well-understood. High levels of selenium are incorporated into selenoproteins, damaging cardiometabolic health mainly by initiating oxidative stress, promoting insulin resistance, and regulating gluconeogenesis and lipid metabolism. Excess selenium increases the production of reactive oxygen species (ROS) by promoting the overexpression of selenocysteine transfer RNA and reducing selenoproteins synthesis. Increased ROS levels cause insulin resistance in the peripheral tissues by affecting insulin receptor signal transduction, ultimately resulting in hyperinsulinemia and cell glucose desensitization (89). In an experimental study, C57BL/6J mice (n = 6 or 7 per group) on a selenium-supplemented diet (0.1 and 0.4 ppm selenium) developed hyperinsulinemia and reduced insulin sensitivity (90). In addition, excessive ROS production can reduce nitric oxide; accelerate endothelial cell apoptosis; induce upregulation of NF-kB; activate intercellular adhesion molecule-1, monocyte chemotactic protein-1, and vascular cell adhesion molecule-1; and trigger diabetic vascular complications and cardiometabolic disorders (91).

Meanwhile, high dietary selenium intake can increase the expression or activity of key proteins related to gluconeogenesis, glycolysis, and lipogenesis by upregulating the selenoproteins glutathione peroxidase family. First, high glutathione peroxidase-1 production due to high selenium intake results in the upregulation of phosphoenolpyruvate carboxykinase and upregulation of fatty acid synthesis. In an experimental study, prolonged high dietary intake of selenium was demonstrated to induce gestational diabetes in rats and hyperinsulinemia in pigs (92). Second, high selenium upregulates protein tyrosine phosphatase 1B (PTP1B), a key enzyme in triggering fatty acid synthesis and in reverse regulation of insulin signaling, eventually inducing lipid disorders and insulin resistance. In one animal study, PTP1B expression were elevated in rats with fructose-rich diets, resulting in the induction of fatty acid synthesis and an increase in liver triglycerides (93). Further studies are warranted to clarify the mechanisms underlying the complex associations between selenium and cardiometabolic risk.

Limitations and Strengths

This study had several limitations. First, we selected whole blood selenium and urinary selenium as the selenium biomarkers to investigate, but serum or plasma selenium may have been more favorable alternatives. Nevertheless, in our study, whole blood selenium and urinary selenium were preferable because of the consistency of concentrations and determinants between the whole blood and serum selenium; in addition, whole blood selenium is reportedly a suitable indicator of medium- to long-term selenium status (94). Moreover, measuring both whole blood and urinary selenium provides an accurate picture of the function and excretion characteristics of selenium (83). Second, information regarding the participants' alcohol consumption, cigarette smoking habits, and disease history was collected using questionnaires. Consequently, unmeasured confounding factors and recall bias likely limit the generalisability of our findings. Third, the modest number of repeat measurements and their short time intervals might have hampered our ability to characterize selenium variability. Despite the limited data, we provide a fresh perspective on the relationships between selenium and cardiometabolic risk factors. Finally, the participants in our study were older adults (age ≥50 years) with high selenium levels; thus, the generalisability of the results to more general populations may be limited.

Despite these limitations, several strengths should also be considered. The three-wave repeated-measures design enabled repeated data collection and inference of causal relationships. Additionally, all three measurements were conducted in winter, minimizing the effect of seasonally relevant confounding factors. To our knowledge, this is the first study to comprehensively evaluate the associations of selenium concentration, VVV, IM with cardiometabolic risk factors. Finally, our findings are derived from observations of older adults with high blood selenium in the Beijing region, elucidating the health effects of high selenium and providing science-based evidence for nutritional guidelines in the region.

Conclusion

We discovered that selenium levels and their IMs were significantly associated with several cardiometabolic factors, namely total cholesterol, LDL-C, and glucose, in older adults with high selenium in the Beijing area. This indicates that selenium affects cardiometabolic risk. However, we do not observe any robust associations between the VVV of the selenium biomarkers and cardiometabolic risk factors after adjustment for IM. The findings suggest that older adults with high selenium should not take dietary selenium supplements to prevent cardiometabolic risk. In the future, longer repeated-measures studies of the general population are warranted to minimize selenium exposure misclassification, explore the associations of the VVV and IM of selenium biomarkers with cardiometabolic risk factors, and determine the relevant underlying mechanisms.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

AL: conceptualization, methodology, data-analysis and interpretation, writing-original draft, writing-review and editing, and funding acquisition. QZ: data-cleaning and interpretation, writing-original draft, and writing-review and editing. YM and JZ: investigation, data cleaning, and writing-review and editing. MZ, JX, and XG: investigation and writing-review and editing. QX: funding acquisition, writing-review and editing, and supervision. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the China Medical Board (Grant No. 15-230), the Fundamental Research Funds for the Central Universities (Grant No. 3332019147), Peking Union Medical College Graduate Innovation Fund (No. 2019-1004-02), the China Prospective cohort study of Air pollution and health effects in Typical areas (C-PAT; Grant No. MEE-EH-20190802), and the Chinese Academy of Medical Science Innovation Fund for Medical Sciences (Grant No. 2017-I2M-1-009).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We are grateful for the kind cooperation of all participants in this study. Also, we would like to thank all the investigators involved in our study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.838613/full#supplementary-material

References

1. Roth GA, Mensah GA, Johnson CO, Addolorato G, Ammirati E, Baddour LM, et al. Global burden of cardiovascular diseases and risk factors, 1990-2019: update from the GBD 2019 study. J Am Coll Cardiol. (2020) 76:2982–3021. doi: 10.1016/j.jacc.2020.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Singh GM, Danaei G, Farzadfar F, Stevens GA, Woodward M, Wormser D, et al. The age-specific quantitative effects of metabolic risk factors on cardiovascular diseases and diabetes: a pooled analysis. PLoS ONE. (2013) 8:e65174. doi: 10.1371/journal.pone.0065174

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Whitlock G, Lewington S, Sherliker P, Clarke R, Emberson J, Halsey J, et al. Body-mass index and cause-specific mortality in 900,000 adults: collaborative analyses of 57 prospective studies. Lancet. (2009) 373:1083–96. doi: 10.1016/S0140-6736(09)60318-4

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wang N, Tan HY, Li S, Xu Y, Guo W, Feng Y. Supplementation of micronutrient selenium in metabolic diseases: its role as an antioxidant. Oxid Med Cell Longev. (2017) 2017:7478523. doi: 10.1155/2017/7478523

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Christensen K, Werner M, Malecki K. Serum selenium and lipid levels: associations observed in the National Health and Nutrition Examination Survey (NHANES) 2011-2012. Environ Res. (2015) 140:76–84. doi: 10.1016/j.envres.2015.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Vinceti M, Filippini T, Rothman KJ. Selenium exposure and the risk of type 2 diabetes: a systematic review and meta-analysis. Eur J Epidemiol. (2018) 33:789–810. doi: 10.1007/s10654-018-0422-8

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Retondario A, Fernandes R, Rockenbach G, Alves MA, Bricarello LP, Trindade E, et al. Selenium intake and metabolic syndrome: a systematic review. Clin Nutr. (2019) 38:603–14. doi: 10.1016/j.clnu.2018.02.021

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang Y, Rijntjes E, Wu Q, Lv H, Gao C, Shi B, et al. Selenium deficiency is linearly associated with hypoglycemia in healthy adults. Redox Biol. (2020) 37:101709. doi: 10.1016/j.redox.2020.101709

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Laclaustra M, Navas-Acien A, Stranges S, Ordovas JM, Guallar E. Serum selenium concentrations and diabetes in U.S. adults: National Health and Nutrition Examination Survey (NHANES) 2003-2004. Environ Health Perspect. (2009) 117:1409–13. doi: 10.1289/ehp.0900704

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Rayman MP. Selenium and human health. Lancet. (2012) 379:1256–68. doi: 10.1016/S0140-6736(11)61452-9

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Iyengar GV. Reference values for the concentrations of As, Cd, Co, Cr, Cu, Fe, I, Hg, Mn, Mo, Ni, Pb, Se, and Zn in selected human tissues and body fluids. Biol Trace Elem Res. (1987) 12:263–95. doi: 10.1007/BF02796686

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zhou L, Luo C, Yin J, Zhu Y, Li P, Chen S, et al. Diverse associations of plasma selenium concentrations and SELENOP gene polymorphism with metabolic syndrome and its components. Oxid Med Cell Longev. (2020) 2020:5343014. doi: 10.1155/2020/5343014

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Lu CW, Chang HH, Yang KC, Chiang CH, Yao CA, Huang KC. Gender differences with dose(-)response relationship between serum selenium levels and metabolic syndrome-a case-control study. Nutrients. (2019) 11:20477. doi: 10.3390/nu11020477

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Fang C, Wu W, Gu X, Dai S, Zhou Q, Deng H, et al. Association of serum copper, zinc and selenium levels with risk of metabolic syndrome: a nested case-control study of middle-aged and older Chinese adults. J Trace Elem Med Biol. (2019) 52:209–15. doi: 10.1016/j.jtemb.2018.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhong Q, Lin R, Nong Q. Adiposity and serum selenium in U.S. adults. Nutrients. (2018) 10:60727. doi: 10.3390/nu10060727

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Centers for Disease Control and Prevention. National Report on Human Exposure to Environmental Chemicals. Centers for Disease Control and Prevention (2019).

Google Scholar

17. Guo YZ-F., Ruan Q-H, Hu S-Y, Zhou E-N, Feng H. Content and distribution of selenium in soil of Beijing. Plain Geosci. (2012) 26:859–64. Available online at: http://www.geoscience.net.cn/CN/Y2012/V26/I5/859

Google Scholar

18. Liu S, Zhang M, Yang L, Li Y, Wang L, Huang Z, et al. Prevalence and patterns of tobacco smoking among Chinese adult men and women: findings of the 2010 national smoking survey. J Epidemiol Community Health. (2017) 71:154–61. doi: 10.1136/jech-2016-207805

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhenping Z, Limin W, Mei Z, Xiao Z, Zhengjing H, Chun L, et al. Geographic distribution of alcohol use among chinese adults — China, 2015. China CDC Weekly. (2020) 2:98–103. doi: 10.46234/ccdcw2020.029

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Sánchez C, López-Jurado M, Aranda P, Llopis J. Plasma levels of copper, manganese and selenium in an adult population in southern Spain: influence of age, obesity and lifestyle factors. Sci Total Environ. (2010) 408:1014–20. doi: 10.1016/j.scitotenv.2009.11.041

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Gerardo B, Cabral Pinto M, Nogueira J, Pinto P, Almeida A, Pinto E, et al. Associations between trace elements and cognitive decline: an exploratory 5-year follow-up study of an elderly cohort. Int J Environ Res Public Health. (2020) 17:176051. doi: 10.3390/ijerph17176051

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Bizerea-Moga TO, Pitulice L, Bizerea-Spiridon O, Moga TV. Evaluation of serum selenium status by age and gender: a retrospective observational cohort study in Western Romania. Nutrients. (2021) 13:51497. doi: 10.3390/nu13051497

PubMed Abstract | CrossRef Full Text | Google Scholar

23. González-Estecha M, Palazón-Bru I, Bodas-Pinedo A, Trasobares E, Palazón-Bru A, Fuentes M, et al. Relationship between serum selenium, sociodemographic variables, other trace elements and lipid profile in an adult Spanish population. J Trace Elem Med Biol. (2017) 43:93–105. doi: 10.1016/j.jtemb.2016.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ju W, Ji M, Li X, Li Z, Wu G, Fu X, et al. Relationship between higher serum selenium level and adverse blood lipid profile. Clin Nutr. (2018) 37:1512–7. doi: 10.1016/j.clnu.2017.08.025

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Fabbri E, Zoli M, Gonzalez-Freire M, Salive ME, Studenski SA, Ferrucci L. Aging and multimorbidity: new tasks, priorities, and frontiers for integrated gerontological and clinical research. J Am Med Dir Assoc. (2015) 16:640–7. doi: 10.1016/j.jamda.2015.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Wu W, Jiang S, Zhao Q, Zhang K, Wei X, Zhou T, et al. Environmental exposure to metals and the risk of hypertension: a cross-sectional study in China. Environ Pollut. (2018) 233:670–8. doi: 10.1016/j.envpol.2017.10.111

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Laclaustra M, Navas-Acien A, Stranges S, Ordovas JM, Guallar E. Serum selenium concentrations and hypertension in the US population. Circ Cardiovasc Qual Outcomes. (2009) 2:369–76. doi: 10.1161/CIRCOUTCOMES.108.831552

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Wang YX, Feng W, Zeng Q, Sun Y, Wang P, You L, et al. Variability of metal levels in spot, first morning, and 24-hour urine samples over a 3-month period in healthy adult Chinese Men. Environ Health Perspect. (2016) 124:468–76. doi: 10.1289/ehp.1409551

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Nwabuo CC, Yano Y, Moreira HT, Appiah D, Vasconcellos HD, Aghaji QN, et al. Association between visit-to-visit blood pressure variability in early adulthood and myocardial structure and function in later life. J Am Med Assoc Cardiol. (2020) 5:795–801. doi: 10.1001/jamacardio.2020.0799

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Wan EYF, Yu EYT, Chin WY, Barrett JK, Mok AHY, Lau CST, et al. Greater variability in lipid measurements associated with cardiovascular disease and mortality: a 10-year diabetes cohort study. Diabetes Obes Metab. (2020) 22:1777–88. doi: 10.1111/dom.14093

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Echouffo-Tcheugui JB, Zhao S, Brock G, Matsouaka RA, Kline D, Joseph JJ. Visit-to-visit glycemic variability and risks of cardiovascular events and all-cause mortality: the ALLHAT study. Diabetes Care. (2019) 42:486–93. doi: 10.2337/dc18-1430

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zeng H. Selenium as an essential micronutrient: roles in cell cycle and apoptosis. Molecules. (2009) 14:1263–78. doi: 10.3390/molecules14031263

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Li A, Pei L, Zhao M, Xu J, Mei Y, Li R, et al. Investigating potential associations between O3 exposure and lipid profiles: a longitudinal study of older adults in Beijing. Environ Int. (2019) 133:105135. doi: 10.1016/j.envint.2019.105135

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Barr DB, Wilder LC, Caudill SP, Gonzalez AJ, Needham LL, Pirkle JL. Urinary creatinine concentrations in the US population: implications for urinary biologic monitoring measurements. Environ Health Perspect. (2005) 113:192–200. doi: 10.1289/ehp.7337

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Hoet P, Deumer G, Bernard A, Lison D, Haufroid V. Urinary trace element concentrations in environmental settings: is there a value for systematic creatinine adjustment or do we introduce a bias? J Expo Sci Environ Epidemiol. (2016) 26:296–302. doi: 10.1038/jes.2015.23

PubMed Abstract | CrossRef Full Text | Google Scholar

36. O'Brien KM, Upson K, Buckley JP. Lipid and creatinine adjustment to evaluate health effects of environmental exposures. Curr Environ Health Rep. (2017) 4:44–50. doi: 10.1007/s40572-017-0122-7

PubMed Abstract | CrossRef Full Text | Google Scholar

37. O'Brien KM, Upson K, Cook NR, Weinberg CR. Environmental chemicals in urine and blood: improving methods for creatinine and lipid adjustment. Environ Health Perspect. (2016) 124:220–7. doi: 10.1289/ehp.1509693

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Bulka CM, Mabila SL, Lash JP, Turyk ME, Argos M. Arsenic and obesity: a comparison of urine dilution adjustment methods. Environ Health Perspect. (2017) 125:087020. doi: 10.1289/EHP1202

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Lee I, Park YJ, Kim MJ, Kim S, Choi S, Park J, et al. Associations of urinary concentrations of phthalate metabolites, bisphenol A, and parabens with obesity and diabetes mellitus in a Korean adult population: Korean National Environmental Health Survey (KoNEHS) 2015-2017. Environ Int. (2021) 146:106227. doi: 10.1016/j.envint.2020.106227

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Yu ZB, Zhu Y, Li D, Wu MY, Tang ML, Wang JB, et al. Association between visit-to-visit variability of HbA(1c) and cognitive decline: a pooled analysis of two prospective population-based cohorts. Diabetologia. (2020) 63:85–94. doi: 10.1007/s00125-019-04986-8

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Yang CD, Shen Y, Ding FH, Yang ZK, Hu J, Shen WF, et al. Visit-to-visit fasting plasma glucose variability is associated with left ventricular adverse remodeling in diabetic patients with STEMI. Cardiovasc Diabetol. (2020) 19:131. doi: 10.1186/s12933-020-01112-6

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Asayama K, Wei FF, Liu YP, Hara A, Gu YM, Schutte R, et al. Does blood pressure variability contribute to risk stratification? Methodological issues and a review of outcome studies based on home blood pressure. Hypertens Res. (2015) 38:97–101. doi: 10.1038/hr.2014.153

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Rothwell PM. Limitations of the usual blood-pressure hypothesis and importance of variability, instability, and episodic hypertension. Lancet. (2010) 375:938–48. doi: 10.1016/S0140-6736(10)60309-1

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Melgarejo JD, Yang WY, Thijs L, Li Y, Asayama K, Hansen TW, et al. Association of fatal and nonfatal cardiovascular outcomes with 24-hour mean arterial pressure. Hypertension. (2021) 77:39–48. doi: 10.1161/HYPERTENSIONAHA.120.14929

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Nishida C, Ko GT, Kumanyika S. Body fat distribution and noncommunicable diseases in populations: overview of the 2008 WHO Expert Consultation on Waist Circumference and Waist-Hip Ratio. Eur J Clin Nutr. (2010) 64:2–5. doi: 10.1038/ejcn.2009.139

PubMed Abstract | CrossRef Full Text | Google Scholar

46. DeBoer MD, Gurka MJ. Clinical utility of metabolic syndrome severity scores: considerations for practitioners. Diabetes Metab Syndr Obes. (2017) 10:65–72. doi: 10.2147/DMSO.S101624

PubMed Abstract | CrossRef Full Text | Google Scholar

47. DeBoer MD, Gurka MJ, Woo JG, Morrison JA. Severity of metabolic syndrome as a predictor of cardiovascular disease between childhood and adulthood: the princeton lipid research cohort study. J Am Coll Cardiol. (2015) 66:755–7. doi: 10.1016/j.jacc.2015.05.061

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Magnussen CG, Cheriyan S, Sabin MA, Juonala M, Koskinen J, Thomson R, et al. Continuous and dichotomous metabolic syndrome definitions in youth predict adult type 2 diabetes and carotid artery intima media thickness: the cardiovascular risk in young Finns study. J Pediatr. (2016) 171:97–103.e1-3. doi: 10.1016/j.jpeds.2015.10.093

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Gurka MJ, Lilly CL, Oliver MN, DeBoer MD. An examination of sex and racial/ethnic differences in the metabolic syndrome among adults: a confirmatory factor analysis and a resulting continuous severity score. Metabolism. (2014) 63:218–25. doi: 10.1016/j.metabol.2013.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Kraus WE, Bhapkar M, Huffman KM, Pieper CF, Krupa Das S, Redman LM, et al. 2 years of calorie restriction and cardiometabolic risk (CALERIE): exploratory outcomes of a multicentre, phase 2, randomised controlled trial. Lancet Diabetes Endocrinol. (2019) 7:673–83. doi: 10.1016/S2213-8587(19)30151-2

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Joint Committee for Guideline Revision. 2016 Chinese guidelines for the management of dyslipidemia in adults. J Geriatr Cardiol. (2018) 15:1–29. doi: 10.11909/j.issn.1671-5411.2018.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Chobanian AV, Bakris GL, Black HR, Cushman WC, Green LA, Izzo JL Jr, et al. Seventh report of the joint national committee on prevention, detection, evaluation, and treatment of high blood pressure. Hypertension. (2003) 42:1206–52. doi: 10.1161/01.HYP.0000107251.49515.c2

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Association AD. Standards of medical care in diabetes-−2014. Diabetes Care. (2014) 37(Suppl.1):S14–80. doi: 10.2337/dc14-S014

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Gu YQ, Xia Y, Zhang SM, Bao X, Wu HM, Bian SS, et al. Method of dietary nutritional status assessment and its application in cohort study of nutritional epidemiology. Zhonghua Liu Xing Bing Xue Za Zhi. (2020) 41:1145–50. doi: 10.3760/cma.j.cn112338-20200110-00027

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Wang Y. China Food Composition Tables. Beijing: Peking University Medical Press (2018).

Google Scholar

56. Alajbeg IZ, Lapić I, Rogić D, Vuletić L, Andabak Rogulj A, Illeš D, et al. Within-subject reliability and between-subject variability of oxidative stress markers in saliva of healthy subjects: a longitudinal pilot study. Dis Markers. (2017) 2017:2697464. doi: 10.1155/2017/2697464

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Cicchetti DV. Guidelines, criteria, and rules of thumb for evaluating normed and standardized assessment instruments in psychology. Psychol Assess. (1994) 6:284–90. doi: 10.1037/1040-3590.6.4.284

CrossRef Full Text | Google Scholar

58. Textor J, Hardt J, Knüppel S. DAGitty: a graphical tool for analyzing causal diagrams. Epidemiology. (2011) 22:745. doi: 10.1097/EDE.0b013e318225c2be

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Textor J, van der Zander B, Gilthorpe MS, Liskiewicz M, Ellison GT. Robust causal inference using directed acyclic graphs: the R package 'dagitty'. Int J Epidemiol. (2016) 45:1887–94. doi: 10.1093/ije/dyw341

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Jafari M, Ansari-Pour N. Why, when and how to adjust your P-values? Cell J. (2019) 20:604–7. doi: 10.22074/cellj.2019.5992

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Li A, Mei Y, Zhao M, Xu J, Seery S, Li R, et al. The effect of ambient ozone on glucose-homoeostasis: a prospective study of non-diabetic older adults in Beijing. Sci Total Environ. (2021) 761:143308. doi: 10.1016/j.scitotenv.2020.143308

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Bulka CM, Scannell Bryan M, Persky VW, Daviglus ML, Durazo-Arvizu RA, Parvez F, et al. Changes in blood pressure associated with lead, manganese, and selenium in a Bangladeshi cohort. Environ Pollut. (2019) 248:28–35. doi: 10.1016/j.envpol.2019.01.129

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Weihrauch-Blüher S, Schwarz P, Klusmann JH. Childhood obesity: increased risk for cardiometabolic disease and cancer in adulthood. Metabolism. (2019) 92:147–52. doi: 10.1016/j.metabol.2018.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Abu-Farha M, Behbehani K, Elkum N. Comprehensive analysis of circulating adipokines and hsCRP association with cardiovascular disease risk factors and metabolic syndrome in Arabs. Cardiovasc Diabetol. (2014) 13:76. doi: 10.1186/1475-2840-13-76

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Mazidi M, Toth PP, Banach M. C-reactive protein is associated with prevalence of the metabolic syndrome, hypertension, and diabetes mellitus in US adults. Angiology. (2018) 69:438–42. doi: 10.1177/0003319717729288

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Qin HB, Zhu JM, Liang L, Wang MS, Su H. The bioavailability of selenium and risk assessment for human selenium poisoning in high-Se areas, China. Environ Int. (2013) 52:66–74. doi: 10.1016/j.envint.2012.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Stranges S, Laclaustra M, Ji C, Cappuccio FP, Navas-Acien A, Ordovas JM, et al. Higher selenium status is associated with adverse blood lipid profile in British adults. J Nutr. (2010) 140:81–7. doi: 10.3945/jn.109.111252

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Jian X. Analysis of the dietary factors on sarcopenia in elderly in Beijing. J Hygiene Res. (2016) 45:388–93.

PubMed Abstract | Google Scholar

69. He H, Pan L, Cui Z, Sun J, Yu C, Cao Y, et al. Smoking prevalence, patterns, and cessation among adults in Hebei Province, Central China: Implications From China National Health Survey (CNHS). Front Public Health. (2020) 8:177. doi: 10.3389/fpubh.2020.00177

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Parascandola M, Xiao L. Tobacco and the lung cancer epidemic in China. Transl Lung Cancer Res. (2019) 8:S21–30. doi: 10.21037/tlcr.2019.03.12

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Ma X, Jiang S, Yan S, Li M, Wang C, Pan Y, et al. Association between copper, zinc, iron, and selenium intakes and TC/HDL-C ratio in US adults. Biol Trace Elem Res. (2020) 197:43–51. doi: 10.1007/s12011-019-01979-x

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Stranges S, Tabák AG, Guallar E, Rayman MP, Akbaraly TN, Laclaustra M, et al. Selenium status and blood lipids: the cardiovascular risk in Young Finns study. J Intern Med. (2011) 270:469–77. doi: 10.1111/j.1365-2796.2011.02398.x

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Stranges S, Galletti F, Farinaro E, D'Elia L, Russo O, Iacone R, et al. Associations of selenium status with cardiometabolic risk factors: an 8-year follow-up analysis of the Olivetti Heart study. Atherosclerosis. (2011) 217:274–8. doi: 10.1016/j.atherosclerosis.2011.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Kohler LN, Foote J, Kelley CP, Florea A, Shelly C, Chow HS, et al. Selenium and type 2 diabetes: systematic review. Nutrients. (2018) 10:121924. doi: 10.3390/nu10121924

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Li Z, Li X, Ju W, Wu G, Yang X, Fu X, et al. High serum selenium levels are associated with impaired fasting glucose and elevated fasting serum glucose in Linyi, China. J Trace Elem Med Biol. (2018) 45:64–9. doi: 10.1016/j.jtemb.2017.09.023

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Lin J, Shen T. Association of dietary and serum selenium concentrations with glucose level and risk of diabetes mellitus: a cross sectional study of national health and nutrition examination survey, 1999-2006. J Trace Elem Med Biol. (2021) 63:126660. doi: 10.1016/j.jtemb.2020.126660

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Wang T, Lv Z, Wen Y, Zou X, Zhou G, Cheng J, et al. Associations of plasma multiple metals with risk of hyperuricemia: a cross-sectional study in a mid-aged and older population of China. Chemosphere. (2021) 287:132305. doi: 10.1016/j.chemosphere.2021.132305

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Pieczyńska J, Płaczkowska S, Sozański R, Orywal K, Mroczko B, Grajeta H. Is maternal dietary selenium intake related to antioxidant status and the occurrence of pregnancy complications? J Trace Elem Med Biol. (2019) 54:110–7. doi: 10.1016/j.jtemb.2019.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Yusuf S, Hawken S, Ounpuu S, Dans T, Avezum A, Lanas F, et al. Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries (the INTERHEART study): case-control study. Lancet. (2004) 364:937–52. doi: 10.1016/S0140-6736(04)17018-9

PubMed Abstract | CrossRef Full Text | Google Scholar

80. Puddu PE, Bilancio G, Terradura Vagnarelli O, Lombardi C, Mancini M, Zanchetti A, et al. Serum uric acid and eGFR_CKDEPI differently predict long-term cardiovascular events and all causes of deaths in a residential cohort. Int J Cardiol. (2014) 171:361–7. doi: 10.1016/j.ijcard.2013.12.029

PubMed Abstract | CrossRef Full Text | Google Scholar

81. de Koning L, Merchant AT, Pogue J, Anand SS. Waist circumference and waist-to-hip ratio as predictors of cardiovascular events: meta-regression analysis of prospective studies. Eur Heart J. (2007) 28:850–6. doi: 10.1093/eurheartj/ehm026

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Wu W, Jiang S, Zhao Q, Zhang K, Wei X, Zhou T, et al. Associations of environmental exposure to metals with the risk of hypertension in China. Sci Total Environ. (2018) 622–3:184–91. doi: 10.1016/j.scitotenv.2017.11.343

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Combs GF Jr. Biomarkers of selenium status. Nutrients. (2015) 7:2209–36. doi: 10.3390/nu7042209

PubMed Abstract | CrossRef Full Text | Google Scholar

84. Chen HG, Chen YJ, Chen C, Tu ZZ, Lu Q, Wu P, et al. Reproducibility of essential elements chromium, manganese, iron, zinc and selenium in spot samples, first-morning voids and 24-h collections from healthy adult men. Br J Nutr. (2019) 122:343–51. doi: 10.1017/S0007114519001193

PubMed Abstract | CrossRef Full Text | Google Scholar

85. Aylward LL, Hays SM, Smolders R, Koch HM, Cocker J, Jones K, et al. Sources of variability in biomarker concentrations. J Toxicol Environ Health B Crit Rev. (2014) 17:45–61. doi: 10.1080/10937404.2013.864250

PubMed Abstract | CrossRef Full Text | Google Scholar

86. Mehdi Y, Hornick JL, Istasse L, Dufrasne I. Selenium in the environment, metabolism and involvement in body functions. Molecules. (2013) 18:3292–311. doi: 10.3390/molecules18033292

PubMed Abstract | CrossRef Full Text | Google Scholar

87. Griffiths NM, Stewart RD, Robinson MF. The metabolism of [75Se]selenomethionine in four women. Br J Nutr. (1976) 35:373–82. doi: 10.1079/BJN19760043

PubMed Abstract | CrossRef Full Text | Google Scholar

88. Bleys J, Navas-Acien A, Guallar E. Serum selenium levels and all-cause, cancer, and cardiovascular mortality among US adults. Arch Intern Med. (2008) 168:404–10. doi: 10.1001/archinternmed.2007.74

PubMed Abstract | CrossRef Full Text | Google Scholar

89. Hurrle S, Hsu WH. The etiology of oxidative stress in insulin resistance. Biomed J. (2017) 40:257–62. doi: 10.1016/j.bj.2017.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

90. Labunskyy VM, Lee BC, Handy DE, Loscalzo J, Hatfield DL, Gladyshev VN. Both maximal expression of selenoproteins and selenoprotein deficiency can promote development of type 2 diabetes-like phenotype in mice. Antioxid Redox Signal. (2011) 14:2327–36. doi: 10.1089/ars.2010.3526

PubMed Abstract | CrossRef Full Text | Google Scholar

91. Yamagishi S, Nakamura K, Matsui T. Role of oxidative stress in the development of vascular injury and its therapeutic intervention by nifedipine. Curr Med Chem. (2008) 15:172–7. doi: 10.2174/092986708783330557

PubMed Abstract | CrossRef Full Text | Google Scholar

92. Zhao Z, Barcus M, Kim J, Lum KL, Mills C, Lei XG. High dietary selenium intake alters lipid metabolism and protein synthesis in liver and muscle of pigs. J Nutr. (2016) 146:1625–33. doi: 10.3945/jn.116.229955

PubMed Abstract | CrossRef Full Text | Google Scholar

93. Anurag P, Anuradha CV. Metformin improves lipid metabolism and attenuates lipid peroxidation in high fructose-fed rats. Diabetes Obes Metab. (2002) 4:36–42. doi: 10.1046/j.1463-1326.2002.00178.x

PubMed Abstract | CrossRef Full Text | Google Scholar

94. Thomson CD. Assessment of requirements for selenium and adequacy of selenium status: a review. Eur J Clin Nutr. (2004) 58:391–402. doi: 10.1038/sj.ejcn.1601800

PubMed Abstract | CrossRef Full Text | Google Scholar

Glossary

ARV, average real variability; ARV-Use, average real variability of urinary selenium; ARV-WBSe, average real variability of whole blood selenium; BMI, body mass index; CI, confidence interval; CV, coefficient of variation; CVD, cardiovascular disease; CV-Use, coefficient of variation of urinary selenium; CV-WBSe, coefficient of variation of whole blood selenium; DBP, diastolic blood pressure; EDF, estimated degrees of freedom; FAS, fatty acid synthesis; FDR, false discovery rate; FFQ, food frequency questionnaire; GAMMs, generalized additive mixed models; HDL-C, high-density lipoprotein cholesterol; hs-CRP, high-sensitivity C-reactive protein; ICCs, intraclass correlation coefficients; ICP-MS, inductively coupled plasma mass spectrometer; LDL-C, low-density lipoprotein cholesterol; LOD, limit of detection; MAP, mean arterial pressure; MSS, metabolic syndrome score; MSS-sex, sex-specific metabolic syndrome score; NHANES, the National Health and Nutrition Examination Survey; OR, odds ratio; PTP1B, protein tyrosine phosphatase 1B; ROS, reactive oxygen species; SBP, systolic blood pressure; SD, standard deviation; SD-USe, standard deviation of urinary selenium; SD-WBSe, standard deviation of whole blood selenium; Use, urinary selenium; VIM, variability independent of the mean; VIM-USe, variability independent of the mean of urinary selenium; VIM-WBSe, variability independent of the mean of whole blood selenium; WBSe, whole blood selenium.