3.1. Features of the participants
The individuals' clinical and demographic details are shown in Table 1. 43.22 ± 14.96 years was the mean age, with 7,615 males (47.58%). During the median follow-up period, 2,246 people (14.03%) developed NAFLD over 2.98 years. Based on the quartiles of RC (Q1: <73.74, Q2: 73.74–90.30, Q3: 90.30-109.10, Q4: ≥109.10 mg/dL), the participants were divided into subgroups. The Q4 group (≥ 109.10 mg/dL) showed higher levels for various covariates compared to the Q1 group (< 73.74 mg/dL), such as UA, age, BUN, TG, ALT, Scr, AST, FPG, TC, weight, BMI, GGT, SBP, LDL-c, and DBP, while HDL-c values and the male proportion were lower. Additionally, the proportion of females was slightly higher than that of males in each RC subgroup, particularly in Q4 compared to Q1. Figure 2 displayed the skewed distribution of RC values with a median of 32.10 mg/dL (23.59–42.92 mg/dL).
Table 1
Baseline characteristics of the participants.
RC quartiles | Q1(< 73.74) | Q2(73.74–90.30) | Q3(90.30-109.10) | Q4(≥ 109.10) | P-value |
participants | 4001 | 4001 | 4001 | 4001 | |
Age(years) | 42.32 ± 14.68 | 42.18 ± 14.62 | 43.36 ± 14.88 | 45.01 ± 15.47 | < 0.001 |
GGT(U/L) | 18.00 (14.00–26.00) | 20.00 (15.00–31.00) | 23.00 (16.00–35.00) | 26.00 (18.00–42.00) | < 0.001 |
ALT(U/L) | 15.00 (11.00-21.34) | 16.00 (11.00–23.00) | 17.00 (12.00–25.00) | 18.00 (13.00-26.87) | < 0.001 |
AST(U/L) | 21.82 ± 9.73 | 22.33 ± 8.90 | 23.11 ± 8.99 | 23.64 ± 10.17 | < 0.001 |
ALB(g/L) | 44.19 ± 2.84 | 44.32 ± 2.67 | 44.42 ± 2.74 | 44.66 ± 2.55 | < 0.001 |
GLB(g/L) | 29.35 ± 4.05 | 29.46 ± 3.86 | 29.58 ± 3.85 | 29.56 ± 3.68 | 0.025 |
TBIL(umol/L) | 12.25 ± 5.43 | 12.03 ± 4.81 | 12.21 ± 4.87 | 12.12 ± 4.78 | 0.205 |
BUN (umol/L) | 4.39 ± 1.39 | 4.44 ± 1.32 | 4.58 ± 1.31 | 4.85 ± 1.39 | < 0.001 |
Scr (umol/L) | 76.93 ± 24.56 | 78.45 ± 27.89 | 79.20 ± 23.75 | 78.95 ± 25.53 | < 0.001 |
UA (umol/L) | 256.60 ± 83.19 | 269.55 ± 80.35 | 284.98 ± 82.88 | 304.75 ± 87.42 | < 0.001 |
FPG(mmol/L) | 5.03 ± 0.67 | 5.10 ± 0.67 | 5.17 ± 0.79 | 5.26 ± 0.91 | < 0.001 |
HDL-c(mg/dL) | 90.31 ± 19.75 | 86.38 ± 20.04 | 82.92 ± 20.89 | 75.66 ± 19.07 | < 0.001 |
TC(mg/dL) | 222.32 ± 29.03 | 253.08 ± 26.84 | 276.02 ± 28.98 | 303.62 ± 29.40 | < 0.001 |
TG(mg/dL) | 77.78 ± 30.47 | 96.24 ± 39.50 | 116.25 ± 51.81 | 157.77 ± 86.99 | < 0.001 |
LDL-c(mg/dL) | 70.98 ± 14.98 | 84.63 ± 13.84 | 93.96 ± 14.03 | 100.96 ± 13.13 | < 0.001 |
Height(cm) | 163.79 ± 7.44 | 164.26 ± 7.78 | 165.03 ± 7.86 | 164.84 ± 8.03 | < 0.001 |
Weight(kg) | 55.75 ± 8.03 | 57.51 ± 8.48 | 59.14 ± 8.54 | 59.77 ± 8.49 | < 0.001 |
BMI(kg/m2) | 20.71 ± 2.05 | 21.23 ± 2.04 | 21.63 ± 1.98 | 21.91 ± 1.92 | < 0.001 |
SBP(mmHg) | 116.53 ± 15.79 | 119.15 ± 16.48 | 121.77 ± 16.23 | 125.14 ± 16.94 | < 0.001 |
DBP(mmHg) | 70.27 ± 9.65 | 71.83 ± 10.08 | 73.61 ± 10.28 | 75.34 ± 10.58 | < 0.001 |
Sex | | | | | 0.051 |
Male | 1972 (49.29%) | 1898 (47.44%) | 1896 (47.39%) | 1849 (46.21%) | |
Female | 2029 (50.71%) | 2103 (52.56%) | 2105 (52.61%) | 2152 (53.79%) | |
Values are mean ± SD or median (quartile) or n (%)
RC, Remnant cholesterol; GGT, γ-glutamyl transpeptidase; ALT, Alanine aminotransferase; AST, Aspartate aminotransferase; ALB, albumin; GLB, globulin; TBIL, total bilirubin; BUN, Serum urea nitrogen; Scr, Serum creatinine; UA, uric acid; FPG, fasting plasma glucose; HDL-c, High-density lipoprotein cholesterol; TC, Total cholesterol; TG, Triglyceride; LDL-c, Low-density lipoprotein cholesterol; BMI, Body mass index; SBP, Systolic blood pressure; DBP, Diastolic blood pressure.
3.2. The incidence rate of NAFLD
In the course of a median follow-up period of 2.98 years, Table 2 reveals that 2246 clients (14.03%) received a diagnosis of NAFLD. 500.78 cases overall, cumulatively, were reported for every 10,000 person-years. For the four subgroups with RC values, the overall prevalence rates were 236.74, 397.71, 578.55, and 766.31 instances per 10,000 person-years, respectively. NAFLD's overall rates of incidence were 14.03% (13.50-14.57%), 6.42% (5.66–7.18%), 10.72% (9.76–11.68%), 15.97% (14.84–17.11%), and 23.02% (21.716–24.32%) for each of the RC subgroups.
Table 2
RC (mg/dL) | Participants(n) | NAFLD events(n) | Incidence rate (95% CI) (%) | Per 10000 person-year |
Total | 16004 | 2246 | 14.03(13.50-14.57) | 500.78 |
Q1(< 73.74) | 4001 | 257 | 6.42 (5.66–7.18) | 236.74 |
Q2(73.74–90.30) | 4001 | 429 | 10.72 (9.76–11.68) | 397.71 |
Q3(90.30-109.10) | 4001 | 639 | 15.97(14.84–17.11) | 578.55 |
Q4(≥ 109.10) | 4001 | 921 | 23.02 (21.716–24.32) | 766.31 |
P for trend | | | < 0.001 | <0.001 |
RC, Remnant cholesterol; NAFLD, non-alcoholic fatty liver disease; CI, confidence interval.
3.3. Single-variable analysis
In the univariate assessment utilizing Cox proportional-hazards regression modeling, it was discovered that NAFLD exhibited no significant link to GLB, ALB, or TBIL, with all P-values exceeding the 0.05 threshold. However, a substantial correlation was identified between NAFLD and various factors, including age, with a hazard ratio (HR) of 1.007 and a 95% confidence interval (CI) ranging from 1.004 to 1.009. Additionally, significant associations were observed with GGT, with an HR of 1.007 and a 95% CI of 1.006 to 1.008; ALT, with an HR of 1.008 and a 95% CI from 1.007 to 1.008; and AST, with an HR of 1.011 and a 95% CI between 1.009 and 1.013. The hazard ratios for BUN, Scr, and UA were 0.932, 1.005, and 1.005, respectively, with a 95% CI that did not overlap with the null value of 1. FPG showed a notably elevated HR of 1.296, with a 95% CI from 1.266 to 1.326. HDL-c was inversely related, with an HR of 0.979 and a 95% CI of 0.977 to 0.981. TC, TG, and LDL-c exhibited positive associations with HRs of 1.004, 1.006, and 1.019, respectively, each accompanied by a 95% CI that confirmed their significance. BMI was strongly correlated with an HR of 1.815 and a 95% CI from 1.765 to 1.867. SBP and DBP also demonstrated significant relationships with HRs of 1.022 and 1.045, respectively, both with 95% CIs that were well below the significance threshold (all P-values were less than 0.05; refer to Table 3 for detailed data).
Table 3
Results of univariate Cox proportional hazards model.
Exposure | Statistics | HR (95%CI) | P-value |
Sex | | | |
Male | 7615 (47.582%) | 1.0 | |
Female | 8389 (52.418%) | 1.178 (1.084, 1.281) | < 0.001 |
Age(years) | 43.217 ± 14.958 | 1.007 (1.004, 1.009) | < 0.001 |
GGT(U/L) | 27.626 ± 30.423 | 1.007 (1.006, 1.008) | < 0.001 |
ALT(U/L) | 19.242 ± 16.515 | 1.008 (1.007, 1.008) | < 0.001 |
AST(U/L) | 22.723 ± 9.486 | 1.011 (1.009, 1.013) | < 0.001 |
ALB(g/L) | 44.399 ± 2.707 | 1.009 (0.993, 1.025) | 0.262 |
GLB(g/L | 29.491 ± 3.866 | 1.004 (0.993, 1.015) | 0.471 |
TBIL(umol/L) | 12.153 ± 4.977 | 1.001 (0.993, 1.009) | 0.843 |
BUN (umol/L) | 4.564 ± 1.364 | 0.932 (0.902, 0.963) | < 0.001 |
Scr (umol/L) | 78.384 ± 25.492 | 1.005 (1.004, 1.005) | < 0.001 |
UA (umol/L) | 278.968 ± 85.400 | 1.005 (1.005, 1.006) | < 0.001 |
FPG(mmol/L) | 5.138 ± 0.769 | 1.296 (1.266, 1.326) | < 0.001 |
HDL-c(mg/dL) | 83.818 ± 20.658 | 0.979 (0.977, 0.981) | < 0.001 |
TC(mg/dL) | 263.758 ± 41.347 | 1.004 (1.003, 1.005) | < 0.001 |
TG(mg/dL) | 112.009 ± 63.778 | 1.006 (1.006, 1.006) | < 0.001 |
LDL-c(mg/dL) | 87.632 ± 17.952 | 1.019 (1.016, 1.021) | < 0.001 |
HEIGHT | 164.480 ± 7.794 | 1.048 (1.043, 1.054) | < 0.001 |
WEIGHT1 | 58.041 ± 8.530 | 1.104 (1.099, 1.110) | < 0.001 |
BMI(kg/m2) | 21.371 ± 2.047 | 1.815 (1.765, 1.867) | < 0.001 |
SBP(mmHg) | 120.647 ± 16.673 | 1.022 (1.020, 1.024) | < 0.001 |
DBP(mmHg) | 72.764 ± 10.328 | 1.045 (1.041, 1.049) | < 0.001 |
Values are presented as the mean ± SD or n (%).
GGT, γ-glutamyl transpeptidase; ALT, Alanine aminotransferase; AST, Aspartate aminotransferase; ALB, albumin; GLB, globulin; TBIL, total bilirubin; BUN, Serum urea nitrogen; Scr, Serum creatinine; UA, uric acid; FPG, fasting plasma glucose; HDL-c, High-density lipoprotein cholesterol; TC, Total cholesterol; TG, Triglyceride; LDL-c, Low-density lipid cholesterol; BMI, Body mass index; SBP, Systolic blood pressure; DBP, Diastolic blood pressure; HR, Hazard ratios; CI, confidence interval.
Figure 3 displays Kaplan-Meier survival curves that illustrate the likelihood of NAFLD-free survival categorized by RC subgroups. P < 0.0001 for the log-rank test indicated significant differences in the NAFLD-free survival probabilities across the RC subgroups. As the RC level increased, the probability of NAFLD-free survival progressively declined, suggesting that the group with the highest RC level also had the highest risk of NAFLD development.
3.4. Multivariate analysis
We developed a trio of Cox proportional hazards regression models to scrutinize the correlation between RC levels and the incidence of NAFLD, as detailed in Table 4. The initial, unadjusted model (referred to as Model I) revealed that for every 10 mg/dL elevation in RC, there was a corresponding 14.2% escalation in NAFLD risk, with an HR of 1.142, a 95% CI ranging from 1.125 to 1.159, and a P-value well below 0.001. Upon introducing minimal adjustments for demographic factors—namely age and sex—in Model II, the risk increment associated with a 10 mg/dL rise in RC was slightly reduced to 13.9%, yielding an HR of 1.139 and a 95% CI from 1.122 to 1.156, with statistical significance persisting (P < 0.001). In the comprehensively adjusted model, Model III, which accounted for a broader spectrum of variables including SBP, AST, sex, BMI, Scr, GGT, ALB, ALT, FPG, UA, age, and DBP, each 10 mg/dL increment in RC corresponded to a 6.1% heightened risk of NAFLD. This was represented by an HR of 1.061 and a 95% CI spanning 1.045 to 1.078, with the relationship maintaining robust statistical support (P < 0.001). The consistency of the confidence intervals across models suggests that the observed association between RC and NAFLD is statistically robust.
Table 4
Association between RC and incident NAFLD in different models.
Exposure | Model I (HR,95%CI, P) | Model II (HR,95%CI, P) | Model III (HR,95%CI, P) | Model IV (HR,95%CI, P) |
RC (per 10mg/dL) | 1.142 (1.125, 1.159) < 0.001 | 1.139 (1.122, 1.156) < 0.001 | 1.061 (1.045, 1.078) < 0.001 | 1.058 (1.041, 1.076) < 0.001 |
RC quartiles | | | | |
Q1 | Ref | Ref | Ref | 1.0 |
Q2 | 1.635 (1.401, 1.908) < 0.001 | 1.634 (1.400, 1.907) < 0.001 | 1.292 (1.106, 1.508) 0.001 | 1.197 (1.024, 1.398) 0.024 |
Q3 | 2.320 (2.007, 2.682) < 0.001 | 2.305 (1.994, 2.665) < 0.001 | 1.576 (1.362, 1.823) < 0.001 | 1.408 (1.216, 1.631) < 0.001 |
Q4 | 2.882 (2.507, 3.312) < 0.001 | 2.840 (2.471, 3.265) < 0.001 | 1.672 (1.452, 1.926) < 0.001 | 1.560 (1.352, 1.799) < 0.001 |
P for trend | < 0.001 | < 0.001 | < 0.001 | < 0.001 |
Model I: We did not adjust for the other covariates.
Model II: Adjusted for age and sex.
Model III: We adjusted for age, sex, SBP, BMI, DBP, ALT, GGT, ALB, FPG, AST, UA, and Scr levels.
Model IV: Adjusted for age(smooth), sex, SBP (smooth), BMI(smooth), DBP(smooth), ALT(smooth), GGT(smooth), ALB(smooth), FPG(smooth), AST(smooth), UA(smooth), and Scr (smooth).
HR, Hazard ratio; CI, confidence interval; Ref, reference; NAFLD, non-alcoholic fatty liver disease.
3.5. Sensitivity analysis
To validate the robustness of our findings, sensitivity analyses were conducted. Initially, the RC levels were divided into quartiles and reintegrated into the regression model. The analysis confirmed that HRs remained stable across quartiles, mirroring the impact observed with the continuous RC variable, and the trend's P-value was congruent with the original analysis. As depicted in Model IV of Table 4, the results were in line with the fully adjusted model, which incorporated smooth terms for Scr, UA, sex, SBP, AST, FPG, DBP, GGT, ALT, ALB, BMI, and age, yielding an HR of 1.058 with a 95% CI of 1.041–1.076 and a P-value less than 0.001. Further sensitivity analyses excluded individuals with FPG exceeding 7.0 mmol/L. Even after controlling for potential confounding factors—comprising BMI, age, ALT, UA, SBP, AST, DBP, sex, GGT, ALB, FPG, and Scr—the relationship between RC and the risk of NAFLD was still positive, with an HR of 1.061 and a 95% CI of 1.044–1.078, achieving statistical significance (P < 0.001) as documented in Table 5. Parallel sensitivity analyses, which omitted participants diagnosed with hypertension (defined by SBP of 140 mmHg or higher or DBP of 90 mmHg or higher), retained the significant association between RC and NAFLD risk after adjustment for the same set of covariates, with an HR of 1.062 and a 95% CI of 1.043–1.082 and a P-value below 0.001. The uniformity of these outcomes across all sensitivity analyses in Table 5 reinforces the reliability of our conclusions.
Table 5
Relationship between RC and NAFLD in different sensitivity analyses
Exposure | Model V (HR,95%CI, P) | Model VI (HR,95%CI, P) |
RC (per 10mg/dL) | 1.061 (1.044, 1.078) < 0.001 | 1.062 (1.043, 1.082) < 0.001 |
RC quartiles | | |
Q1 | Ref | Ref |
Q2 | 1.300 (1.109, 1.524) 0.001 | 1.279 (1.071, 1.528) 0.007 |
Q3 | 1.562 (1.344, 1.814) < 0.001 | 1.528 (1.292, 1.808) < 0.001 |
Q4 | 1.664 (1.439, 1.925) < 0.001 | 1.672 (1.418, 1.973) < 0.001 |
P for trend | < 0.001 | < 0.001 |
Model V was a sensitivity analysis for participants without FPG > 7.0 mmol/L (N = 15683). We adjusted for age, sex, SBP, BMI, DBP, ALT, GGT, ALB, FPG, AST, UA, and Scr.
Model VI was a sensitivity analysis in participants without hypertension, SBP ≥ 180 mmHg, or DBP ≥ 90 mmHg (N = 13622). We adjusted for age, sex, SBP, BMI, DBP, ALT, GGT, ALB, FPG, AST, UA, and Scr.
HR, Hazard ratio; CI, confidence interval; Ref, reference; eGFR, evaluated glomerular filtration rate (mL/min·1.73 m2); NAFLD, non-alcoholic fatty liver disease.
3.6. Nonlinear relationship analysis
Utilizing a Cox proportional hazards regression model complemented by cubic spline functions, we discerned a nonlinear association between RC levels and the prevalence of NAFLD, as depicted in Fig. 4. To dissect this nonlinearity, a piecewise binary logistic regression modeling was engaged to delineate two separate gradients. This approach was juxtaposed with a conventional binary logistic regression model through the log-likelihood ratio test, as presented in Table 6, with the outcome indicating statistical significance (P < 0.001). Employing a recursive strategy, the critical inflection point of RC within the data was pinpointed at 98.29 mg/dL. Subsequently, a two-piecewise Cox proportional hazards regression modeling was applied, taking into account covariates such as UA, DBP, age, AST, SBP, sex, ALT, BMI, ALB, GGT, FPG, and Scr. This model was utilized to gauge the HR and its CI on either side of the identified inflection point. Our findings revealed a positive link between RC and NAFLD on the lower side of the curve's inflection point, as illustrated in Fig. 4. Specifically, the HR was 1.150 on the left (95% CI: 1.106–1.194), contrasting with an HR of 1.009 on the right side (95% CI: 0.982–1.037), showcasing a differential impact of RC levels on NAFLD risk across various thresholds.
Table 6
The result of the two-piecewise Cox regression model
Incident NAFLD | HR (95%CI) P-value |
Fitting model by standard Cox regression | 1.061 (1.045, 1.078) < 0.001 |
Fitting model by two-piecewise Cox regression | |
The inflection point of RC | 98.29 |
< 98.29(per 10mg/dL) | 1.150 (1.106, 1.194) < 0.001 |
≥ 98.29 (per 10mg/dL) | 1.009 (0.982, 1.037) 0.521 |
P for log-likelihood ratio test | < 0.001 |
HR, Hazard ratio; CI, confidence interval; RC, remnant cholesterol. We adjusted for age, sex, SBP, BMI, DBP, ALT, GGT, ALB, FPG, AST, UA, and Scr levels.