Article Text

Download PDFPDF

Association between non-alcoholic fatty liver disease-associated hepatic fibrosis and bone mineral density in postmenopausal women with type 2 diabetes or impaired glucose regulation
  1. Xiaopeng Zhu1,2,
  2. Hongmei Yan1,2,
  3. Xinxia Chang1,2,
  4. Mingfeng Xia1,2,
  5. Linshan Zhang1,2,
  6. Liu Wang1,2,
  7. Xiaoyang Sun1,2,
  8. Xinyu Yang1,2,
  9. Xin Gao1,2,
  10. Hua Bian1,2
  1. 1Department of Endocrinology and Metabolism, Zhongshan Hospital Fudan University, Shanghai, China
  2. 2Fudan Institute for Metabolic Disease, Fudan University, Shanghai, China
  1. Correspondence to Dr Hua Bian; zhongshan_bh{at}126.com; Dr Xin Gao; zhongshan_endo{at}126.com

Abstract

Introduction To evaluate the association of non-alcoholic fatty liver disease (NAFLD)-associated hepatic fibrosis with bone mineral density (BMD) in postmenopausal women with type 2 diabetes mellitus (T2DM) or impaired glucose regulation (IGR).

Research design and methods Two cohorts including 46 subjects with biopsy-proven NAFLD and 445 subjects with proton magnetic resonance spectrum-proven NAFLD were enrolled in this study. All subjects were postmenopausal women with T2DM or IGR. BMD at the lumbar spine L1–L4 and hip was measured using dual-energy X-ray absorptiometry. NAFLD fibrosis stage and NAFLD fibrosis score were used to evaluate the severity of liver fibrosis.

Results In subjects with liver biopsy-proven NAFLD, BMD (T-score, Z-score and BMD value) in the advanced fibrosis group were significantly lower than that in the non-advanced fibrosis group (p<0.05). Fibrosis stage was negatively associated with T-score, Z-score and BMD value after adjusting for age, body mass index (BMI) and fasting plasma glucose (FPG). Additionally, fibrosis stage was independently associated with T-score, Z-score and BMD value after adjusting for age, BMI and FPG. These results were validated in a large cohort of 445 subjects. Additionally, bone metabolism-associated factors, including calcium and phosphate, were associated with liver fibrosis, indicating that bone metabolism may play a critical role in the association between liver fibrosis and BMD. Mechanically, parathyroid hormone and biomarkers of bone formation (osteocalcin and procollagen type 1 N-terminal propeptide) and bone resorption (procollagen type I carboxy terminal peptide β special sequence) were increased in subjects with advanced liver fibrosis than in subjects without advanced liver fibrosis, indicating that liver fibrosis decreased BMD probably via increasing bone turnover.

Conclusions NAFLD-associated hepatic fibrosis was negatively associated with decreased BMD in postmenopausal women with T2DM or IGR. Liver fibrosis decreased BMD probably via increasing bone turnover. Severe liver fibrosis may represent high risk for osteoporosis in postmenopausal women with T2DM or IGR.

  • bone density
  • liver disease
  • biopsy
  • type 2 diabetes
http://creativecommons.org/licenses/by-nc/4.0/

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, appropriate credit is given, any changes made indicated, 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.

Significance of this study

What is already known about this subject?

  • There is an association of non-alcoholic fatty liver disease (NAFLD) with bone mineral density (BMD) in children, adult males and females; however, the relation between NAFLD-associated hepatic fibrosis and BMD remains unknown.

What are the new findings?

  • NAFLD-associated hepatic fibrosis was negatively associated with decreased BMD in postmenopausal women with type 2 diabetes mellitus (T2DM) or impaired glucose regulation (IGR).

  • Bone metabolism-associated factors, including calcium and phosphate, were associated with liver fibrosis.

  • Biomarkers of bone formation and resorption were increased in subjects with advanced liver fibrosis than in subjects without advanced liver fibrosis.

  • Liver fibrosis decreased BMD probably via increasing bone turnover.

How might these results change the focus of research or clinical practice?

  • Severe NAFLD-associated liver fibrosis may represent high risk for osteoporosis in postmenopausal women with T2DM or IGR; screening of reduced BMD in this population may be needed.

Introduction

Non-alcoholic fatty liver disease (NAFLD), characterized by excessive fat deposition in hepatocytes without overt alcohol consumption or other liver injury, is becoming a major public health issue worldwide.1 A recent study shows that the global prevalence of NAFLD is up to 25.2% and the proportion of fibrosis progression in non-alcoholic steatohepatitis (NASH) is 40.8%.2 NAFLD, the pathology spectrum of which ranges from simple steatosis to NASH and cirrhosis, is considered as a hepatic manifestation of metabolic syndrome.3 In addition, NAFLD is a multiple system disease, which is usually accompanied by extrahepatic complications, such as type 2 diabetes mellitus (T2DM), cardiovascular disease, chronic kidney disease, osteoporosis and neoplasm.4–7 It seems that NAFLD probably plays a significant role in the pathogenesis of other organs or systems.

Osteoporosis is defined as a bone mineral density (BMD) >2.5 standard deviations (SDs) below the young normal mean at the spine, hip or mid-radius by World Health Organization (WHO). Osteoporosis is also a significant public health issue worldwide like NAFLD, especially in elderly population.8 Previous studies illustrated that there was a complex network between adipose tissue, liver and bone.9 10 Consistent with these findings, several studies demonstrated there was an association of NAFLD with BMD in children, adult males and females.11–20 However, the relation between NAFLD and BMD remains controversial. Our previous study including 1659 subjects found that there was a negative correlation between liver fat content (LFC) and BMD in middle-aged and elderly Chinese men and postmenopausal women.21 However, the method of NAFLD diagnosis in these studies is ultrasonography,12 14–16 18–20 which is subjective and is greatly influenced by operator. Liver biopsy is still the gold standard to diagnose liver inflammation and fibrosis at present. Furthermore, whether NAFLD-associated hepatic fibrosis is correlated with BMD remains unknown. Previous study showed chronic inflammation and cirrhosis were risk factors for reduced BMD.22 Therefore, we hypothesized that NAFLD-associated hepatic fibrosis was correlated with low BMD. Previous studies showed that postmenopausal women had a high risk of reduced BMD, and patients with T2DM had a high risk of NAFLD and fracture.23–25 Therefore, we chose postmenopausal women with T2DM or impaired glucose regulation (IGR) as our target population.

In this study, we evaluated the association between NAFLD-associated liver fibrosis and BMD (T-score, Z-score and BMD value) in 46 postmenopausal women with biopsy-proven NAFLD and validated this correlation in a large cohort of 445 postmenopausal women with proton magnetic resonance spectrum (1H-MRS)-proven NAFLD. We also explored the role of bone metabolism-associated factors in NAFLD-associated liver fibrosis, in order to uncover the potential mechanism of low BMD in patients with NAFLD-associated liver fibrosis.

Methods

Study subjects

Two cohorts including 46 subjects with biopsy-proven NAFLD (mean age 58±5 years) and 445 subjects with 1H-MRS-proven NAFLD (mean age 64±9 years) from Zhongshan Hospital, Fudan University, were included in this study. All subjects had IGR (impaired fasting glucose or impaired glucose tolerance) or T2DM. The inclusion criteria for participants were as follows: i) postmenopausal women; ii) fatty liver diagnosed by biopsy or 1H-MRS; iii) no known acute or chronic disease other than obesity, hypertension, diabetes, osteoporosis or NAFLD, based on history and physical examination. Exclusion criteria included i) autoimmune, viral or drug-induced liver disease; ii) excessive use of alcohol (>20 g/day for men, >10 g/day in women26); iii) celiac disease.

Dual-energy X-ray absorptiometry measurements of BMD

T-score, Z-score and BMD value at the lumbar spine L1–L4 (lumbar total) and hip (femoral neck, inter, troch and total) were measured using dual-energy X-ray absorptiometry (DXA) (Hologic Discovery DXA system).

Evaluation of liver steatosis, inflammation and fibrosis in subjects with biopsy

Liver biopsy was performed to determine the severity of liver steatosis, inflammation and fibrosis. NASH was diagnosed by the joint presence of steatosis, ballooning and lobular inflammation.27 NAFLD activity score (NAS) and NAFLD fibrosis stage (S) were applied to evaluate the severity of NAFLD.28–30 NAS, ranging from 0 to 8, consists of steatosis (0–3), lobular inflammation (0–3) and ballooning (0–2). Briefly, steatosis was scored as 0 (<5%), 1 (5%–33%), 2 (34%–66%) and 3 (>66%); lobular inflammation (200× field) was scored as 0 (no foci), 1 (<2 foci), 2 (2–4 foci) and 3 (>4 foci); ballooning was scored as 0 (none), 1 (few balloon cells) and 2 (many cells/prominent ballooning).28 S comprises four stages, including S0 (none), S1 (perisinusoidal or periportal), S2 (perisinusoidal and portal/periportal), S3 (bridging fibrosis) and S4 (cirrhosis).28

Evaluation of liver steatosis, inflammation and fibrosis in subjects without biopsy

1H-MRS was used to quantify LFC. Serum alanine aminotransferase (ALT) level was applied as a simple marker of liver inflammation. NAFLD fibrosis score (NFS) was calculated to evaluate the severity of hepatic fibrosis.31 NFS=−1.675+0.037×age (years)+0.094×body mass index (BMI, kg/m2)+1.13×impaired fasting glucose/diabetes (yes=1, no=0)+0.99×aspartate aminotransferase (AST)/ALT ratio–0.013×platelet (109/L)–0.66×albumin (g/dL).31 NFS ≥0.676 was considered as advanced fibrosis.31

Anthropometric and laboratory measurements

All subjects performed routine anthropometric measurements, serum biochemical examinations and history collections. A standard meal glucose tolerance test was performed in patients with a history of T2DM; an oral glucose tolerance test (OGTT, 75 g glucose) was performed in patients without a history of T2DM. Standing height, waist circumference and bodyweight were measured without shoes and outer clothing, and the BMI was calculated as weight (kg) divided by standing height squared (m2). Blood samples were collected after a fasting period of at least 12 hours. Total serum cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL-ch), low-density lipoprotein cholesterol (LDL-ch), calcium, phosphate, uric acid, parathyroid hormone (PTH) and liver enzyme (ALT, AST, gamma-glutamyl transferase (γ-GT)) levels were measured using a model 7600 automated bio-analyzer (Hitachi, Tokyo, Japan). The fasting plasma glucose (FPG) and the 2-hour OGTT glucose following an OGTT or a standard meal were measured through glucose oxidase method. The fasting C-peptide and 2-hour OGTT C-peptide were measured via electrochemiluminescence. Hemoglobin A1c (HbA1c) was measured via high performance liquid chromatography. Plasma osteocalcin and 25-hydroxyvitamin D (25(OH)D) were measured using electrochemiluminescence immunoassay.

Statistical analysis

All statistical analyses were carried out using SPSS software V.20.0 (SPSS, Chicago, Illinois, USA). All data were texted for normality and homogeneity of variance. Normally distributed data are shown as the mean±SD, whereas non-normally distributed data are shown as median followed by the 25th and 75th percentiles. For normal and non-normal distribution of data, the unpaired Student’s t-test and the Mann-Whitney U test were used to compare mean values between two groups, respectively. Pearson’s correlation coefficient was calculated for normal distribution of data and Spearman’s rank correlation coefficient for non-normal distribution of data. Multiple linear regression analysis was performed to investigate the independent associated factors with BMD. All tests were two-sided, and p<0.05 was considered as statistical significance.

Results

Clinical characteristics of subjects with biopsy-proven NAFLD

A total of 46 subjects (mean age 58±5 years) were divided into non-advanced fibrosis and advanced fibrosis groups based on the fibrosis stage (S≥2 is considered as advanced fibrosis).28 The proportion of T2DM in the advanced fibrosis group was higher than that in the non-advanced fibrosis group (100.0% vs 78.9%, p<0.05, table 1). No significant difference was found in serum liver enzymes (AST and γ-GT), FPG, 2-hour postprandial plasma glucose (PG), HbA1c, TG, HDL-ch, LDL-ch, TC, osteocalcin, 25(OH)D and LFC between the two groups.

Table 1

Clinical characteristics of the subjects

BMD in different stage of liver fibrosis in subjects with biopsy-proven NAFLD

We next compared the BMD (T-score, Z-score and BMD value) in groups divided by fibrosis stage based on liver biopsy. In subjects with liver biopsy, only T-score (femoral inter) in the advanced fibrosis group was significantly lower than that in the non-advanced fibrosis group (p<0.05, table 2). In addition, subjects in the advanced fibrosis group had a trend toward lower T-score on hip at femoral neck, troch and total (table 2). After adjusting for age and BMI, T-score (femoral inter) and T-score (femoral total) in the advanced fibrosis group were significantly lower than those in the non-advanced fibrosis group (all p<0.05, table 2). The Z-score is a value that compares with the age-matched normal. Z-score (femoral total) in the advanced fibrosis group was significantly lower than that in the non-advanced fibrosis group after adjusting for age and BMI (p<0.05, table 2). We also analyzed the value of BMD in non-advanced fibrosis and advanced fibrosis groups. BMD value (femoral total) in the advanced fibrosis group was significantly lower than that in the non-advanced fibrosis group after adjusting for age and BMI (p<0.05, table 2).

Table 2

The difference of BMD in the different fibrosis groups in postmenopausal women

We further explored BMD in the non-NASH and NASH groups. No significant difference was found in T-score, Z-score and BMD value between the two groups (online supplementary table 1).

Supplemental material

Association of multiple metabolic factors with BMD or liver fibrosis in subjects with biopsy-proven NAFLD

The association of multiple metabolic factors with BMD (T-score, Z-score and BMD value) or liver fibrosis in postmenopausal women with biopsy-proven NAFLD was determined. Serum HbA1c, phosphate, osteocalcin, PTH, β-CTX and PINP levels were negatively correlated with T-score, Z-score and BMD value after adjusting for BMI, age and FPG (all p<0.05, online supplementary table 2, online supplementary table 3 and online supplementary table 4). In addition, fibrosis stage was negatively associated with T-score, Z-score and BMD value on hip at femoral inter and total after adjusting for BMI, age and FPG (figure 1A and B, online supplementary table 2, online supplementary table 3 and online supplementary table 4).

Supplemental material

Supplemental material

Supplemental material

Figure 1

Association between liver fibrosis and bone mental density in postmenopausal women. (A) The association between fibrosis stage and T-score (femoral inter). (B) The association between fibrosis stage and T-score (femoral total). (C) The association between NFS and T-score (femoral troch). NAFLD, non-alcoholic fatty liver disease; NFS, NAFLD fibrosis score.

The association between multiple metabolic factors relative to bone metabolism and liver fibrosis in postmenopausal women with liver biopsy was explored. However, no association between metabolic factors (except PTH) relative to bone metabolism and liver fibrosis was found in these subjects (table 3).

Table 3

The association between multiple metabolic parameters and liver fibrosis in postmenopausal women after adjusting for BMI, age and FPG

We further used multiple linear regression analysis to investigate the independent association between liver fibrosis and BMD (T-score, Z-score and BMD value) in postmenopausal women with liver biopsy. After adjusting for age, BMI and FPG, age and fibrosis stage were independently associated with T-score (femoral inter) (table 4); age, BMI and fibrosis stage were independently associated with T-score (femoral total) (table 4). Additionally, BMI, FPG and fibrosis stage were independently associated with Z-score (femoral inter) (online supplementary table 5); BMI and fibrosis stage were independently associated with Z-score (femoral total) after adjusting for age, BMI and FPG (online supplementary table 5). The linear regression was also performed for BMD value. Age, BMI and fibrosis stage were independently associated with BMD value on hip at femoral inter and total after adjusting for age, BMI and FPG (online supplementary table 6).

Supplemental material

Supplemental material

Table 4

Multiple linear regression analysis for the association between multiple metabolic parameters and BMD in postmenopausal women (T-score)

Clinical characteristics of subjects with 1H-MRS-proven NAFLD

We validated the results above in a large cohort of 445 postmenopausal women (mean age 64±9 years) with 1H-MRS-proven NAFLD. The subjects were divided into non-advanced fibrosis and advanced fibrosis groups based on the NFS (NFS ≥0.676 was considered as advanced fibrosis). Age, BMI, waist circumference, and serum PTH level in the advanced fibrosis group were higher than those in the non-advanced fibrosis group (all p<0.01, table 1), whereas serum total cholesterol, calcium, and phosphate in the advanced fibrosis group were lower than those in the non-advanced fibrosis group (all p<0.01, table 1). No significant difference was found in serum FPG, 2-hour PG, HbA1c, TG, HDL-ch, LDL-ch, osteocalcin, 25(OH)D, LFC and the proportion of diabetes between the two groups.

BMD in the different stage of liver fibrosis in subjects with 1H-MRS-proven NAFLD

In agreement with the result from subjects with liver biopsy, T-score on hip at femoral neck, troch, inter and total in the advanced fibrosis group were significantly lower than those in the non-advanced fibrosis group (all p<0.001, table 2). No significant difference was found in T-score (lumbar total) between the two groups. After adjusting for age and BMI, T-score on hip at femoral troch and inter in the advanced fibrosis group were still significantly lower than those in the non-advanced fibrosis group (all p<0.05, table 2). Additionally, Z-score on hip at femoral troch and total in the advanced fibrosis group were significantly lower than those in the non-advanced fibrosis group after adjusting for age and BMI (all p<0.05, table 2). BMD value on hip at femoral troch, inter and total in the advanced fibrosis group were significantly lower than those in the non-advanced fibrosis group after adjusting for age and BMI (all p<0.05, table 2).

Association of multiple metabolic factors with BMD or liver fibrosis in subjects with 1H-MRS-proven NAFLD

The association of multiple metabolic factors with BMD (T-score, Z-score and BMD value) or liver fibrosis in postmenopausal women with 1H-MRS-proven NAFLD was determined. NFS was negatively associated with T-score (femoral troch), Z-score (femoral troch), BMD value (femoral troch) and BMD value (femoral total) after adjusting for BMI, age and FPG (figure 1C, online supplementary table 2, online supplementary table 3 and online supplementary table 4). Additionally, serum osteocalcin, PTH, β-CTX, PINP and HDL-ch were also negatively associated with T-score, Z-score and BMD value, whereas serum ALT was positively associated with T-score and Z-score (all p<0.05, online supplementary table 2, online supplementary table 3 and online supplementary table 4) after adjusting for BMI, age and FPG.

We also explored the association between multiple metabolic factors relative to bone metabolism and liver fibrosis in this cohort. Serum calcium, phosphate, magnesium and 25(OH)D levels were negatively associated with NFS after adjusting for BMI, age and FPG (table 3).

We further performed multiple linear regression analysis to investigate the independent association between liver fibrosis and BMD in subjects with1H-MRS-proven NAFLD. In line with previous results, age, BMI and NFS were independently associated with T-score on hip at femoral troch and total (table 4) after adjusting for age, BMI, FPG, TC, HDL-ch and LDL-ch. Additionally, BMI and NFS were independently associated with Z-score on hip at femoral troch and total; BMI, FPG and NFS were independently associated with Z-score (femoral inter) after adjusting for age, BMI, FPG, TC, HDL-ch and LDL-ch (online supplementary table 5). Age, BMI and NFS were independently associated with BMD value on hip at femoral inter, troch and total after adjusting for age, BMI, FPG, TC, HDL-ch and LDL-ch (online supplementary table 6).

Discussion

In this study, based on liver biopsy, we found that NAFLD-associated hepatic fibrosis had an independent association with low BMD (T-score, Z-score and BMD value) in postmenopausal women with T2DM or IGR. The results were verified in a large cohort of postmenopausal women with T2DM or IGR, showing that NAFLD-associated liver fibrosis evaluated by NFS had an independent association with low BMD (T-score, Z-score and BMD value) in postmenopausal women with T2DM or IGR. Additionally, bone metabolism-associated factors, including calcium and phosphate, were associated with liver fibrosis, indicating that bone metabolism probably plays a critical role in the association between liver fibrosis and BMD. Based on these results, we tried to investigate the potential mechanism underlying the relation between BMD and liver fibrosis. We found that PTH and biomarkers of bone formation (osteocalcin and PINP) and bone resorption (β-CTX) were increased in subjects with advanced liver fibrosis compared with those in subjects without advanced liver fibrosis, indicating that liver fibrosis decreased BMD probably via increasing bone turnover. These data suggest that severe liver fibrosis may represent high risk for osteoporosis in postmenopausal women with T2DM or IGR.

Few studies had focused on the relation of liver inflammation and fibrosis with BMD based on liver histology. Recent study including 35 children with histologically proven NAFLD has showed that children with NASH had a lower lumbar Z-score as well as a lower whole body Z-score, compared with those with simple steatosis.17 Another study enrolling 38 children with obesity with biopsy-proven NAFLD and 38 children with obesity without NAFLD demonstrated that children with obesity with NAFLD had lower Z-score than subjects without NAFLD, and children with NASH had a lower Z-score than subjects with NAFLD but without NASH.13 These two studies were conducted in children and had not reported the association between liver fibrosis and BMD. Recent study has revealed significant liver fibrosis diagnosed via liver transient elastography was independently associated with decreased BMD in patients with NAFLD.32 However, liver transient elastography is not the gold standard for evaluation of liver fibrosis. In current study, we found that advanced liver fibrosis diagnosed by the gold standard, liver biopsy, was independently associated with decreased BMD in postmenopausal women with NAFLD with T2DM or IGR.

Due to the small size (46 subjects) of postmenopausal women with liver biopsy, we validated the results in another large cohort of 445 postmenopausal women with 1H-MRS-proven NAFLD. As expected, advanced liver fibrosis was also independently associated with decreased BMD in this large cohort.

Although there is an association between NAFLD fibrosis and decreased BMD, the potential molecular mechanisms underlying this relation are still unknown. Study showed that NAFLD, especially NASH, may participate in the pathophysiology of decreased BMD via the systemic release of multiple pro-inflammatory, procoagulant, pro-oxidant and profibrogenic mediators and/or via the direct effect on hepatic and systemic insulin resistance.33 In present study, 25(OH)D, calcium, phosphate and magnesium were negatively associated with NFS. This might be because the synthesis of 25(OH)D was reduced in liver fibrosis condition. Reduced 25(OH)D led to decreased calcium absorption, which resulted in decreased BMD. As to phosphate, the explanation was that LFC was decreased as liver fibrosis became more serious, and phosphate was positively associated with LFC.22 Therefore, phosphate was negatively associated with liver fibrosis.

Osteocalcin, a major non-collagenous hydroxyapatite-binding protein produced by osteoblasts, which is considered as a biomarker of bone formation,34 35 was reported to be negative associated with BMD both in premenopausal women with obesity and in postmenopausal women.36 37 PINP is also a biomarker of bone formation. In current study, we found that osteocalcin and PINP were positively associated with NFS. β-CTX, a biomarker of bone resorption, had a trend toward positively associated with fibrosis stage/NFS in this study. Therefore, we speculated that liver fibrosis increased bone turnover, thus leading to decreased BMD. Moreover, in this study, increased PTH could be explained by compensation for low calcium, and increased PTH indirectly supported the increased bone turnover.

Previous studies showed that age was negatively associated with BMD, and BMI was positively associated with BMD.38 39 We also found the same results in the present study. Study showed that patients with diabetes were prone to develop osteoporosis.40 In current study, we adjusted for FPG and found that NAFLD fibrosis was still an independent risk factor for decreased BMD.

Several limitations exist in this study. First, this is a cross-sectional study, and the association of NAFLD fibrosis with BMD should be further tested and verified in a prospective study. Second, the sample size of subjects with liver biopsy is small, the results of which are limited. Third, most subjects in this study have NAFLD accompanied by T2DM or IGR. Although we adjusted for FPG, the conclusion should be further investigated in non-diabetes population.

In conclusion, NAFLD-associated hepatic fibrosis was independently associated with decreased BMD in postmenopausal women with T2DM or IGR. Severe liver fibrosis may represent high risk for osteoporosis in postmenopausal women with T2DM or IGR.

Acknowledgments

The authors would like to thank Zhongshan Hospital for supporting this study.

References

Footnotes

  • Contributors HB and XG designed the study. XZ, LZ, LW, XS and XY collected the data. HY, XC and MX assessed the data. XZ and HB performed statistical analyses, and wrote the manuscript. XG edited the manuscript. All authors have approved the final version of the article.

  • Funding This work was supported by the financial support of National Natural Science Foundation of China (81471073 to HB); National Key R&D Program of China (2017YFC1309800 and 2017YFC1309801 to HB); Shanghai Municipal Population and Family Planning Commission (201740092 to HY); Special Project of Integrating Traditional Chinese and Western Medicine in Shanghai General Hospital from Shanghai Municipal Population and Family Planning Commission and Shanghai TCM Development Office (ZY(2018-2020)-FWTX-3019 to HY).

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval This study (approval number: B2013-132) was in accordance with the Helsinki Declaration of 1975 and was approved by the ethics committee of Zhongshan Hospital, Fudan University, and each subject provided written informed consent before taking part.

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

  • Data availability statement All data relevant to the study are included in the article or uploaded as supplementary information.