Data resource and participants
NHANES, a major program of the National Center for Health Statistics (NCHS), is a program of studies designed to estimate the health and nutritional status of participants in the United States. The survey combines interviews and physical examinations. We utilized 5 continuous cycles of NHANES data from year 2007 to 2016. This cross-sectional study included 29201 adult American participants aged 20 years or older. 317 subsequently pregnant participants were excluded. Furthermore, 2483 participants with incomplete kidney stones history information or other factors unrelated to drinking were excluded. Finally, we selected a sample of 10 people based on questions in the Alcohol Use Questionnaire (ALQ), named Q1-Q10. For a more in-depth analysis, we processed certain questions from different cycles of the ALQ to get samples (Q1, Q4-Q6, Q8-Q10) that are more helpful for the analysis, shown in Supplementary Figure S1. Detailed inclusion and exclusion standards are shown in Fig. 1.
Results and exposure factors
Results of the questionnaire determines whether the participant has a history of nephrolithiasis. If the answer to the question “Have you ever had kidney stones?” is “yes”, we assume that the participant has a history of kidney stones. The veracity of self-reported was confirmed in a previous study[15].
The main exposure factor was the lifetime and current (during the past 12 months) alcohol consumption of the participants, as determined by the ALQ in NHANES 2007–2016. The type of alcohol consumed was not specified.
Covariates
In order to increase accuracy and credibility, we included the following covariates: age, gender, race, marital status, education, recreational activities, smoking, asthma, overweight, gout, congestive heart failure, coronary heart disease, angina, stroke, cancer, diabetes, hypertension, and BMI. According to the recommendations of NHANES, we devided the age into three groups: 20–39, 40–59, and ≥ 60 years. Races included Mexican American, other Hispanic, non-Hispanic white, non-Hispanic black, and other races. Marital status was categorized as married, widowed, divorced, separated, never married, and living with partner. Education level included less than 9th grade, 9-11th grade (includes 12th grade with no diploma), high school graduate/GED or equivalent, some college or AA degree, and college graduate or above. BMI was divided into < 25, 25–30, and > 30 kg/m2. The comorbidities included asthma, overweight, gout, congestive heart failure, coronary heart disease, angina, stroke, cancer, diabetes, and hypertension.
Statistical methods
We used the NHANES recommended weighting data and the merging method. In the baseline characteristics table, categorical variables are expressed as proportions, and all continuous variables are treated as categorical variables, also expressed as proportions. The difference between those with and without kidney stones history was tested by survey-weighted Chi-square test.
To investigate the relationship between alcohol consumption and kidney stones prevalence, we applied three logistic regression models, including unadjusted, slightly adjusted, and fully adjusted covariates. Crude model was unadjusted. Model I was adjusted for age, gender and race. Model II was adjusted for age, gender, race, marital status, education, recreational activities, smoking, asthma, overweight, gout, congestive heart failure, coronary heart disease, angina, stroke, cancer, diabetes, hypertension, and BMI. The same method was used to examine ten different population samples separately to explore the association between alcohol consumption and kidney stones prevalence. To better assess the relationship, we treated the continuous variables in the ten samples as categorical variables for analysis and estimated trends. In addition, we performed a univariate analysis of the association of all confounders listed in the baseline tables with kidney stones.
We merged the sample weights of 5 continuous cycles according to the recommended method on the NHANES website (https://www.cdc.gov/nchs/nhanes/index.htm). According to the suggestions, we used a weight that was appropriate for the variable of interest that was collected on the smallest number of respondents. All analyses were performed using R packages (http://www.R-project.org; The R Foundation) and EmpowerStats (www.empowerstats.com, X&Y solutions Inc., Boston, MA.). A 2-tailed P < 0.05 was considered statistically significant.