Article Text

Download PDFPDF

Fatty liver index, albuminuria and the association with chronic kidney disease: a population-based study in China
  1. Kan Sun1,
  2. Diaozhu Lin1,
  3. Feng Li1,
  4. Yiqin Qi1,
  5. Wanting Feng1,
  6. Li Yan1,
  7. Chaogang Chen2,
  8. Meng Ren1,
  9. Dan Liu1
  1. 1 Department of Endocrinology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, Republic of China
  2. 2 Department of Nutrition, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, Republic of China
  1. Correspondence to Professor Dan Liu; liudanye66{at}126.com

Abstract

Objectives The effects of lipid metabolism disorder on renal damage have drawn much attention. Using the fatty liver index (FLI) as a validated indicator of hepatic steatosis, this study aims to provide insight about the possible links between fatty liver and the development of chronic kidney disease (CKD).

Setting Hospital.

Participants We performed a population-based study on 9436 subjects aged 40 years or older.

Primary and secondary outcome measures FLI is calculated using an algorithm based on body mass index, waist circumference, triglycerides and γ-glutamyltransferase. Increased urinary albumin excretion was defined according to the urinary albumin to creatinine ratio ranges ≥30 mg/g. CKD was defined as estimated glomerular filtration rate (eGFR) <60 mL/min/1.73 m² or presence of albuminuria.

Results There were 620 (6.6%) subjects categorised to have increased urinary albumin excretion and 753 (8.0%) subjects categorised to have CKD. Participants with higher FLI had increased age, blood pressure, low-density lipoprotein cholesterol, fasting plasma glucose, fasting insulin and decreased eGFR level. Prevalence of increased urinary albumin excretion and CKD tended to increase with the elevated FLI quartiles. In logistic regression analysis, compared with subjects in the lowest quartile of FLI, the adjusted ORs in the highest quartile were 2.30 (95% CI 1.36 to 3.90) for increased urinary albumin excretion and 1.93 (95% CI 1.18 to 3.15) for CKD.

Conclusion Hepatic steatosis evaluated by FLI is independently associated with increased urinary albumin excretion and prevalence of CKD in middle-aged and elderly Chinese.

  • fatty liver index
  • hepatic steatosis
  • increased urinary albumin excretion
  • chronic kidney disease

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • The study was performed in a large population-based cohort of 9436 Chinese subjects.

  • Findings of the study may be applied to the majority of patients in general practice with suspected hepatic steatosis.

  • Results should be interpreted cautiously due to the observational design of the current study.

Introduction 

Chronic kidney disease (CKD) has become one of the leading public health problems worldwide.1 A recent national survey conducted between 2007 and 2010 reports that the prevalence of CKD was 10.8%, representing an estimated 119.5 million patients in China with chronic kidney damage.2 In addition to CKD, an increasing number of studies have provided substantial evidence of albuminuria as a risk factor for future cardiovascular events.3 Both renal and cardiovascular diseases sharing similar traditional risk factors, such as lipid metabolism disorder, could have particularly broad implications for the outcome of cardiovascular morbidity and mortality.

Association of hepatic steatosis with CKD development and its impact on the reduction of estimated glomerular filtration rate (eGFR) have been extensively investigated over the past decade.4 The substantial evidence linked hepatic steatosis to the increased risk and severity of CKD, which may be a target in the prevention and treatment of the disease.5 As a convenient scoring system for the presence of hepatic lipid deposits, the fatty liver index (FLI) is a surrogate steatosis biomarker developed in a cohort of patients from the general population.6 Compared with other techniques for evaluating hepatic steatosis, FLI is simple to obtain as body mass index (BMI), waist circumference (WC), triglycerides (TG) and γ-glutamyltransferase (γ-GGT) are routine measurements in clinical practice. Previous studies have demonstrated that FLI could determine fatty liver disease, incident type 2 diabetes and incident hypertension with considerable accuracy.6–8 Moreover, FLI is associated with insulin resistance, early atherosclerosis and risk of coronary heart disease, which could help physicians early detect subjects of greater cardiovascular risk and select patients for intensified lifestyle counselling.9 10

Clarifying the association of FLI with albuminuria and prevalent CKD would probably shed light on the prevention and pre-emptive treatment of related diseases. Recently, a cross-sectional study was conducted to investigate the association between FLI and CKD by recruiting adults undergoing a health check-up.11 However, by including only 731 subjects, the study did not evaluate the association between FLI and albuminuria, either. Therefore, we analysed data from a community-based Chinese population to comprehensively look into the relationship of FLI with both increased urinary albumin excretion and CKD.

Subjects and methods

Study population and design

We performed a cross-sectional study in a community in Guangzhou, China, from June to November 2011. The study population was from the Risk Evaluation of cAncers in Chinese diabeTic Individuals: a lONgitudinal (REACTION) study and the details of this study have been published previously.12–14 During the recruitment phase, a total of 10 104 residents aged 40 years or older were invited to participate by examination notices or home visits. In total, 9916 subjects signed the consent form and agreed to participate in the survey. The participation rate was 98.1%. Subjects who failed to provide information (BMI, n=206; WC, n=62; TG, n=23; γ-GGT, n=38; or urinary albumin to creatinine ratio (ACR), n=149) were excluded from the analyses. Accordingly, a total of 9438 eligible individuals were included in the final data analyses. Written informed consent was obtained from each participant before data collection.

Clinical and biochemical measurements

We collected information on lifestyle factors, sociodemographic characteristics and family history using a standard questionnaire. Smoking or drinking habit was classified as ‘never’, ‘current’ (smoking or drinking regularly in the past 6 months) or ‘ever’ (cessation of smoking or drinking more than 6 months).15 A short form of the International Physical Activity Questionnaire was used to estimate physical activity at leisure time by adding questions on frequency and duration of moderate or vigorous activities and walking.16 Separate metabolic equivalent hours per week (MET-hour/week) were calculated to evaluate total physical activity.

All participants completed the anthropometrical measurements with the assistance of trained staff using standard protocols. Blood pressure measurements were obtained three times consecutively by the same observer at a 5 min interval using an automated electronic device (OMRON, Omron Company, Dalian, China). The average of three blood pressure measurements was used for analysis. Body height and body weight were recorded to the nearest 0.1 cm and 0.1 kg while participants were wearing light indoor clothing without shoes. BMI was calculated as weight in kilograms divided by height in metres squared (kg/m2). Obesity was defined as BMI ≥28, and overweight was defined as BMI ≥24 and <28.17 WC was measured at the umbilical level with the participant in standing position, at the end of gentle expiration.

Venous blood samples were collected for laboratory tests after an overnight fasting of at least 10 hours. Measurement of fasting plasma glucose (FPG), fasting serum insulin, TG, total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), creatinine, γ-GGT, aspartate aminotransferase (AST) and alanine aminotransferase (ALT) was done using an autoanalyser (Beckman CX-7 Biochemical Autoanalyzer, Brea, California, USA).

As surrogate marker of hepatic steatosis, FLI was analysed based on BMI, WC, TG and γ-GGT, and has been validated against liver ultrasound in the general population and proven accurate in detecting fatty liver.6 10 FLI is calculated as the following: FLI = (e0.953 × loge(TG) + 0.139 × BMI + 0.718 × loge(GGT) + 0.053 × WC − 15.745) / (1 + e0.953 × loge (TG) + 0.139 × BMI + 0.718 × loge (GGT) + 0.053 × WC − 15.745) × 100. The abbreviated Modification of Diet in Renal Disease formula recalibrated for Chinese population was used to calculate eGFR expressed in mL/min/1.73 m2 using the following formula: eGFR=175 × (serum creatinine × 0.011)−1.234 × (age)−0.179 × (0.79 if female), where serum creatinine was expressed in μmol/L.18 Diabetes was diagnosed according to the 1999 WHO diagnostic criteria.19

Definition of increased urinary albumin excretion, CKD and NAFLD

Definitions of abnormalities in albumin excretion were according to the latest guidelines of the American Diabetes Association’s Standards of Medical Care.20 The first morning spot urine samples were collected to assess the ACR. Urine albumin and creatinine were measured by chemiluminescence immunoassay (Siemens Immulite 2000, Orlando, Florida, USA) and the Jaffe’s kinetic method (Biobase-Crystal, Jinan, China) on the automatic analyser, respectively. ACR was calculated by dividing the urinary albumin concentrations by the urinary creatinine concentrations, and is expressed in mg/g. The primary and secondary outcome measures were increased urinary albumin excretion and CKD, respectively. Increased urinary albumin excretion was defined according to the ACR ranges ≥30 mg/g. CKD was defined as eGFR <60 mL/min/1.73 m² or presence of albuminuria (ACR ≥30 mg/g). The optimal cut-off value of FLI in predicting non-alcoholic fatty liver disease (NAFLD) was 30 in Asian populations.21 Therefore, we classified the study population in non-current drinking group into NAFLD group (FLI ≥30) and non-NAFLD group (FLI <30).

Statistical analysis

Statistical analysis was performed using SAS V.9.2. Continuous variables were presented as mean±SD except for skewed variables, which were presented as median (IQR). Categorical variables were expressed as numbers (proportions). FLI, FPG, TG, ALT, AST, γ-GGT and MET-hour/week were logarithmically transformed before analysis due to a non-normal distribution. FLI was presented as quartiles, and linear regression analysis was used to test for trend across groups. Differences among groups were tested by one-way analysis of variance, and post hoc comparisons were performed using Bonferroni correction. Comparisons between categorical variables were performed using the χ2 test.

Pearson’s correlations were performed to test the correlations between FLI and the risk factors for kidney disease. Variables significant at P<0.20 in Pearson’s correlations were put into the multivariate stepwise linear regression models to identify factors that independently associated with FLI. We analysed the impact of FLI on the prevalence of increased urinary albumin excretion and CKD. Unadjusted and multivariate-adjusted logistic regression analyses were used to assess the risk of prevalent increased urinary albumin excretion and CKD in relation to each quartile increase in FLI level. Variables considered as potential covariates and significant in the stepwise linear regression were put into multivariate-adjusted logistic regression analysis. Model 1 is unadjusted. Model 2 is adjusted for age. Model 3 is adjusted for age, sex, current smoking status, current drinking status, physical activity, systolic blood pressure (SBP), diastolic blood pressure (DBP), LDL-C, fasting insulin, ALT and AST. Model 4 is adjusted for age, sex, BMI, WC, current smoking status, current drinking status, physical activity, SBP, DBP, TG, LDL-C, fasting insulin, ALT, AST and γ-GGT. OR and the corresponding 95% CI were calculated. The relationship of FLI level with albuminuria and CKD was also explored in subgroups stratified by gender (men/women), age (≥60/<60 years), degree of obesity (normal/overweight/obesity), current smoking (yes/no), current drinking (yes/no), hypertension (yes/no) and diabetes (yes/no). Tests for interaction were performed by including simultaneously each strata factor, the quartiles of FLI level and the respective interaction terms (strata factor multiplied by quartiles of FLI level) in the models.

All statistical tests were two-sided, and a P value <0.05 was considered statistically significant.

Results

Clinical characteristics of the study population

Among the 9436 enrolled individuals, the mean age was 55.9±8.0 years. The median FLI was 19.1 with IQR of 8.6–37.4. There were 620 (6.6%) subjects categorised to have increased urinary albumin excretion and 753 (8.0%) subjects with CKD, respectively. Table 1 shows the clinical and biochemical characteristics of the participants according to FLI quartiles. Participants with higher FLI level had increased age, BMI, WC, SBP, DBP, TG, TC, LDL-C, FPG, fasting insulin, ALT, AST and γ-GGT, and higher proportions of current smokers and current drinkers (all P for trend <0.0001). Higher FLI level also associated with decreased HDL-C and eGFR (all P for trend <0.0001).

Table 1

Characteristics of study population by FLI quartiles

Associations between FLI and metabolic risk factors

Analysis of Pearson’s correlation showed that age, sex, BMI, WC, SBP, DBP, TG, TC, HDL-C, LDL-C, FPG, fasting insulin, ALT, AST, γ-GGT and eGFR were significantly correlated with FLI level. Further multivariate stepwise linear regression showed that age, sex, BMI, WC, SBP, DBP, TG, LDL-C, fasting insulin, ALT, AST and γ-GGT were independent determinants of FLI level (table 2).

Table 2

Pearson’s correlation and stepwise regression analysis of determinants of FLI

Associations of FLI with increased urinary albumin excretion and CKD

As shown in figure 1A, from the lowest quartile to the highest quartile of FLI level, the prevalence of increased urinary albumin excretion was 3.64%, 4.83%, 6.23% and 11.57%, respectively (P for trend <0.0001). Strikingly, the prevalence of CKD also tended to increase with the elevated FLI quartile (figure 1B, P for trend <0.0001). As shown in table 3, compared with participants in quartile 1 of FLI, univariate logistic regression analysis showed that participants in quartile 2, quartile 3 and quartile 4, respectively, have a significant correlation with increased odds of increased urinary albumin excretion and CKD (all P for trend <0.0001). In multivariate logistic regression analyses (model 3), the ORs of increased urinary albumin excretion for increasing FLI quartiles were 1.00 (reference), 0.96 (95% CI 0.66 to 1.39), 1.17 (95% CI 0.77 to 1.77) and 2.30 (95% CI 1.36 to 3.90). Similarly, the ORs of CKD for increasing FLI quartiles in model 3 were 1.00 (reference), 1.00 (95% CI 0.71 to 1.40), 1.03 (95% CI 0.70 to 1.51) and 1.93 (95% CI 1.18 to 3.15), respectively (table 3). The prevalence of increased urinary albumin excretion was 51.6% and 29.6% in FLI-established NAFLD and non-NAFLD group (P<0.0001). Similar trends were detected in the prevalence of CKD (NAFLD group: 49.9%; non-NAFLD group: 31.5%; P<0.0001). Compared with the participants in the non-NAFLD group, those in the NAFLD group had higher prevalence of increased urinary albumin excretion (OR 1.58, 95 % CI 1.18 to 2.13) and CKD (OR 1.39, 95 % CI 1.05 to 1.82) in multivariate logistic regression analyses.

Figure 1

Prevalence of increased urinary albumin excretion and chronic kidney disease (CKD) in different quartiles of fatty liver index levels: (A) increased urinary albumin excretion and (B) CKD.

Table 3

Risk of prevalent albuminuria and CKD according to quartiles of FLI

Subgroup analysis of FLI with increased urinary albumin excretion and CKD

As shown in figures 2 and 3, the associations of FLI level with increased urinary albumin excretion and CKD were not consistently the same in the subgroup analyses. A significant relationship of FLI level with both increased urinary albumin excretion and CKD was detected in women, younger subjects (age less than 60 years), overweight subjects, non-current smokers, non-current drinkers, and in those with hypertension or diabetes (all P<0.05). In the subgroup analysis, no statistically significant interaction term between quartiles of FLI and each strata factor was detected.

Figure 2

Risk of prevalent increased urinary albumin excretion with each quartile increase of fatty liver index levels in different subgroups. BMI, body mass index.

Figure 3

Risk of prevalent chronic kidney disease with each quartile increase of fatty liver index levels in different subgroups. BMI, body mass index.

Discussion

We evaluated the association between hepatic steatosis and kidney disease in a large population of middle-aged Chinese subjects from the REACTION study. The presence of fatty liver assessed by FLI was associated with increased urinary albumin excretion and reduction of eGFR in the present study. The association was independent of potential confounding risk factors. To our current knowledge, this is the largest population-based study to explore the association of FLI with both albuminuria and CKD in Asian population. Early intervention is of great importance for albuminuria and CKD, and the present findings may just give insights into lipid metabolism for the prevention and early detection of the diseases.

The problem of obesity and NAFLD is now increasingly recognised in the Asian population. Prevalence of obesity was 7.9% (8.4% in men and 7.6% in women) in southern China, which has increased dramatically over the past several decades.22 There is a strong correlation between established obesity and incidence of NAFLD. Pooled prevalence of NAFLD diagnosed by ultrasound, CT scan and MRI was estimated to be 27.4% in subjects from Asian countries aged over 30 years.23 Even among the non-obese Chinese, 8.9% developed NAFLD in 5 years from 2006 to 2011.24 Therefore, early and accurate diagnosis of NAFLD is of great importance. The best method for an accurate assessment and diagnosis of hepatic steatosis is histological analysis of biopsies.25 However, it is uneconomical to conduct liver biopsies especially because of our large sample population. Hepatic ultrasonic examination is widely used in clinical practice and epidemiological studies in detecting fatty infiltration of the liver.26 27 However, the non-invasive technique is not sensitive enough to detect mild steatosis and does not allow precise quantification of severity of steatosis in hepatic tissue.28

As another surrogate marker of histological fatty liver, FLI is defined as the accumulation of excessive liver fat.29 Based on previous research, FLI has been proven accurate in detecting fatty liver against liver ultrasound and in demonstrating the presence of hepatic fat against magnetic resonance spectroscopy.6 9 10 21 The superiority of this non-invasive assessment techniques is that a higher score will indicate a higher degree of steatosis in hepatic tissue. However, an optimal cut-off point of the FLI in evaluating liver fatty infiltration should be considered as it varied according to the study population.21 30 Originally, FLI>60 was suggested to rule in NAFLD in Caucasian subjects. However, the optimal cut-off value of FLI in predicting NAFLD was different in Asian populations. In one recent study, Huang et al 21 found that FLI could accurately identify NAFLD and the optimal cut-off point was 30 in middle-aged and elderly Chinese. FLI could also accurately identify ultrasonography fatty liver in a large-scale population in Taiwan but with different optimal cut-off values, while an FLI>35 for men and >20 for women rule in NAFLD in their study.30 Based on the results of our research on Chinese subjects, further studies are needed to externally discuss the optimal cut-off point of the FLI in predicting hepatic steatosis.

Detection and prevention of kidney disease progression and urinary albumin excretion are difficult to process in the early stage. Dyslipidaemia is increasingly recognised as an important pathogenic mechanism in the deterioration of renal function. Recently, we conducted a clinical investigation to assess the associations of routine lipid measures with kidney disease in the same cohort. In the study, discordant associations of lipid parameters with renal insufficiency were detected, while the TG to HDL-C ratio is a better marker for evaluating increased urinary albumin excretion and CKD.31 As one of the phenotypes of dyslipidaemia, the pathogenesis of hepatic steatosis is closely related to kidney disease with regard to insulin resistance and chronic inflammation.32 Based on recent studies, hepatokines, which are proteins secreted by hepatocytes, have been found to be linked to the induction of metabolic phenotypes through interorgan communication.33 Because of the high prevalence and burden of the fatty liver disease, it is important to identify which patients are most likely to be exposed to early-stage renal injury.23 Consequently, we closely monitor the association of hepatic steatosis predicted by FLI with prevalent increased urinary albumin excretion and CKD.

Consistent with our findings, a previous study reported that hepatic steatosis evaluated by FLI might contribute to CKD development.11 Elevated albuminuria is well known to be associated with increased risk for early diabetes renal damage; however, the identification and classification of kidney disease were assessed only by eGFR in that study. Moreover, 731 adults who underwent routine health evaluations were included in that study and the small sample size cannot better represent the whole population. By totally including 9438 subjects and adopting both albuminuria and eGFR for renal damage assessment, data in our study demonstrated that the FLI is associated with kidney disease, which might be an efficient screening indicator for the early prevention of related diseases in Chinese subjects. Recently, an interesting study by Giorda et al 34 reported that NAFLD is a dynamic condition in subjects with type 2 diabetes, and about 5% Italian patients with diabetes entering or leaving FLI assessed NAFLD status every year. They found that male sex and established organ damage, especially kidney function, were independent risk predictors for the dynamic NAFLD condition in a longitudinal 3-year analysis. Because of the similarity in traditional risk factors for both NAFLD and CKD, the relationship between the prevalence of earlier stages of kidney damage and the incidence of NAFLD is complex. Longitudinal observation of our cohort needs to be carried out to determine whether such dynamic condition existed in Chinese, especially in those with type 2 diabetes.

Alcohol consumption can profoundly disturb lipid metabolism, which has prominent effects on hepatic tissue steatosis and insulin sensitivity.35 However, potential health effects regarding alcohol consumption in this field is also worth attaching attention. A meta-analysis of intervention studies by Schrieks et al 36 showed that moderate alcohol intake could improve insulin sensitivity by decreasing fasting insulin level in women. Recently, a prospective cohort study found that alcohol consumption was consistently inversely associated with urinary albumin excretion and the risk of developing CKD.37 Therefore, advice on alcohol consumption to subjects with low-grade hepatic tissue steatosis should consider the full range of benefits and risks, especially among those who drink moderately.

Some limitations of the study must be noted. First, owing to the observational design of the current study, we should cautiously interpret the present findings as no causal inference can be drawn. Further prospective studies are therefore needed to determine the precise relationship between FLI and risk of renal diseases. Second, by including only Chinese subjects, the results of the present study might not be representative of other ethnic groups, especially those in the developed or undeveloped countries. To some extent, however, the present study of Chinese population was still a convenience sample and selection bias is inevitable. Third, when evaluating the findings of the present study, the results should be interpreted cautiously due to possible bias from using the indirect indicator FLI to assess fatty liver disease. The calculated FLI may relate to various liver diseases with associated steatosis and not only NAFLD, despite the fact that metabolic disturbances make obesity-related steatosis likely. The internal accuracy of FLI in evaluating hepatic steatosis should also be validated using other techniques before it can be employed for these purposes. Fourth, we observed that FLI seems to play a different efficiency for kidney disease assessment in different stratifications. A significant association of FLI with increased urinary albumin excretion and CKD was only detected in subjects without current alcohol consumption. Average daily alcohol intake influences FLI, and missing such data in the present study does not permit comparisons between and within alcoholic and non-alcoholic fatty liver disease groups. To better discriminate alcoholic fatty liver disease and NAFLD, further studies need to clearly describe the precise exposure of alcohol use by collecting histories of alcohol intake in a quantitative manner. Fifth, viral hepatitis infection is one of the most serious infectious diseases worldwide, and can be associated with both liver and kidney diseases. Recent survey data showed that the hepatitis B surface antigen and anti-hepatitis C virus-positive rates were already 6.1% and 3.0% in China. Epidemiology of viral hepatitis infection by hepatitis B virus and hepatitis C virus serological testing, therefore, should also be evaluated to strengthen the findings of the present study.38 Sixth, although a spectrum of covariates was included in the adjustment, other potential mediators such as daily energy and protein intake and medicine that influence the renin-angiotensin system of the subjects should also be considered in the present study.

In conclusion, by including a large population-based cohort, the present study provides evidence that increased FLI is independently associated with prevalence of albuminuria and CKD. Findings of the present study suggested we should pay more attention to albuminuria and eGFR variation in patients with dyslipidaemia and fatty liver disease. Further prospective studies are necessary to verify our findings in external populations.

Acknowledgments

We are indebted to the participants of the present study for their persistent outstanding support and to our colleagues for their valuable assistance.

References

  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.

Footnotes

  • KS and DL contributed equally.

  • Contributors Conceived and designed the experiments: LY and KS. Performed the experiments: FL, YQ, WF, CC, KS and DL. Analysed the data: KS and MR. Wrote the manuscript: KS and DL. All authors believe that the manuscript represents valid work and have reviewed and approved the final version. The work has not been published previously and is not under consideration for publication elsewhere, in part or in whole.

  • Funding This work was supported by grants from (1) the National Natural Science Foundation of China (81600642, 81471034, 81370910, 81300675); (2) the Natural Science Foundation of Guangdong Province, China (2015A030310433, 2017A030313831); (3) Sun Yat-sen University Medical 2016 Youth Teacher Research Funding Project (16ykpy27); (4) major project of the people’s livelihood science and technology in Guangzhou (201300000102); (5) the 863 Project of Young Scientist (SS2015AA020927); (6) the Zhujiang Star of Science and Technology Foundation in Guangzhou (2014J2200046); (7) grants from the Chinese Society of Endocrinology and National Clinical Research Center for Metabolic Diseases; (8) State Key Clinical Specialty Construction Project (2011); and (9) the Science and Technology Planning Project of Guangdong Province, China (2014A020212161).

  • Disclaimer The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

  • Competing interests None declared.

  • Patient consent Obtained.

  • Ethics approval The study protocol was approved by the Institutional Review Board of the Sun Yat-sen Memorial Hospital affiliated to Sun Yat-sen University and was in accordance with the principle of the Declaration of Helsinki II.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement Main document data and additional unpublished data from the study are available by sending email to lizyhenu@163.com with proper purposes.