Study Design and Participants
The survey data utilized in this study originated from individuals who enrolled in the U.S. National Health and Nutrition Examination Survey (NHANES) during the cycle of January 2017 to March 2020. The NHANES is an ongoing epidemiological survey program conducted by the Centers for Disease Control and Prevention (CDC) to assess the health and nutritional conditions of the U.S. population [24]. It relies on a complex multi-stage probability sampling strategy to collect information from a nationally representative sample through standardized household interviews, health screenings, and laboratory tests [25]. The survey was approved by the Ethics Review Board of the National Center for Health Statistics (NCHS), and all participants voluntarily participated and signed an informed consent form. More information is available free of charge on the NHANES website (https://www.cdc.gov/nchs/nhanes/index.htm).
A total of 15,560 participants were enrolled in the 2017–2020 cycle of NHANES. We sequentially excluded individuals who were younger than 20 years, pregnant, lacked information on nine body measurements (BMI, WC, LAP, VAI, BRI, WWI, CMI, AIP, TYG), missing information on gallstones, and missing information on covariates. A total of 3698 subjects were ultimately selected for analysis. The detailed process of participant screening is illustrated in Fig. 1.
Definition of 9 obesity-related screening indicators and gallstones:
In this study, the 9 screening measures associated with obesity include BMI, WC, LAP, VAI, BRI, WWI, CMI, AIP, and TyG. The formulas for their calculation are summarized in Fig. 2. Data related to body measurements (BMI, weight, height, and WC) were collected by trained professional technicians at the mobile examination center according to standardized protocols. In addition, HDL-C, TG, and FPG data, which are the other indicators required in the formula, were analyzed by the University of Minnesota laboratory under strict surveillance and quality management. HDL-C, FPG, and TG were measured by turbidimetric immunoassay, hexokinase (HK) assay, and automated direct chemiluminescence analyzer, respectively.
For the definition of gallstones, participants were interviewed by a well-trained specialized interviewer who asked the question, "Has a doctor or other health professional ever diagnosed you with gallstones?" The answer to this question was the basis for determining whether or not a person suffered from gallstones. This convenient and rapid approach has been utilized in previous studies [23, 26].
Covariates:
To explore the association between the 9 measures and gallstones, we considered several confounders to be incorporated as covariates. Demographic information includes: age (≤ 40, over 40), sex, race (Mexican American, non-Hispanic white, non-Hispanic black, etc.), marital status, and education level (less than high school, high school or equivalent, and college or above). For smoking status, we categorized them into three groups depending on the self-reported questionnaire: never smokers, former smokers, and current smokers. Hypertension was recognized when any of the following conditions were satisfied: 1. Three consecutive diastolic pressures ≥ 90 mmHg and/or systolic pressures ≥ 140 mmHg. 2. self-reported or informed by a healthcare provider that they had hypertension. Diabetes mellitus was diagnosed according to the following criteria: glycosylated hemoglobin ≥ 6.5% or fasting blood glucose level ≥ 7.0 mmol/L or self-reported diagnosis of diabetes by a healthcare provider. Hypercholesterolemia was based on laboratory-measured cholesterol levels (≥ 6.2 mmol/L). More detailed descriptions of the above-mentioned covariates are publicly accessible online at the NHANES website.
Statistical analysis
All data analyses for this study were performed in strict accordance with CDC guidelines. To provide nationally representative effect estimates, we applied the recommended sample weights (WTSAFPP). Different presentation methods were adopted depending on the type of data, with continuous variables expressed as weighted means ± standard error (SE) and categorical variables as weighted proportions ± SE. In the comparison of patients with non-gallstones and gallstones, weighted linear regression was employed for the continuous variables, and weighted chi-square tests were applied to the categorical variables. To explore the relationship between 9 obesity-related screening indicators and gallstones, weighted logistic regression was applied to analyze their performance in 3 different models. Model 1 remained unadjusted, model 2 was corrected for age, sex, and race, and model 3 was further adjusted for marital status, educational attainment, smoking, hypertension, diabetes, and hypercholesterolemia. To determine the robustness of the results, we analyzed the independent variables in both continuous and categorical forms for their association with gallstone disease. We also plotted subjects' work characteristics (ROC) curves and computed the area under the curve (AUC) values to evaluate the predictive power of the different indicators. Furthermore, bootstrap resampling (carried out 500 times) was conducted as a sensitivity analysis to assess the stability of AUC. After the identification of the two most predictive indicators, restrictive cubic bars were developed to visualize their association with the risk of gallstones. A comparison of the variability of AUC between the different screening parameters was performed using the Delong test. Statistical differences existed when the p-value obtained from the significance test was less than 0.05. All analyses were performed using EmpowerStats and R4.3.2.