Nomogram for the prediction of diabetic nephropathy risk among patients with type 2 diabetes mellitus based on a questionnaire and biochemical indicators: a retrospective study

Purpose: Develop a diabetic nephropathy incidence risk nomogram in a Chinese population with type 2 diabetes mellitus. Results: Predictors included systolic blood pressure, diastolic blood pressure, fasting blood glucose, glycosylated hemoglobin A1c, total triglycerides, serum creatinine, blood urea nitrogen and body mass index. The model displayed medium predictive power with a C-index of 0.744 and an area under curve of 0.744. Internal verification of C-index reached 0.737. The decision curve analysis showed the risk threshold was 20%. The value of net reclassification improvement and integrated discrimination improvement were 0.131, 0.05, and that the nomogram could be applied in clinical practice. Conclusion: Diabetic nephropathy incidence risk nomogram incorporating 8 features is useful to predict diabetic nephropathy incidence risk in type 2 diabetes mellitus patients. Methods: Questionnaires, physical examinations and biochemical tests were performed on 3489 T2DM patients in six communities in Shanghai. LASSO regression was used to optimize feature selection by running cyclic coordinate descent. Logistic regression analysis was applied to build a prediction model incorporating the selected features. The C-index, calibration plot, curve analysis, forest plot, net reclassification improvement, integrated discrimination improvement and internal validation were used to validate the discrimination, calibration and clinical usefulness of the model.


INTRODUCTION
Over the past 20 years, due to the increase in the obesity rate and the prevalence of sedentary lifestyles, the number of people worldwide diagnosed with diabetes has been increasing [1], which is rapidly becoming a public health problem in both developed and developing countries [2]. According to the report released by the International Diabetes Federation in 2015, the prevalence of diabetes among adults worldwide is 9.1% [3]. According to this figure, 415 million adults worldwide have diabetes, and 318 million adults have impaired blood glucose regulation, composing a group with a high risk of diabetes [4]. It was estimated that the number of diabetes patients among adults aged 20 and over will increase from 171 million in 2000 to 552 million in 2030 [5], and this figure is expected to exceed 693 million in 2045 [6]; the most dramatic increase is estimated to occur in developing countries [7].
The lifestyle of Chinese people is also changing rapidly, leading to a continuous increase in diabetes incidence and the obesity rate. The diabetes prevalence rate was lower than 1% in 1980 [8], reaching 5.5% in 2001 [9], 9.7% in 2008 [10], and 10.9% in 2013 [11]. According to the latest international diagnostic criteria, 11.6% of Chinese adults had diabetes in 2010, the prevalence of prediabetes was estimated at 50.1% [12], and approximately 114 million people had type 2 diabetes mellitus (T2DM).
The incidence of DN in China has surpassed that of primary glomerulonephritis and has become the first cause of chronic renal insufficiency, the first cause of hospitalization for chronic kidney disease in urban regions and the second leading cause of ESRD [22]. It is estimated that the number of patients with diabetesrelated chronic kidney disease in China has reached 24.3 million [22]. The number of DN-related deaths in the world also increased by 2.8-fold from 463000 in 1990 to 1731000 in 2013, which is the fastest increase among chronic diseases [23].
Most Chinese studies on DN have focused on inpatients in secondary and tertiary hospitals. For example, an investigation of 3469 inpatients with T2DM in Beijing, Shanghai, Tianjin, and Chongqing showed that the incidence of DN complications was 39.7% from 1991 to 2000. The prevalence of diabetic microvascular complications in large Chinese cities has reached or even exceeded the average level in developed countries.
Early diagnosis and treatment can reduce the mortality and disability rate of DN and improve the quality of life of patients. However, due to the insidious onset of DN, renal function damage is irreversible at the time of diagnosis of most patients. Therefore, the purpose of this study was to develop a valid and simple prediction tool using a questionnaire, physical examination, biochemical tests, and estimation for patients with type 2 diabetes mellitus (T2DM) to assess the risk of DN incidence [24].

Patient characteristics
All 3489 community-based patients with T2DM were divided into nondiabetic nephropathy (NDN) and DN groups according to the ratio of urinary microalbumin to uric creatinine (ACR), and 701 patients were diagnosed with DN (prevalence 20.10%). Among all participants, there were 1562 males (44.8%) and 1927 females (55.2%), the mean age was 64.54 ± 6.78 years [range 28-87 years], and the duration of disease was 9.48 ± 6.93 years. All data of patients, including demographics, lifestyle habits, physical examination results and biochemical test results, in the two groups are given in Table 1.

Feature selection
Of the demographics, lifestyle habits, physical examination results and biochemical test results, 19 features were reduced to 8 potential predictors on the basis of the results of 3489 patients (~2:1 ratio; Figures 1,2) and had nonzero coefficients in the least absolute shrinkage and selection operator (LASSSO) regression model. The coefficient of lambda was 0.002. These features included systolic blood pressure (SBP), diastolic blood pressure (DBP), fasting blood glucose (FBG), glycosylated hemoglobin A1c(HbA1c), total triglycerides (TGs), serum creatinine (SCR), blood urea nitrogen (BUN) and body mass index (BMI) ( Table 2).

Development of an individualized prediction model
The results of the logistic regression analysis among SBP, DBP, FBG, HbA1c, TGs, SCR, BUN and BMI are given in Table 3. A model that incorporated the above independent predictors was developed and presented as the nomogram (Figure 3).

Apparent performance of the DN incidence risk nomogram in the array
For the 3489 T2DM patients in the array, the calibration curve of the nomogram to predict DN risk showed medium agreement (Figure 4). The C-index of the nomogram was 0.744 (95% CI: 0.724-0.764), and the Cindex of internal verification was 0.737, which indicated that the model had medium prediction accuracy. The area under the receiver operating characteristic curve (AUC) was 0.744 ( Figure 5), and the cutoff AGING was 0.14, which indicates moderate performance. To summarize the results from the above verification, the nomogram of the model has medium prediction ability.

Clinical use
The decision curve analysis for the DN incidence risk nomogram is presented in Figure 6. The decision curve showed that if the threshold probability of a patient and a doctor is >20%, using this DN incidence nomogram to predict DN incidence risk for T2DM patients is more beneficial than the existing scheme. Within this range, the net benefit was comparable with several overlaps on the basis of the DN incidence risk nomogram.
In this study, according to AUC 0.744, cutoff 0.14, calibration, and decision curve analysis, the predictive   Figure 8). Therefore, the team developed a dynamic nomogram online APP to predict the DN of T2DM, the URL: https://doctorhu.shinyapps.io/DN_DynNomapp/. The APP is developed based on the characteristic indicators derived from the model. It can assist clinicians in diagnosing DN and can also be used in communitybased diabetes DN prevention [25].

DISCUSSION
Currently, nomograms are widely used as prognostic devices in oncology and medicine [26]. Nomograms depend on user-friendly digital interfaces, have improved accuracy, and provide more easily understood prognoses to aid better clinical decision making [27]. Our study was the first to apply a nomogram to the incidence of DN in T2DM patients.

Prevalence of DN among T2DM patients is consistent with other studies
In our study, the prevalence of DN in community-based patients with T2DM was 20.10%, which was consistent with some other studies [15].
According to a prospective survey of 66726 patients with T2DM in 28 countries, the prevalence of DN was 27.9%, with the highest prevalence in the Middle East and Russia and the lowest prevalence in South Asia [28]. In the United States, the National Health and Nutrition Survey  The decision curve showed that if the thresh old probability of a patient and a doctor is >20%, respectively, using this DN incidence risk nomogram in the current study to predict DN incidence risk adds more benefit than the intervention-all-patients scheme.
showed that the prevalence of DN in patients with diabetes is 23.7% [29,30]

Main risk factors for DN in T2DM patients
In the risk factor analysis, SBP, DBP, FBG, HbA1c, TGs, SCR, BUN and BMI were associated with DN incidence in T2DM patients. This nomogram suggested that lower SBP, DBP, FBG, HbA1c, TGs, SCR, BUN and BMI may be the key factors that can reduce the DN incidence risk of T2DM patients.
Systolic blood pressure and diastolic blood pressure SBP and DBP are the main diagnostic indicators of hypertension. This study found that SBP and DBP are risk factors for DN in patients with T2DM. It has been proven that among T2DM patients, hypertension is an independent risk factor for the occurrence and development of microvascular complications [34]. The prevalence of hypertension is higher in patients with DN. SBP and DBP were significantly higher in DN patients than in NDN patients (P <0.001). Boguslawa et al [35] found that for every 10 mmHg increase in blood pressure, the risk of microalbuminuria (MAU) in T2DM patients would increase by 23%. In Thailand, a nationwide crosssectional study including 55797 T2DM patients found that patients with hypertension had a 1.32-fold higher risk of microangiopathy than patients with normal blood pressure [36]. A European study of T2DM patients [37] found that the SBP and DBP in patients with MAU were significantly higher, and multivariate analysis showed that SBP was an independent risk factor for MAU in male patients with diabetes. In the UKPDS, a randomized clinical trial including 1148 T2DM patients [38], it was found that strict control of blood pressure could reduce the risk of microangiopathy by 37% (P=0.0092). A study conducted in Turkey with 202 T2DM patients in Zonguldak Atatürk State Hospital found that SBP (P=0.043) was independently related to logarithmically converted 24-hour UAER [39]. Leehey et al [40] confirmed that DN patients with mean SBP ≥140 mmHg had poor renal prognosis.  The cause of hypertension in diabetes mellitus patients may be associated with DN pathogenesis with small artery pressure causing renal arteriolar sclerosis, hemal wall thickening, luminal narrowing, increased glomerular filtration (increased plasma osmotic pressure during hyperglycemia, which promotes the entry of interstitial fluid into the blood vessels), which can directly lead to the presence of hyperfiltration (increased the glomerular filtration rate (GFR), and there is no accepted GFR standard yet) and therefore polyuria and hyperperfusion (increased plasma flow). The major early changes in DN include glomerulosclerosis, and hyperglycemia puts blood vessels in a hyperperfusion state. Moreover, capillary blood pressure is increased, and filtration pressure is increased, resulting in increased GFR, glomerular hyperperfusion and excessive filtration. Mesangial cells expand, epithelial cell foot processes fuse and produce dense droplets, and glomerular epithelial cells fall off the basement membrane, leading to glomerular fibrosis [41]. Small and medium-sized artery endothelial cells and smooth muscle cells undergo secondary hardening and damage. Hyperplasia after repair causes thickening of the basement membrane and high levels of glucagon in people with diabetes lead to the activation of the renin-angiotensin system (RAS) to speed up diabetic microvascular and major vascular complications, eventually leading to kidney damage and renal dysfunction [42]. Therefore, effective blood pressure reduction can significantly reduce the occurrence and development of DN [43], reduce the excretion rate of urinary albumin [44], and delay the occurrence of decompensated nephropathy. DN is occult, and there is no obvious decline in renal function in the early stage, but there has been a state of high perfusion and high filtration of the glomeruli, and damage to the glomeruli has begun [45,46].

Serum creatinine and blood urea nitrogen
SCR is the most commonly used index to assess whether the kidney is damaged. When the glomerulus is severely damaged and the GFR is decreased, SCR is increased [47]. However, because of the kidney's strong ability to compensate and form reserves, SCR remains at a normal level in the early stage of kidney injury; then, SCR increases rapidly and significantly, but SCR is also affected by sex, muscle mass and diet [48].

AGING
As the main nitrogen-containing substance excreted by the kidney, BUN can reflect the filtration function of the glomerulus to some extent, which is of great significance to the condition, course and prognosis of kidney disease. It was clearly proven that glomerular damage, tubular injury, inflammatory responses, and oxidative stress contribute to the development of DN.

Body mass index and total triglycerides
High BMI and TGs indicate that patients are obese, and such patients tend to have high blood lipid levels.
Increased blood TGs will increase blood viscosity, thus reducing blood flow velocity and inducing microvascular lesions. Under the condition of high glucose, cells will become hyperosmotic, and a large amount of extracellular fluid infiltration will cause cell edema, destroy cell structure, and directly or indirectly inhibit the function of glomerular and tubular cells. The hyperfiltration state of the glomerulus is a direct factor leading to proteinuria. Moreover, the metabolic demand of the kidney in obese people is higher than that of nonobese people. Both anthropometric and biochemical indices of adiposity and, to a smaller extent, elevated blood pressure were implicated as major determinants of glomerular hyperfiltration [49,50]. To meet the higher metabolic demand, renal glomerular arterioles continue to expand, blood perfusion continues to increase, GFR increases, and the kidney's automatic regulation ability is impaired. Simultaneously, the pressure in glomerular capillaries increases, the tension of the tube wall increases, and the capillaries become thicker, making the distribution of podia cells relatively sparse, leading to the decrease in renal function [51]. Obesity further increases the risk of DN and ESRD [52].
It was found in the DCCT/EDIC study that lower TGs were associated with a reduced risk for progression from MAU to clinical albuminuria (CAU) or ESRD [53]. A cross-sectional epidemiological study in Singapore [54] showed that obesity significantly increased the risk of DN, with the risk ratios for DN in groups with a BMI 23-24.9, 25-29.9 and ≥ 30 being 3.12, 2.49 and 3.70 (95% CI: 2.13~6.42), respectively, compared with the group with a BMI <23. Cohen et al [55] further confirmed that overweight and obesity are risk factors for chronic kidney disease (CKD). The multivariate model analysis of Kwan JM et al [56] found that obese patients had a significantly increased risk of adverse transplant outcomes, including delayed graft function recovery, graft failure, urinary protein and acute rejection. Hyperleptinemia caused by obesity is also involved in the occurrence of DN: leptin can stimulate the proliferation of glomerular endothelial cells and glomerular mesangial cells, as well as the proliferation of tubulointerstitial cells and matrix components, by inducing oxidative stress of glomerular endothelial cells, leading to the occurrence of glomerular sclerosis. In addition, obesity increases sympathetic excitability and hyperactivity of the RAS, which is also associated with the occurrence of DN. The Look AHEAD Research Group [57] found that weight loss, namely, calorie restriction and increased exercise, could significantly reduce the occurrence and development of DN in obese or overweight T2DM patients. A randomized controlled study in Hong Kong showed that effective comprehensive management of T2DM patients could effectively reduce the risk of ESRD and death [58].

Fasting blood glucose and glycosylated hemoglobin A1c
FBG is the most basic indicator of the development of T2DM. HbA1c can basically reflect the average blood glucose level in the past two to three months, and longterm hyperglycemia results in protein nonenzymatic glycosylation, leading to systemic microvascular damage and the promotion and development of microangiopathy [59]. It was confirmed that strict glycemic control can significantly reduce the risk of diabetic microvascular complications [60]. Epidemiological analysis of the Diabetes Control and Complications Trail [61] showed that for every 10% reduction in HbA1c, there was a 25% and 44% reduction in the risk of MAU and CAU, respectively, in patients with type 1 diabetes mellitus. Similarly, the UKPDS showed that for every 1% decrease in HbA1c, there was a 37% reduction in the risk of T2DM microvascular complications. The UKPDS was a clinical study that included 10 years of blood glucose intervention in patients with newly diagnosed T2DM [62] and showed that HbA1c (7.0%) in the intervention group was 11% lower than that in the traditional treatment group (7.9%), and the risk of microvascular complications was reduced by 25%. The incidence of DN was reduced by 33% at the 12-year follow-up. Ten years after the trial, the benefits of intensive glycemic control persisted after treatment [35]. Moreover, a study in Taiwan found that HbA1c significantly promoted DN progression in T2DM patients with high levels of normal albuminuria [63]. Moreover, a meta-analysis showed that patients with HbA1c<7% had less MAU than patients with HbA1c<7.9% [64]. Chronic hyperglycemia in diabetic patients accelerates nonenzymatic glycosylation, leading to increased hypoxia of tissues and the release of vasoactive substances, which result in hemodynamic changes and are important risk factors for DN. The effect of blood glucose on DN is reflected in three aspects. The first is poor long-term blood glucose control. Ficociello et al [65] conducted a 10-year follow-up study and found that diabetic patients with persistent hyperglycemia or poor blood glucose control had a 3 to 8 times higher risk of DN development or progression. Persistent hyperglycemia can induce glomerular mesangial hyperplasia, resulting in basement AGING membrane thickening and glomerular capillary wall sclerosis, ultimately leading to the occurrence of proteinuria. On the other hand, persistent hyperglycemia causes the glomerulus to be in a state of hyperperfusion and hyperfiltration, resulting in capillary sclerosis and changes in vascular permeability, which ultimately leads to plasma protein deposition on capillary walls, progressive vascular hyperplasia, hyaline degeneration and intravascular thrombosis [66]. Second, blood glucose fluctuates greatly. The UKPDS report confirmed that blood glucose fluctuation was closely related to microvascular complications and was a risk factor for the occurrence or progression of DN, independent of blood glucose and HbA1c, even exceeding the risk associated with hyperglycemia alone [67]. A newly published study confirmed that fasting glucose variation and HbA1c variation are risk factors for progression to ESRD among patients with diabetes [68]. Moreover, studies have shown that not only hyperglycemia but also insulin resistance or hyperinsulinemia before the onset of diabetes begins to result in glomerular hypertrophy and thickened capillary basement membranes [69].

Limitations
There are several limitations in our current study. First, our acquired data might have a low representation of males and only a partial representation of T2DM patients. The cohort was not representative of all Chinese patients with T2DM. Our study mainly focused on community-based T2DM patients in Shanghai, and patients who did not participate in community T2DM management were not included. The number of males and females included in our study was quite different. Second, risk factor analysis did not include all potential factors that affect DN incidence. Some possible aspects of DN incidence risk were not thoroughly assessed, such as drug treatment, exercise and other factors. Third, although the robustness of our nomogram was examined extensively with internal validation in the same population, there is a lack of external verification in other T2DM populations in other regions and countries. The nomogram needs to be externally evaluated in wider T2DM populations.

CONCLUSIONS
This study developed a novel nomogram with a relatively moderate accuracy to help clinicians access the risk of DN incidence in T2DM patients. With an estimate of individual risk, clinicians and patients can take more necessary measures in regard to lifestyle monitoring and medical interventions. This nomogram requires broader external clinical validation, and further research is needed to determine whether individual interventions based on this nomogram will reduce DN incidence risk and improve treatment outcomes.
The prevalence rate of DN was 20.10% among community-based patients with T2DM in Shanghai, which was consistent with previous research. We found that HbA1c was one of the most important risk factors for DN, which suggests that the government should establish more effective and extensive DN screening and T2DM management programs at the community level, with special attention to interventions and monitoring of risk indicators, such as HbA1c, FBG, and hypertension. Moreover, patients should strengthen their self-management of T2DM, especially those with hypertension and obesity, and T2DM patients need to pay close attention to the occurrence of DN.

Patients
The project started in September 2014 and was led by The patients included in the study were diagnosed with T2DM in the hospital through medical records and medical history and participated in community diabetes management. The criteria defined by the World Health Organization in 1999 were used as the diagnostic criteria for diabetes [70]. The baseline data of this paper were collected from September 2014 to April 2015. The follow-up data of this paper were collected from October 2018 to April 2019. We obtained written informed consent from all participants before they were included in the study. A total of 3661 T2DM patients were enrolled in the study in the baseline survey. The inclusion criteria were as follows: patients with T2DM who were no more than 75 years old, had a district household registration or were a permanent resident (living in the community for at least 6 months), did not have serious chronic diseases (according to Chinese medical insurance regulations: malignant tumor, uremia dialysis, anti-rejection immunomodulator after organ transplant, chronic pulmonary heart disease, active tuberculosis, neurological deficit caused by sequelae of cerebrovascular disease, myocardial infarction/plug, chronic moderate/severe viral hepatitis, stage III highrisk and very-high-risk hypertension, chronic aplastic AGING anemia, systemic lupus erythematosus, schizophrenia, hemophilia), and participated in the voluntary screening of diabetic complications. Among all participants, 43 were excluded due to incomplete data, and 129 were lost, the missed interview rate was 3.52%; finally, 3489 participants were included in the analysis. From September 2014 to April 2019, a total of 3489 T2DM patients in 6 communities were obtained baseline and follow-up data, and the final DN was used as the gold standard ( Figure 9).

Investigation methods
We determined the number of patients with T2DM in each community through the health center hospital and community health service center of Shanghai University of Traditional Chinese Medicine and invited the patients to the community service center for examination one week in advance and examined 100 patients per day. We completed the questionnaire survey, physical examination, and the collection of fasting blood, postprandial blood and urine for each patient. At seven o'clock in the morning, the patients' fasting blood was collected, and breakfast was provided. Two hours later, the patients' postprandial blood and urine were collected. During the waiting period, faceto-face questionnaires and physical examinations were administered, and participants were prompted by welltrained investigators to sign patient consent and provide information. The questionnaire survey included basic sociodemographic information, diabetes mellitus status (course of disease, complications, DN), and daily lifestyle habits (smoking, alcohol consumption). The physical examination included height, weight, waist circumference, hip circumference, SBP and DBP. BMI and waist-to-hip ratio (WHR) were calculated. On the day of the examination, blood and urine samples of all patients involved were sent to the laboratory for further processing. The blood samples were separated from the plasma and serum by centrifugation and then sent to a testing center for biochemical testing to obtain corresponding data. The laboratory examination included FBG, postprandial blood glucose (PBG), HbA1c, low-density lipoprotein (LDL)high-density lipoprotein (HDL), TGs, total cholesterol (TC), SCR, BUN, uric acid (UA), uric creatinine (UCR), and urinary microalbumin (UMA). ACR was calculated.
The diagnosis of DN was divided into three stages according to ACR (ACR less than 30 mg/g was defined as nonalbuminuria (NAU), 30 mg/g < ACR < 300 mg/g was defined as MAU, and ACR greater than or equal to 300 mg/g was defined as CAU [71]. Both MAU and CAU indicated DN in this study.

Statistical analysis
All data, including demographics, lifestyle habits, and data from physical examinations and biochemical tests, are expressed as frequencies and counts (%). Statistical analysis was performed using R software (version 3.6.2; https://www.R-project.org). LASSO method is a shrinkage and variable selection method for linear regression models and is used to simplify high dimensional data [72]. To obtain the subset of predictors, LASSO regression analysis minimizes the prediction error for a quantitative response variable by imposing a constraint on the model parameters that cause regression coefficients for some variables to shrink toward zero [73]. Variables with a regression coefficient equal to zero after the shrinkage process are excluded from the model, while variables with nonzero regression coefficient variables are most strongly associated with the response variable. Based on the type measure of -2 log-likelihood and binomial family, the LASSO regression analysis running in R software runs 10 times K cross-validation for the centralization and normalization of included variables and then picks the best lambda value [74]. "Lambda.lse" gives a model with good performance but the least number of independent variables. Therefore, the LASSO method was used to select the optimal predictors from the present risk factors from the T2DM patients in the study. Then, multivariable logistic regression analysis was used to build a prediction model by introducing the features selected in the LASSO regression model [75]. Additionally, the forest plot was drawn to describe the P-value, OR and 95% CI of selected validation visually. The statistical significance levels were all two-sided. By introducing all the selected features and analyzing the statistical significance levels of the features, the statistically significant predictors were applied to develop a model of DN incidence risk prediction for patients with T2DM [76].
Calibration curves were used to evaluate the calibration of the DN incidence risk nomogram, accompanied by the Hosmer-Lemeshow test. A significant test statistic implies that the model does not calibrate perfectly [77]. To quantify the predictive ability of the nomogram, the Cindex measurement was performed, and internal validation was performed using bootstrapping. Bootstrap samples are random samples drawn with replacement from the original sample, and 1000 bootstraps were repeatedly fitted in the model to evaluate the performance [78]. The receiver operating characteristic curve and the area under the receiver operating characteristic curve was used to perform medium discrimination of the quality of the risk nomogram to separate true positives from false positives [79]. Decision curve analysis was used to determine the clinical practicability of nomograms based on the net benefit according to different threshold probabilities in T2DM patients [80].
Both the NRI and IDI were used to compare the diagnostic capabilities of the two indicators by assessing whether one indicator improved the diagnostic accuracy of the other indicator.
NRI quantitatively indicates to what extent the diagnostic accuracy of one indicator is better than the accuracy of another indicator. , and the IDI test idea is similar to AUC and considers the overall situation of the indicator under different boundaries, allowing it to represent the overall improvement in the indicator (model). The application of NRI and IDI quantitatively solves the problem of the comparison of the diagnostic efficacy of the two models and can quantitatively evaluate the improvement in the diagnostic efficacy of one indicator over the other. The NRI and IDI indicators can clearly show the proportion of research subjects that are correctly reclassified (discriminated), and they are widely used in research as clinical prediction models.

Ethics statement
All participants were carefully informed about the protocol and provided written informed consent before inclusion in the study. The study was performed in accordance with the principles of the Declaration of Helsinki.

AUTHOR CONTRIBUTIONS
Fan Hu and Yuhong Hu designed the study. Rong Shi acquired data. Yuhong Hu and Ruohui Mo collated the data. Fan Hu and Yuhong Hu performed the analysis, interpretation and discussion of the data. Yuhong Hu wrote the article. All the authors approved the final version of the manuscript.