Abstract

Background. The systemic immune-inflammation index (SII) is used as an indicator of prognosis for a wide range of diseases. Thyroid function has been found to be strongly associated with inflammation. The purpose of this investigation was to analyze the correlation between SII and various thyroid functions. Methods. This study utilized data from the National Health and Nutrition Examination Survey (NHANES) 2007–2012. The association between SII and thyroid function was analyzed using weighted univariate and multivariate linear regression analyses. Subgroup analyses, interaction tests, and weighted restricted cubic spline (RCS) regression analyses were also employed to test this correlation. Results. Of the 6,875 participants (age ≥ 20 years), the mean age was 46.87 ± 0.40 years. The adjusted model showed that lnSII was negatively correlated with FT3 (β = −0.0559, 95% CI −0.1060 to −0.0059,) and FT3/FT4 (β = −0.0920, 95% CI −0.1667 to −0.0173,). There was a positive correlation between lnSII and TT4 (β = 0.1499, 95% CI 0.0722–0.2276,). In subgroup analyses, lnSII still independently affected a wide range of thyroid functions. Weighted RCS analysis showed a nonlinear relationship between FT3 and lnSII. Conclusion. Close relationships exist between SII and a variety of thyroid functions. SII can be used as an indicator to predict thyroid dysfunction. Control of inflammatory activity may be a protective measure against thyroid dysfunction. More large-scale prospective studies are necessary to further explore the correlation between SII and thyroid function and the role of obesity in this.

1. Introduction

Triiodothyronine (T3) and tetraiodothyronine (T4) are the two major forms of thyroid hormones that are released by the thyroid gland and are regulated by the hypothalamic–pituitary–thyroid axis [1]. The synthesis and secretion of thyroid hormones are stimulated by thyroid stimulating hormone (TSH), which is synthesized by the adenohypophysis. T3 and T4 secreted into the circulation are tightly bound to plasma proteins, a small proportion of these that are biologically active are known as free triiodothyronine (FT3) and free thyroxine (FT4), respectively [2]. Thyroglobulin antibodies (TgAb) and thyroid peroxidase antibodies (TPOAb) are related to autoimmune thyroid disease and are commonly used to identify patients at high risk of progressing to hypothyroidism requiring treatment [3].

Thyroid hormones are crucial in various physiological processes, including growth, normal development, neural differentiation, and metabolic regulation. It can regulate energy expenditure and thermogenesis, the composition of lipids, and ensure the maintenance of normal reproductive function [4, 5]. Fluctuations in thyroid hormone have been correlated with a variety of diseases, an excess of thyroid hormone may lead to atrial fibrillation and heart failure [6, 7]. Insufficient thyroid hormone is associated with a high risk of hypertension, dyslipidemia, and type 2 diabetes [8, 9]. Due to the significant negative impact of adverse clinical outcomes from thyroid dysfunction, a high degree of clinical attention should be paid to changes in thyroid function.

Previous studies have found that abnormal thyroid function is closely related to the inflammatory response. Large population-based studies have revealed a potential association between systemic inflammation and the formation and progress of thyroid nodules, and inflammatory markers are strongly correlated with the size of thyroid nodules in all genders [10]. COVID-19-mediated overproduction of inflammatory cytokines like IL-6, IL-1β, and interferon-γ increased the incidence of subacute thyroiditis, Graves’ disease, and Hashimoto’s thyroiditis, which affect thyroid function [11]. Animal studies have found that inflammatory vesicle activation in obese mice upregulates deiodinase activity by downregulating ubiquitin-mediated protein degradation. This contributes to thyroid hormone resistance at the pituitary level and has an effect on thyroid hormones [12, 13].

The systemic immune-inflammation index (SII) is widely recognized as a stable and comprehensive indicator of the body’s local immune response and systemic inflammatory state [14, 15]. SII is a composite index based on peripheral blood lymphocyte, neutrophil, and platelet counts (SII = platelet counts × neutrophil/lymphocyte counts). Several previous studies have shown that SII can predict the prognosis of tumor patients by reflecting the balance between the body’s inflammatory and immune responses and has high prognostic value for patients with prostate cancer [16], hepatocellular carcinoma [14], cervical cancer [17], and biliary tract cancer [18]. SII provided more clinical information than previous inflammatory factors. Patients with high levels of SII were associated with thrombocytosis, neutrophilia, or lymphocytopenia [19]. It is now widely recognized that SII accurately assesses the inflammatory balance of the body [20], which reflects the degree of systemic inflammatory activity, with a higher SII representing a more pronounced inflammatory response. A higher SII was found to be correlated with an increased risk of type 2 diabetes mellitus combined with depression [21], urinary albumin excretion [22], osteoporosis in postmenopausal women [23], peripheral arterial disease [24], ulcerative colitis [25], and diabetic nephropathy [20]. However, there are no studies on SII and thyroid function at this time.

To fill this gap, the objective of this study was to explore the correlation between SII and thyroid function using a nationally representative sample of adults (age ≥ 20 years) in the United States.

2. Materials and Methods

2.1. The Subjects in the Study

Subject data for this study were obtained from the National Health and Nutrition Examination Surveys (NHANES). NHANES is a nationwide, population-based, and cross-sectional survey dedicated to collecting information on the nutritional status and health risk factors of U.S. nationals. NHANES incorporates nationally representative subjects through a well-designed, complex, stratified, and multistage sampling design. Informed permission was obtained in writing from all participants before they were enrolled in the investigation. The study design and written informed consent were examined and confirmed by the Research Ethics Review Board of the National Center for Health Statistics (NCHS). Specific research programs and information are available at https://www.cdc.gov/nchs/nhanes/.

Three cycles of data from the database, 2007–2008, 2009–2010, and 2011–2012, were included in this study. We excluded from the 30,442 subjects 20,039 participants missing relevant thyroid function information, 44 participants missing SII data, 1,362 participants younger than 18 years of age, and 86 pregnant patients. This was followed by the exclusion of 430 participants missing information on marital status and education level, 759 missing information on family poverty income ratio (PIR), 100 missing information on body mass index (BMI), 575 missing information on alcohol use and smoking status, 23 missing information on glycosylated hemoglobin and alanine aminotransferase (ALT), and 149 missing information on urine iodine concentration (UIC). Since NHANES only collected data related to smoking and alcohol consumption for people 20 years of age and older, we excluded people under 20 years of age from the subject exclusion process. Ultimately, 6,875 U.S. adult participants aged 20 years or older were enrolled in the study (Figure 1).

2.2. Systemic Immune-Inflammation Index

As the independent variable in this study, SII was calculated on the basis of a complete blood count test. According to Hu et al. [14], the definition of SII is as follows: SII = platelet counts × neutrophil/lymphocyte counts, which was presented as ×109 cells/μL [22]. Complete blood counts were analyzed using automated hematology analysis instruments (Coulter® DxH 800 analyzer). Detailed laboratory methodology can be found on the NHANES website.

2.3. Thyroid Function

Thyroid function was the dependent variable in this study, which in this study included FT3, FT4, TT3, TT4, FT3/FT4, TSH, TG, TGAb, and TPOAb. Competitive binding immunoenzyme assays were used to measure serum TT4, TT3, and FT3, while two-step enzyme immunoassays were used to measure serum FT4. Using a two-site immunoenzymatic test, serum TSH levels were determined. The FT3/FT4 value was calculated to reflect the activity of peripheral T4 to T3 conversion [26]. More detailed collection and processing information can be viewed in the NHANES Laboratory/Medical Technologists Procedures Manual.

2.4. Covariates

A number of covariates potentially associated with SII and thyroid function were preselected and collected from the NHANES. Demographic covariates included age, sex, race, marital status, education level, and family PIR, which were collected through in-person family interviews. Lifestyle habit-related covariates included smoking status, alcohol use, and BMI, which were obtained from health questionnaires. Clinically relevant index covariates included UIC, glycosylated hemoglobin (HbA1c), ALT, aspartate aminotransferase (AST), creatinine (Cr), total cholesterol (TC), and the previously mentioned thyroid function, which were derived from analyses of laboratory instruments. Preexisting disease-related information, including the presence of diabetes, was obtained through questionnaires and laboratory tests. The diagnosis of diabetes in this study was made through the following criteria: (1) the participant had been diagnosed with diabetes by a physician, (2) glycosylated hemoglobin was greater than 6.5%, (3) fasting blood glucose was greater than 7 mmol/L, (4) random blood glucose or 2-hr postprandial blood glucose was greater than or equivalent to 11.1 mmol/L, and (5) insulin or diabetic medication was being used. According to the guideline published by the American Diabetes Association [27], a diagnosis of diabetes is confirmed when any of these criteria are met.

For some continuous variables, we did some regrouping. The family PIR was classified as <1, ≥3, and ≥1, <3 to assess participants’ socioeconomic status. Alcohol use was categorized according to the frequency of drinking: never: had <12 drinks in their entire lives; former: drank 12 or more times in 1 year but abstained from alcohol last year, or abstained from alcohol last year but drank 12 or more times in their entire lives; mild: <2 drinks/day; moderate: ≥2 drinks/day for women, ≥3 drinks/day for men, and ≥2 day/month for those who binge drink; and heavy: 3 drinks/day for women, ≥4 drinks/day for men, or binge drinking ≥5 days/month. Smoking status was categorized according to the frequency of smoking: never: smoked fewer than 100 cigarettes in his or her lifetime; former: consumed over 100 cigarettes in his or her lifetime and has not smoked any until now; and now: consumed over 100 cigarettes in his or her lifetime and now smoke sometimes or daily. According to WHO [28], BMI was categorized into four levels: (1) ≤18.5 kg/m2, (2) >18.5 and ≤25.0 kg/m2, (3) >25.0 and ≤30.0 kg/m2, and (4) >30 kg/m2. Urinary iodine concentrations were categorized into three groups: (1) <100 µg/L, (2) ≥100 and <300 µg/L, and (3) ≥300 µg/L.

2.5. Statistical Analysis

The statistical analyses in this study followed the recommendations of the CDC guidelines. Sampling weights were employed in the statistical analyses given the complex, multistage sampling design of the NHANES data collection. Means and standard error of mean were reported for continuous variables, whereas percentages and 95% confidence intervals (CIs) were presented for categorical variables. Utilizing survey-weighted linear regression analysis to test the differences between three SII tertile groups for continuous variables and survey-weighted chi-square testing for categorical variables. Due to the skewed distribution of the SII data, we performed a natural log transformation on SII and divided it into three groups by thirds tertiles. We first analyzed the relationship between lnSII and each thyroid function using weighted univariate linear regression analysis. Thyroid functions with statistically significant associations with lnSII were selected. We then used weighted multivariate linear regression to analyze the connection between lnSII and selected thyroid functions. To further analyze the impact of covariates on these associations, we established three models according to the recommendations of Strengthening Reporting of Observational Studies in Epidemiology (STROBE). Model 1 did not make any covariate adjustments, Model 2 adjusted for age, sex, race, marital status, family PIR, and education level, and Model 3 adjusted for all covariates. Considering the effects of BMI, UIC, and sex on thyroid function as well as conducting the sensitivity analyses [2931], we performed subgroup analyses with BMI, UIC, and sex subgroups. A weighted analysis was used in the subgroup analyses and adjusted for all covariates. To test for heterogeneity between subgroups, we performed interaction term tests. Last, we analyzed the nonlinear relationship between lnSII and thyroid function and performed sensitivity analyses using weighted restricted cubic spline (RCS) regression analyses with four slice points for modeling. R software was used for all statistical analyses, and a two-tailed p-value of less than 0.05 was considered significant.

3. Results

3.1. Baseline Characteristics of Participants

This study included 6,875 participants from NHANES 2007–2012, of which 3,424 (49.81%) were male and 3,451 (50.19%) were female. Among all participants, the mean age was 46.87 ± 0.40 years, and 4,932 (71.74%) were non-Hispanic White. After using sample weights, the 6,875 participants represented 9,104,2939 noninstitutionalized citizens in the United States. As demonstrated in Table 1, thyroid functions with statistically significant differences () between SII tertiles included FT3, FT4, TT4, FT3/FT4, and TGAb. The ranges of the three SII tertiles were, respectively, 13.750–387.529, 387.545–596.077, and 596.111–4,085. Among the covariates, the variables with statistically significant differences were sex, race, marital status, family PIR, smoking status, BMI, AST, and TC ().

3.2. The Association between Systemic Immune-Inflammation Index and Thyroid Function

We first employed the weighted univariate linear regression analysis to analyze the connection between lnSII and thyroid function (Table 2). The association between FT3, TT4, FT3/FT4, TG, TPOAb, and lnSII was found to be statistically significant () when each thyroid function was included in a univariate regression model separately.

We then included FT3, TT4, FT3/FT4, TG, and TPOAb in the weighted multivariate linear regression for analysis, respectively. The independent connection between lnSII and thyroid function was also evaluated using a variety of models. In Table 3, the outcomes of the weighted multivariate linear regression analysis are displayed. In the unadjusted model, lnSII was negatively correlated with FT3 (β = −0.0603, 95% CI −0.1078 to −0.0128, ) and FT3/FT4 (β = −0.0918, 95% CI −0.1658 to −0.0178, ). There was a positive correlation observed between lnSII and TT4 (β = 0.1837, 95% CI 0.0974–0.2700, ) and TPOAb (β = 5.1330, 95% CI 0.4263–9.8398,). The relationships between FT3 (β = −0.0559, 95% CI −0.1060 to −0.0059, ), FT3/FT4 (β = −0.0920, 95% CI −0.1667 to −0.0173, ), TT4 (β = 0.1499, 95% CI 0.0722–0.2276,), and lnSII also remained significant in the fully adjusted model.

After grouping lnSII by tertiles, lnSII remained negatively correlated with FT3 and FT3/FT4 and positively correlated with TT4. There was no statistically significant link between lnSII and TG (β = 1.6395, 95% CI −0.1964 to 3.4754, ) when presented as a continuous variable, but its effect size was more evident and statistically significant when presented as a categorical variable (Tertiles 3, β = 3.1108, 95% CI, 0.2400–5.9817, ). Continuous and categorical variables of lnSII and TPOAb were only statistically significant in the crude model, and the association between TPOAb and lnSII was not statistically significant after adjusting for confounders.

3.3. Subgroup Analysis

We conducted three subgroup analyses to examine the association between SII and thyroid function in various populations that were categorized according to different gender, BMI, and UIC, and tested the significance of the interaction term with exposure. Analysis grouped according to gender showed that the relationship between lnSII and FT3, TT4, FT3/FT4, and TG remained consistent in different genders (Table 4). The relationship between lnSII and FT3, FT3/FT4, and TG was more pronounced in the female population compared to the male and general population. There was an interaction between the association of lnSII and FT3 in the sex subgroups (). In subgroup analysis grouped according to BMI, the relationships between FT3, TT4, FT3/FT4, and lnSII remained generally consistent between the different BMI groups but were all statistically significant only in the group of people with their BMI > 25.0 and ≤30.0 kg/m2 (Table 5). TG was positively associated with lnSII in those with their BMI >30 kg/m2 and the effect size was more significant (β = 5.1916, 95% CI 1.9946–8.3887,). Interaction tests suggested an interaction between FT3, FT3/FT4, and TG in the BMI subgroups (). In subgroup analysis stratified by UIC, the relationships between FT3, TT4, FT3/FT4, and lnSII remained consistently the same (Table 6). Nevertheless, these relationships were statistically significant only in the population with their UIC ≥100 and <300 µg/L. There was an interaction between the association of lnSII and FT3 in the UIC subgroups (). As part of the sensitivity analysis, the results of the subgroup analyses indicated that the results of this study were very robust.

3.4. Weighted Restricted Cubic Spline Regression Analysis

To further discuss the nonlinear relationship between SII and thyroid function, we performed a weighted RCS analysis (Figure 2). After adjusting for confounders, the results of the weighted RCS analysis showed a negative correlation with FT3 with lnSII, and this correlation followed a nonlinear pattern (p for nonlinear = 0.008). A threshold effect can be observed. As shown in Figure 2(a), when lnSII is less than about 6.17, the curve declines rapidly with the increase of lnSII, and when lnSII is greater than about 6.17, the decline of the curve is relatively smooth. Weighted RCS analysis showed no statistically significant nonlinear relationship between TT4, FT3/FT4, TG, and lnSII (), suggesting that these relationships did not depart from linearity. Sensitivity analysis proved our results to be robust.

4. Discussion

The purpose of this study was to examine the connection between SII and thyroid function using weighted multivariable linear regression. In this cross-sectional study, we included 6,875 adults (age ≥ 20 years) in which negative correlations of lnSII with FT3 and FT3/FT4 and positive correlations of lnSII with TT4 and TG were found to be statistically significant (). We performed subgroup analyses grouped by sex, BMI, and UIC. The results showed that in female participants, FT3 and FT3/FT4 decreased more with increasing lnSII. In the overweight and obese population, the relationship between lnSII and various thyroid functions was more pronounced and statistically significant. The positive correlation between TG and lnSII was only significant in those with a BMI ≥ 30 kg/m2. The association between lnSII and thyroid function was only significant in the population with normal iodine nutritional status (UIC ≥ 100 and <300 µg/L). Interaction tests showed that the association between FT3 and lnSII interacted across subgroups. Multiple thyroid functions interacted in the obesity subgroup. Weighted RCS analysis showed a nonlinear relationship between FT3 and lnSII, with no deviation from linearity between TT4, FT3/FT4, TG, and lnSII. Multigroup sensitivity analyses demonstrated the robustness of our results. Our results indicated that SII has an important effect on several thyroid functions.

To our knowledge, this research is the first to investigate the degree of inflammation and thyroid function in vivo. Numerous studies have shown that the inflammatory response has a tremendous impact on the thyroid gland and that inflammation affects the structure of the thyroid gland. Liu et al. [10] found that thyroid nodule formation and progression may be influenced by systemic inflammation, and inflammatory markers are strongly related to the existence and size of thyroid nodules. Inflammation affected thyroid function more. Cross-sectional studies have found that proinflammatory dietary patterns are associated with higher TT4 levels [32]. C-reactive protein (CRP) is a commonly used clinical marker of nonspecific inflammation, and higher CRP levels are independently associated with lower FT3 [29, 33]. Lania et al. [34] reported that the proinflammatory cytokine IL-6 was positively associated with the risk of thyrotoxicosis in a multivariate analysis. These two studies suggest a strong link between systemic inflammation and thyroid hormone levels. Inflammatory responses due to infection are also closely linked to fluctuations in thyroid hormones. A meta-analysis based on 20 studies found that the higher the severity level of COVID-19, the greater the decrease in FT3, FT4, and TSH levels, which could be caused by the “cytokine storm” brought on by the SARS-COV-2 infection [11, 35]. At the same time, infections are environmental triggers for subsequent autoimmune thyroid diseases (AITD) and thyroid dysfunction [36]. IL-21 and IL-21R expression is upregulated in AITD, predominantly Graves’ disease and Hashimoto’s thyroiditis, and may be involved in the pathogenesis of AITD by enhancing the aberrant immune cascade, which indirectly affects thyroid function [37].

Fluctuations in thyroid hormones also have a significant impact on the inflammatory response in the body. According to Perrotta et al. [38], in vivo experiments demonstrated a significant protective effect of T3 against LPS-induced endotoxemia in mice, and in vitro experiments demonstrated this negative effect on inflammation. Thyroid hormones have a protective effect on the inflammatory process by controlling macrophage maturation and function. Contreras-Jurado et al. [39] found that thyroid hormone receptors act as potent modulators of inflammatory responses and immune homeostasis during sepsis and can directly inhibit cellular responses to inflammatory mediators. Thyroid hormone T3 inhibits IL-6 signaling in macrophages and hepatocellular carcinoma cells, suppresses STAT3 activation, and prevents the deleterious effects of excess hormone signaling during infection [39]. Various pieces of evidence have demonstrated the interaction between inflammation and thyroid function. The SII, as an indicator reflecting systemic inflammatory activity, may be able to predict changes in thyroid hormones to some extent.

The results of subgroup analyses showed that the effect of lnSII on FT3, TT4, and FT3/FT4 was more pronounced in overweight and obese populations with a BMI > 25.0 and <30.0 kg/m2 than in the general population. In accordance with our findings, numerous prior studies have demonstrated that obesity and overweight are associated with fluctuations in thyroid function. The values of FT3/FT4 reflect the activity of peripheral T4 to T3 conversion. Zupo et al. [40] found that increased skeletal muscle mass in obesity led to increased T4 to T3 conversion, especially in men. This difference may be mediated by leptin. Obese patients have been found to have high levels of circulating leptin [41, 42]. It has been found that leptin is secreted by adipocytes and promotes the expression of the gene for the thyrotropin-releasing hormone directly in the paraventricular nucleus, eventually stimulating the synthesis and secretion of TSH [43, 44]. Leptin also promotes the conversion of T4 to T3 via deiodinase, which in turn promotes an increase in FT3 [43]. Obesity or overweight affects peripheral T4 to T3 conversion activity and increases the effect of inflammation on conversion activity. Thyroglobulin is the proteogen of thyroid hormones [45], which is strongly linked to diseases such as differentiated thyroid cancer, subacute thyroiditis, and Graves’ disease. Our results suggested that the positive correlation between TG and lnSII was significant in those with their BMI ≥30 kg/m2 and that obesity may play a role in this positive correlation. Consistent with our results, Guo et al. [46] found that obesity may exacerbate the effects of autoimmunity on TSH and that there is an interaction between obesity and autoimmunity on thyroid function. Fierabracci et al. [47] reported a high prevalence of autoimmune thyroiditis in obese patients. Increased BMI is associated with adipocyte and lymphocyte infiltration of the thyroid gland, which may provide an explanation for the impairment of thyroid function [48]. Obesity is characterized by systemic tissue-specific chronic inflammation and oxidative stress [49, 50], and this environmental alteration may promote the development of AITD and exacerbate the effects of inflammation on the thyroid. The role that obesity plays in the association between inflammation and TG in vivo needs to be demonstrated by more studies.

A close relationship between SII and thyroid function in adults was demonstrated, and SII had a higher sensitivity as a predictor in obese and overweight populations. Interaction tests in the BMI subgroup analyses showed that there was an interaction of FT3, FT3/FT4, and TG with lnSII in the obesity subgroup. This indicated significant differences in the effects of SII on various thyroid functions in obese populations compared to the populations with other BMIs. Several previous studies have proved that obesity and overweight have a significant impact on inflammation and thyroid function. Our study suggested that obesity may play a mediating role in the link between SII and thyroid function, but this hypothesis still needs to be proved by more prospective experiments.

Subgroup analyses based on gender showed more significant associations between FT3, FT3/FT4, and lnSII in women compared with men and the general population. Mendelian randomization studies have found that higher levels of testosterone and estradiol in women are linked to the genetic risk score of FT4, genetically predicting that thyroid function is associated with sex hormones [51]. Estrogen therapy increases the concentration of thyroid hormone-binding globulin. Differentiated thyroid cancer is more prevalent in women than in men, and estrogen may be involved in thyroid changes by regulating enzymatic mechanisms in thyroid cells and inflammatory processes associated with tumor growth [52]. These studies provided a possible explanation for the differences in the effects of SII on thyroid function in different sexes.

There are some strengths and shortcomings in this study. We used a sizable, nationally representative sample for our investigation and considered a variety of potential confounders. The conclusions are highly generalizable. A multigroup sensitivity analysis proved the robustness of the conclusions. However, we recognized some limitations of this study. A variety of evidence has confirmed the interaction between inflammation and thyroid function, but as a cross-sectional study, we cannot draw inferences about causality. We excluded pregnant women and minor participants due to uncertainty about the effects of pregnancy and childhood on thyroid function, and future studies with pregnant women and minors should be conducted. The NHANES database lacks information on medications that could potentially impact thyroid function, and future studies should incorporate this influence. Our model has adjusted for a variety of potential confounders, but it is still unable to completely exclude the effects of other potential confounding variables. More large-scale prospective studies exploring the causal connection between SII and thyroid function are needed in the future.

5. Conclusions

In summary, our study identified significant negative correlations between SII and FT3 and FT3/FT4, and a significant positive correlation between SII and TT4. Furthermore, the association between SII and thyroid function was consistent across gender, BMI, and iodine nutritional status, but was more pronounced in women and overweight or obese populations. There were some interactions between this association in obese and overweight populations. Thus, inflammatory activity in the body may have a complex impact on thyroid disease risk, in which obesity may play an important role. Our study suggested that some measures to control the degree of inflammation in the body may be essential to maintaining normal thyroid function in our lives. Besides, to confirm this inference and validate the causal association between SII and thyroid function, additional large-scale prospective investigations are required.

Data Availability

This study used data from a publicly available datasets. These data can be found at https://www.cdc.gov/nchs/nhanes.

Conflicts 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.

Authors’ Contributions

Xin-Yu Hu and Ying-Chao Liang contributed equally to this study. All the authors contributed to the article and approved the final version.

Acknowledgments

The authors thank the staff at the National Center for Health Statistics of the Centers for Disease Control for designing, collecting, and collating the NHANES data and creating the public database. This study was funded by the Shenzhen Science and Technology Program (no. JCYJ20210324111212033), National Natural Science Foundation of China (no. 82104759) and the Natural Science Foundation of Guangdong Provincial (no. 2019A1515110108).