DOI: https://doi.org/10.21203/rs.3.rs-22525/v1
Evidence regarding the relationship between normal low-density lipoprotein cholesterol (LDL-c) levels and non-alcoholic fatty liver disease (NAFLD) was limited. Therefore, this dissertation seeks to investigate the relationship between LDL-c and NAFLD in the non-obese Chinese population after adjusting for other covariates.
The present study was a cross-sectional study. A total of 183903 non-obese participants were involved in a Wenzhou Medical Center of Wenzhou People’s Hospital from 2010 to 2014. The target independent variable and the dependent variable were LDL-c measured at baseline and NAFLD respectively. Covariates involved in this study included SEX, AGE, BMI, SBP, DBP, FPG, ALB, ALT, AST, BUN, Cr, TG, TC, HDL-c, UA. It was noted that the entire study was completed by Dan-Qin Sun et al., and uploaded the data to the DATADRYAD website. The author only used this data for secondary analysis.
The average age of 183903 selected participants was 41.0 ± 14.1 years old,and about 49.6% of them was male. After adjusting potential confounders (SEX, AGE, BMI, FPG, ALB, GLB, ALT, AST, GGT, BUN, Cr, TG, TC, HDL-c, UA), non-linear relationship was detected between normal LDL-c levels and NAFLD, whose point was 1.51. The effect sizes and the confidence intervals on the left and right sides of inflection point were 0.87 (0.64, 1.18) and 1.79 (1.67, 1.92), respectively.
The relationship between normal LDL-c levels and NAFLD is non-linear. Normal LDL-c levels was positively correlated with NAFLD when LDL-c was more than 1.51.
Non-alcoholic fatty liver disease (NAFLD) results from hepatic fat accumulation in the absence of quantities of alcohol and any secondary cause[1–5]. Histologically, It comprises a range of pathologic conditions including simple nonalcoholic steatosis, nonalcoholic steatohepatitis (NASH) and hepati cirrhosis[6, 7]. NAFLD is a major, worldwide public health problem and its incidence among adults and adolescents has increased rapidly in recent years. Salvoza etal[8]found that the prevalence of NAFLD is 30–40% in the United States, 2–44% in Europe and 15–45% in Asia, respectively.
One meta-analysis reports that it is the relatively high prevalence of NAFLD found in the Asian population (27%)[9–11]. The prevalence rate of NAFLD in the general population of China varies from 24.77–43.91% in recent years[11, 12].
NAFLD is associated with an increased risk of developing to severe and decompensated liver diseases such as cirrhosis and hepatocellular carcinoma[13–15]. In addition, NAFLD also can increase the risk of cardiovascular disease, type 2 diabetes and chronic kidney disease[16–19].
Obesity is a well-known risk factor for the development of NAFLD, and NAFLD is now considered to be one of the most prevalent manifestations of obesity-related metabolic syndrome (MS)[20]. However, there is growing evidence that the prevalence of NAFLD in non-obese individuals is not uncommon. Especially in Asia, where people are less obese than in Western countries and the prevalence of NAFLD is increasing with time[21]. Studies have shown that approximately 15–21% of the Asian patients with NAFLD are non-obese.
Dyslipidemia is a comorbidity of NAFLD, which can lead to hypertriglyceridaemia, reductions in high-density lipoprotein cholesterol (HDL-c), an increase in the size of very low-density lipoprotein and low-density lipoprotein cholesterol (LDL-c)[22]. In addition, more and more studies have shown that LDL-c is involved in the development of NAFLD and non-alcoholic steatohepatitis (NASH)[23–25]. This mechanism is still unclear, but insulin resistance and disorders of fat and sugar metabolism may explain some reasons[26–28]. In addition, damage to the LRP6 receptor is also responsible[29, 30].
In any case, so far, no studies have explored the curves relationship between LDL-c and NAFLD. In biomedical research, connection between exposures and outcomes may be non-linear. In that case, researchers need a more effective method to deal with non-linear relationship. In this study, we use curves to fit LDL-c and NAFLD.
We obtained data from ‘DATADRYAD’ database (www.Datadryad.org). This website permitted users to freely download the raw data. According to Drya Terms of Service, we cited Dryad data package in the present study. (Sun, Dan-Qin; Wu, Sheng-Jie; Liu, Wen-Yue; Wang, Li-Ren; Chen, Yi-Ran; Zhang, Dong-Chu; et al. (2016): Data from: Association of low-density lipoprotein cholesterol within the normal range and NAFLD in the non-obese Chinese population: a cross-sectional and longitudinal study. Dryad. Collection. https://doi.org/10.5061/dryad.1n6c4). Variables included in the database file were as follows: age, sex, body mass index (BMI), fasting plasma glucose(FPG), albumin (ALB), globulin (GLB), alanine aminotransferase (ALT), aspartate transaminase (AST), γ-glutamyltranspeptidase (GGT), blood urea nitrogen (BUN), creatinine (CR), triglyceride (TG), total cholesterol (TC), high-density lipoprotein cholesterol (HDL-c), uric acid (UA), total bilirubin (TBIL), low density lipoprotein (LDL-c), non-alcoholic fatty liver disease (NAFLD).
Sun D-Q, Wu S-J, Liu W-Y, et al[31] completed the entire study. In order to allow to understand the entire research process more clearly, we have outlined the steps of the study here. The specific details are described in the original Reported by Sun D-Q, Wu S-J, Liu W-Y, et al. They conducted a cross-sectional study at Wenzhou Medical Center of Wenzhou People’s Hospital, zhejiang city, china from 2010 to 2014. Participants received a medical health examination including the biochemical measurements and abdominal ultrasonography. A total of 183903 non-obese individuals were enrolled and selected according to exclusion standards. Exclusion standards: (1) Individuals were reported excess alcohol consumption (> 140 g/week for men and > 70 g/week for women), (2) those taking antihypertensive agents, antidiabetic agents, lipid-lowing agents, (3) other known causes of chronic liver disease, (4) LDL-c of > 3.12 mmol/L, (5) BMI of ≥ 25 kg/m2, (6) Missing data. Researchers obtained information (values) of NAFLD, LDL-c and other covariants at baseline. This study was performed by Chinese researcher Sun D-Q, Wu S-J,Liu W-Y,et al. at Wenzhou Medical Center of Wenzhou People’s Hospital, zhejiang city, china. In the previously published article [31, 32], Sun D-Q, Wu S-J,Liu W-Y, et al. has clearly stated that: Verbal informed consent was obtained from each participant before their participation in the study. The personal information of participants was erased and replaced by the health examination number. The research protocol of the study was approved by the ethics committee of Wenzhou People’s Hospital. Informed consent was obtained from all Participants.
The biochemical measurements included albumin (ALB), alanine aminotransferase (ALT), aspartate aminotransferase (AST), fasting plasma glucose (FPG), globulin (GLB), γ-glutamyltranspeptidase (GGT), blood urea nitrogen (BUN), creatinine (Cr), uric acid (UA), total bilirubin (TBIL), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-c) and low-density lipoprotein cholesterol (LDL-c). All values were measured by an automated analyser (Abbott AxSYM) using standard methods. Medical history and a health habit inventory were taken by a physician. BMI (kg/m2), used as an index of body fat, was calculated as weight in kilograms divided by height in m2. The ultrasound criteria for the diagnosis of fatty liver were based on those suggested by the Chinese Liver Disease Association. The assessment of NAFLD described in detail in the original.
Continuous variables were expressed as mean ± standard deviation (normal distribution) or median (quartile) (skewed distribution), and categorical variables were expressed in frequency or as a percentage. The One-Way Anova (normal istribution), Kruscal Whallis H (skewed distribution) test and chi-square tests (categorical variables) were used to determine any statistical differences between the means and proportions of the groups. Univariate linear regression model was used to evaluate the associations between LDL-c and NAFLD. Both non-adjusted and multivariate adjusted models were listed in the paper. According to the recommendation of STROBE statement, we simultaneously showed the results of unadjusted, minimally adjusted analyses and those from fully adjusted analyses. Whether the covariances were adjusted determined by the following principle: when added to this model, changed the matched odds ratio by at least 10% [33]. Besides, we also used generalized additive model (GAM) to identify the non-linear relationship. If the non-linear correlation was observed, a two-piecewise linear regression model was performed to calculate the threshold effect of the LDL-c and NAFLD in terms of the smoothing plot. When the ratio between LDL-c and NAFLD appears obvious in smoothed curve, recursive method calculates automatically the inflection point, where the maximum model likelihood will be used[34]. All of the analyses were performed with the statistical software packages R (http://www. R-project.org, The R Foundation) and EmpowerStats (http://www.empowerstats.com, X&Y Solutions, Inc., Boston, MA). P values less than 0.05 (two-sided) were considered statistically significant.
A total of 183903 participants were sclected for the final data analysis based on the inclusion and exclusion criteria. We showed baseline characteristics of these selected participants in Table 1 according to Quartile of LDL-c. In general, the average age of the 183903 selected participants was 41.0 ± 14.1 years old, and about 49.6% of them were male. statistically significant differences were detected in AGE, BMI, HDL, ALT, AST, GGT, BUN, Cr, TG, FPG, TC, UA, GLB and ALB, among different LDL-c groups (all p values < 0.01). Participants wih the highest group of LDL-c (Q4) had the higher values in AGE, BMI, ALT, AST, GGT, BUN, CR, TG, FPG, TC, UA, GLB and ALB than those of the other groups. The opposite pattern were observed in HDL-c.
LDL-c | Q1 (0.41–1.91) | Q2 (1.92–2.26) | Q3 (2.27–2.60) | Q4 (2.61–3.12) | P-value |
---|---|---|---|---|---|
Number | 45352 | 45917 | 45679 | 46952 | |
Age, years | 38.0 ± 13.9 | 39.3 ± 13.4 | 41.7 ± 13.7 | 44.8 ± 14.2 | < 0.001 |
BMI, kg/m2 | 20.8 ± 2.2 | 21.2 ± 2.1 | 21.7 ± 2.1 | 22.1 ± 1.9 | < 0.001 |
FPG log2 transform | 2.3 ± 0.2 | 2.3 ± 0.2 | 2.4 ± 0.2 | 2.4 ± 0.2 | < 0.001 |
ALB, g/L | 44.4 ± 2.9 | 44.6 ± 2.8 | 44.7 ± 2.8 | 44.7 ± 2.8 | < 0.001 |
GLB, g/L | 29.3 ± 4.0 | 29.3 ± 3.8 | 29.3 ± 3.8 | 29.5 ± 3.8 | < 0.001 |
ALT log2 transform | 4.0 ± 0.7 | 4.1 ± 0.7 | 4.2 ± 0.8 | 4.3 ± 0.8 | < 0.001 |
AST log2 transform | 4.4 ± 0.4 | 4.4 ± 0.4 | 4.4 ± 0.4 | 4.5 ± 0.4 | < 0.001 |
GGT log2 transform | 4.3 ± 0.8 | 4.5 ± 0.8 | 4.6 ± 0.9 | 4.7 ± 0.9 | < 0.001 |
BUN log2 transform | 2.0 ± 0.4 | 2.1 ± 0.4 | 2.1 ± 0.4 | 2.1 ± 0.4 | < 0.001 |
Cr log2 transform | 6.2 ± 0.3 | 6.2 ± 0.3 | 6.3 ± 0.3 | 6.3 ± 0.3 | < 0.001 |
TG log2 transform | -0.1 ± 0.8 | 0.1 ± 0.7 | 0.3 ± 0.7 | 0.4 ± 0.7 | < 0.001 |
TC, mmol/L | 3.8 ± 0.5 | 4.3 ± 0.5 | 4.8 ± 0.5 | 5.3 ± 0.5 | < 0.001 |
HDL-c, mmol/L | 1.5 ± 0.4 | 1.4 ± 0.4 | 1.4 ± 0.3 | 1.4 ± 0.3 | < 0.001 |
UA, µmol/L | 264.5 ± 85.9 | 277.1 ± 87.6 | 289.5 ± 88.7 | 300.9 ± 89.1 | < 0.001 |
TBIL, µmol/L | 12.4 ± 5.4 | 12.3 ± 5.0 | 12.3 ± 5.0 | 12.2 ± 4.8 | 0.815 |
Sex,n(%) | < 0.001 | ||||
female | 27221 (60.0%) | 24576 (53.5%) | 21431 (46.9%) | 19412 (41.3%) | |
male | 18131 (40.0%) | 21341 (46.5%) | 24248 (53.1%) | 27540 (58.7%) | |
NAFLD,n(%) | < 0.001 | ||||
no | 41736 (92.0%) | 40988 (89.3%) | 38582 (84.5%) | 37091 (79.0%) | |
yes | 3616 ( 8.0%) | 4929 (10.7%) | 7097 (15.5%) | 9861 (21.0%) | |
ALB, albumin; GLB, globulin, GGT, γ-glutamyltranspeptidase, ALT, alanine aminotransferase; AST, aspartate aminotransferase; BMI, body mass index; BUN, blood urea nitrogen; Cr, creatinine; FPG, fasting plasma glucose; HDL-c, high-density lipoprotein cholesterol; LDL-c, low-density lipoprotein cholesterol; NAFLD, non-alcoholic fatty liver disease; TC, total cholesterol; TG, triglyceride; UA, uric acid; TBIL, total bilirubin. |
We listed the results of univariate analyses in Table 2. The results of univariate analysis showed that AGE, BMI, ALT, AST, GGT, BUN, CR, TG, FPG, TC, UA, TBIL, GLB and ALB were correlated with NAFLD. We found that AGE 1.03 (1.03, 1.03), BMI 2.08 (2.06, 2.11), FPG log2 transform 2.53 (2.48, 2.59), ALB 1.07 (1.06, 1.07), GLB 1.01 (1.01, 1.01), ALT log2 transform 3.05 (2.99, 3.11), AST log2 transform 3.18 (3.09, 3.28), GGT log2 transform 2.72 (2.68, 2.77), BUN log2 transform 1.87 (1.81, 1.93), Cr log2 transform 3.33 (3.20, 3.47), TG log2 transform 4.67 (4.57, 4.76), TC 1.74 (1.71, 1.77), and UA 1.01 (1.01, 1.01) are risk factors for NAFLD, however HDL-c 0.11 (0.10, 0.11) is a protective factor for NAFLD.
Covariate | Statistics | OR (95%CI) | P-value |
---|---|---|---|
Age, years | 40.98 ± 14.06 | 1.03 (1.03, 1.03) | < 0.001 |
BMI, kg/m2 | 21.43 ± 2.13 | 2.08 (2.06, 2.11) | < 0.001 |
FPG log2 transform | 5.08 ± 0.53 | 2.53 (2.48, 2.59) | < 0.001 |
ALB, U/L | 44.59 ± 2.84 | 1.07 (1.06, 1.07) | < 0.001 |
GLB, U/L | 29.34 ± 3.86 | 1.01 (1.01, 1.01) | < 0.001 |
ALT log2 transform | 4.12 ± 0.76 | 3.05 (2.99, 3.11) | < 0.001 |
AST log2 transform | 4.43 ± 0.42 | 3.18 (3.09, 3.28) | < 0.001 |
GGT log2 transform | 4.54 ± 0.86 | 2.72 (2.68, 2.77) | < 0.001 |
BUN log2 transform | 2.08 ± 0.40 | 1.87 (1.81, 1.93) | < 0.001 |
Cr log2 transform | 6.26 ± 0.33 | 3.33 (3.20, 3.47) | < 0.001 |
TG log2 transform | 0.19 ± 0.76 | 4.67 (4.57, 4.76) | < 0.001 |
TC, mmol/L | 4.54 ± 0.74 | 1.74 (1.71, 1.77) | < 0.001 |
HDL-c, mmol/L | 1.45 ± 0.36 | 0.11 (0.10, 0.11) | < 0.001 |
UA, µmol/L | 283.13 ± 88.90 | 1.01 (1.01, 1.01) | < 0.001 |
TBIL, µmol/L | 12.30 ± 5.04 | 1.01 (1.01, 1.01) | < 0.001 |
LDL-c, mmol/L | 2.25 ± 0.47 | 2.54 (2.46, 2.61) | < 0.001 |
LDL-c (Quartile) | < 0.001 | ||
Q1 | 45352 (24.66%) | 1.0 | |
Q2 | 45917 (24.97%) | 1.39 (1.33, 1.45) | < 0.001 |
Q3 | 45679 (24.84%) | 2.12 (2.04, 2.22) | < 0.001 |
Q4 | 46952 (25.53%) | 3.07 (2.95, 3.20) | < 0.001 |
Sex,n(%) | < 0.001 | ||
female | 92643 (50.38%) | 1.0 | |
male | 91260 (49.62%) | 4.22 (4.09, 4.35) |
In this study, We used Multivariate linear regression model to evaluate the associations between LDL-c and NAFLD. Meanwhile,we showed the non-adjusted and adjusted modles in Table 3. In crude model, LDL-c showed positive correlation with NAFLD (OR = 2.54, 95% confidence interval (CI): 2.46 to 2.61, P < 0.001). In minimally adjusted model (adjusted age, sex), the result have obvious changes (OR = 2.01, 95% confidence interval (CI): 1.94 to 2.07, P < 0.001). However,we detect the Independent effect LDL-c on NAFLD in fully adjusted model (OR = 1.70, 95% confidence interval (CI): 1.59 to 1.82, P < 0.001). For the purpose of sensitivity analysis, we also handled LDL-c as Categorical variable (Quartile), and found that the same trend was observed as well (p for trend was less than 0.001). Besides,we also found the trend of the effect size in different model was equidistant.
Exposure | Non-adjusted | Adjust I | Adjust II |
---|---|---|---|
LDL-c | 2.54 (2.46, 2.61) < 0.001 | 2.01 (1.94, 2.07) < 0.001 | 1.70 (1.59, 1.82) < 0.001 |
LDL-c (quartile) | |||
Q1 | Ref | Ref | Ref |
Q2 | 1.39 (1.33, 1.45) < 0.001 | 1.28 (1.23, 1.35) < 0.001 | 1.18 (1.10, 1.25) < 0.001 |
Q3 | 2.12 (2.04, 2.22) < 0.001 | 1.78 (1.70, 1.86) < 0.001 | 1.47 (1.37, 1.58) < 0.001 |
Q4 | 3.07 (2.95, 3.20) < 0.001 | 2.33 (2.23, 2.43) < 0.001 | 1.76 (1.62, 1.91) < 0.001 |
P for trend | < 0.001 | < 0.001 | < 0.001 |
Crude model: we did not adjust other covariants | |||
Minimally adjusted Model : we adjusted age and sex | |||
Fully adjusted Model : we adjusted AGE, SEX, BMI, ALT, AST, GGT, BUN, Cr, TG, FPG, TC, UA, HDL-c, GLB and ALB. CI, confidence interval; OR, odds ratio |
Because LDL-c was continuous variable, the analyses of non-linear relationship are necessary. In the present study (Fig. 1), we found that the relationship between LDL-c and NAFLD was non-linear (after adjusting AGE, SEX, BMI, ALT, AST, GGT, BUN, Cr, TG, HDL-c, FPG, TC, UA, GLB and ALB). By two-piecewise linear regression model, we calculated that the inflection point was 1.51. On the left side of the inflection point, the effect size, 95% CI and P value were 0.87, 0.64 to 1.18 and 0.3760, respectively. However, On the right side of the inflection point, the effect size, 95% CI and P value were 1.79, 1.67 to 1.92 and p < 0.0001, respectively (Table 4).
Table 4 The results of two-piecewise linear regression model
Inflection point of LDL-c |
OR |
95%CI |
P value |
<1.51 |
0.87 |
(0.64, 1.18) |
0.3760 |
≥1.51 |
1.79 |
(1.67, 1.92) |
<0.0001 |
Effect: NAFLD Cause: LDL-c Adjusted: SEX; AGE; BMI; ALB; GLB; ALT; AST; GGT; BUN; Cr; FPG; HDL-c; TG; TC; UA.
Previous studies showed that LDL-c was associated with NAFLD, but the relationship between the risk factor LDL-c and NAFLD was not fully described. We conducted a PubMed search and three scientific papers were retrieved from the database as of the end of March 2020. All of these studies showed that LDL-c was associated with NAFLD[24, 31, 32]. Sun D-Q, Wu S-J,Liu W-Y, et al. found that LDL-c was associated with NAFLD in the non-obese Chinese population in a cross-sectional and longitudinal study[31, 32]. In their studies, the multivariable logistic regression models were used to calculate the OR of LDL-c on NAFLD. After adjusting potential confounders (SEX, AGE, BMI, DBP, SBP, ALB, ALT, AST, BUN, Cr, FPG, HDL-c, TC, TG and UA), the OR gradually increased in Q1 to Q4 of LDL-c quartile, and the P for trend was less than 0.05. This suggests that the connection between LDL-c and NAFLD is non-linear. However, none of them discussed the non-linear connection between LDL-c and NAFLD. To our knowledge, this is the first study to investigate the non-linear relationship between LDL-c and NAFLD.
In the present study, we used GLM and GAM models to elucidate the relationship between LDL-c and NAFLD among participants. As is shown in the fully adjusted model, LDL-c was associated with NAFLD. When we handled LDL-c as a categorical variable, the same trend was observed. However, the results obtained from GAM and two-piecewise linear regression model showed that the relationship between LDL-c and NAFLD was non-linear, and the correlations between LDL-c and NAFLD were different on the left and right sides of the inflection point (LDL-c = 1.51). LDL-c, as assessed at baseline, was not statistically significant on the left side of the inflection point, but LDL-c was positively associated with NAFLD on the right of the inflection point.
Our study has a number of strengths
(1) we not only use the generalized linear model to evaluate the linear relationship between LDL-c and NAFLD, but also use the generalized additive model to clarify the nonlinear relationship. GAM has obvious advantages in dealing with non-linear relations and it can handle the non-parametric smoothing and will fit a regression spline to the data. The use of GAM will help us to better discover the real relationships between exposure and outcome. we detect this nonlinear relationship in our study after adjusting confounding factors which were not founded by previous study. (2) this study is an observational study including unavoidable potential confounding, so we used strict statistical adjustment to minimize residual confounding.
There are some limitations in our study
(1) this study is a analytical cross-sectional study and therefore only provides weak evidence between exposure and outcome, and it is difficult to distinguish the cause and effect. (2) because the study population contains only Chinese, it may be not generalisable to other biographic ethic groups. (3) the study is lack of anthropometric parameters regarding central obesity, lifestyle and dietary factors. In addition,we cannot determine the severity of NAFLD diagnosed by ultrasonographic. Due to the limitation of the original data, we cannot observe the correlation between insulin resistance and NAFLD, although insulin resistance may be closely associated with NAFLD in non-obese individuals[26].
Increased normal LDL-c levels are independently associated with an elevated risk of NAFLD in non-obese individuals. The relationship between normal LDL-c levels and the risk of NAFLD is non-linear. LDL-c is positively correlated with NAFLD when normal LDL-c levels is larger than 1.51 mmol/L.
Acknowledgments
We appreciate Sun D-Q, Wu S-J, Liu W-Y et al. for their spirit of scientific sharing. We also thank Chang-zhong Chen, Xin-Lin Chen and Chi-Chen of Yi-er college.
Authors’ contributions
XZG contributed to the drafting of the manuscript. LMY and ZJ analysed and interpreted the data. XBH contributed to the conception and critical revision of the manuscript, analysis and interpretation of the data and approved the final version of the submitted manuscript. All authors read and approved the final manuscript.
Funding
No funding was received.
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The research ethics were not required anymore, which had been authorized in the former study according to public policy statements of the dataset.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no conflict of interest.
Author details
1Department of Gastroenterology, LuHe Hospital of Capital Medical University, 82# Xinhuanan Road, Beijing 101149, China.