Enteral nutrition feeding in Chinese intensive care units: a cross-sectional study involving 116 hospitals

Background

There is a lack of large-scale epidemiological data on the clinical practice of enteral nutrition (EN) feeding in China. This study aimed to provide such data on Chinese hospitals and to investigate factors associated with EN delivery.

Methods

This cross-sectional study was launched in 118 intensive care units (ICUs) of 116 mainland hospitals and conducted on April 26, 2017. At 00:00 on April 26, all patients in these ICUs were included. Demographic and clinical variables of patients on April 25 were obtained. The dates of hospitalization, ICU admission and nutrition initiation were reviewed. The outcome status 28 days after the day of investigation was obtained.

Results

A total of 1953 patients were included for analysis, including 1483 survivors and 312 nonsurvivors. The median study day was day 7 (IQR 2–19 days) after ICU entry. The proportions of subjects starting EN within 24, 48 and 72 h after ICU entry was 24.8% (84/352), 32.7% (150/459) and 40.0% (200/541), respectively. The proportion of subjects receiving > 80% estimated energy target within 24, 48, 72 h and 7 days after ICU entry was 10.5% (37/352), 10.9% (50/459), 11.8% (64/541) and 17.8% (162/910), respectively. Using acute gastrointestinal injury (AGI) 1 as the reference in a Cox model, patients with AGI 2–3 were associated with reduced likelihood of EN initiation (HR 0.46, 95% CI 0.353–0.599; p < 0.001). AGI 4 was significantly associated with lower hazard of EN administration (HR 0.056; 95% CI 0.008–0.398; p = 0.004). In a linear regression model, greater Sequential Organ Failure Assessment scores (coefficient – 0.002, 95% CI – 0.008 to − 0.001; p = 0.024) and male gender (coefficient – 0.144, 95% CI – 0.203 to − 0.085; p < 0.001) were found to be associated with lower EN proportion. As compared with AGI 1, AGI 2–3 was associated with lower EN proportion (coefficient – 0.206, 95% CI – 0.273 to − 0.139; p < 0.001).

Conclusions

The study showed that EN delivery was suboptimal in Chinese ICUs. More attention should be paid to EN use in the early days after ICU admission.

Patients requiring intensive care unit (ICU) admission are at increased risk of death, owning to a variety of life-threatening conditions such as respiratory failure, shock, severe infection and multiple organ dysfunction [1, 2]. Metabolic response to critical illness is characterized by accelerated catabolism, resulting in wasting and negative nitrogen balance [3]. Nutritional support during the catabolic phase will not only lead to positive nitrogen balance but also prevent weakness and, eventually, multiple organ failure and death [4, 5]. Regarding the route of energy administration, enteral feeding has been proven to be superior to the parenteral nutrition in terms of the outcomes such as nosocomial infection, medical cost saving and even mortality rate [6,7,8,9]. Thus, enteral nutrition (EN) feeding is the primary choice of nutrition therapy in current clinical practice; a careful and well-monitored approach is proposed based on the risk of poor tolerance [10]. Critically ill patients usually have compromised gastrointestinal function, making it difficult to increase EN to a target. Other factors that can delay achieving a target EN proportion include compromised gastric function, abdominal surgery and unstable hemodynamics [11]. The clinical practices of EN feeding vary substantially across different regions and hospitals.

China has the world’s largest population of critically ill patients, and the clinical practice of EN feeding varies substantially across regions and hospitals. However, to the best of our knowledge, there is a lack of large-scale data on the practice of enteral nutrition in ICUs. Previous small studies have shown that the proportion of EN was as low as 40% on day 2, which could probably be improved with implementation of an EN feeding protocol [12, 13]. Compliance to the clinical practice guidelines was also found to be suboptimal in neurological ICUs [14]. Due to the lack of epidemiological data on clinical practice of EN in Chinese ICUs, a nationwide cross-sectional study was performed. The objective of the study was to provide epidemiological data on the clinical practice of EN feeding in Chinese ICUs. Also, factors associated with the initiation and proportion of EN were investigated. We hypothesized that the severity of illness and gastrointestinal function would influence the timing and quantity of EN feeding.

Study design

This cross-sectional study was conducted in 118 ICUs of 116 hospitals covering most provinces of mainland China on April 26, 2017. The types of ICUs in this study vary, including general ICUs, surgical ICUs, medical ICUs, respiratory ICUs, coronary ICUs, emergency ICUs and so on. Demographic and clinical variables of patients on April 25, 2017 (hereafter referred to as “the study day”) were recorded. The outcome of the study subjects on the 28th day after investigation was obtained. Data collection was performed via a customized website. The study was approved by the ethics committee of Jinling Hospital (approval no. 2017NZKY-010-01).

Study population

At 00:00 on April 26, 2017, all patients in the study ICUs were included. Exclusion criteria were: patients or their surrogate refused to participate in the study; and patients whose 28-day outcomes were not recorded.

Variables

Demographics such as age, gender, source of admission, weight (past weight and temporal weight) and height were obtained. The past weight was defined as the body weight when patients were healthy or before 1 week, whereas the temporal weight was the body weight on the study day or within 1 week. The dates of hospital admission and ICU admission were reviewed.

Routine laboratory tests of blood samples were measured, for example, C-reactive protein (CRP), percentage of lymphocyte, albumin, maximum and minimum blood glucose, and arterial lactate. The level of consciousness was estimated using the Glasgow Coma Scale (GCS). If a patient was under sedation, the sedative was reduced or discontinued to make an appropriate assessment. Vital signs such as mean blood pressure and some particular treatments such as using vasopressors were recorded for Sequential Organ Failure Assessment (SOFA) and Acute Physiologic and Chronic Health Evaluation II (APACHE II) scores. The maximum and minimum blood glucose (BG) values on the study day were also reported. If a patient had more than one measurement of these variables on the study day, the worst one was chosen.

The date for the start of EN was retrospectively recorded in the study, but details regarding EN delivery were only recorded on the study day. Gastrointestinal function was evaluated as normal or mild injured, moderate to severe injured and failure referring to the acute gastrointestinal injury (AGI) grading. Enteral feeding volume and energy density of EN were recorded in the units of milliliters and kilocalories per milliliter, respectively. The EN proportion of target was computed as the ratio of the energy provided by EN to the estimated energy target (= past body weight × 25 kcal/kg). The nutrition could be warmed and diluted, and the dilution effect was considered in calculating the targets–achievement ratio. The enteral feeding routes included gastric feeding, postpyloric feeding, percutaneous gastrostomy/jejunostomy (PEG/J), jejunostomy and others [15]. The EN administration styles included continuous pump, gravity, intermittent feeding and other styles. The body positions were recorded as elevation of the bed head by 30°, sitting or other positions. The frequency of gastric residual volume (GRV) measurement and the largest GRV were recorded. In a case when gastrointestinal decompression was used, GRV was not measured and the gastrointestinal decompression volume was measured. Abdominal pressure was estimated using the bladder pressure. Duration of EN feeding was defined as the hours of the study day during which EN was administered. The times and causes of EN discontinuation during the study day were recorded. Intolerance to EN included nausea, vomiting, aspiration, abdominal pain and distension, and diarrhea. The presence of nausea and vomiting were recorded as yes, no, unable to judge and so on. The severity of abdominal pain and distension were grade-evaluated. Stool frequency and total volume on the study day were recorded.

The outcome status at 28 days after the study day was recorded as the endpoint. The outcome variable was labeled as alive, dead and lost, corresponding to vital status of a subject. If a patient’s condition on day 28 was not available, he/she was considered lost to follow-up.

Patients with missing outcomes

Missing values were present in the dataset. We performed single imputation for missing values [16]. For patients with missing values on outcome, their baseline characteristics were also reported and compared with those with a complete dataset.

Statistical analysis

Continuous variables were expressed as mean and standard error, or median and interquartile range, as appropriate. Categorical data were expressed as the number and proportion. Because there were two groups in the study, analysis of variance was performed for normally distributed data and the Kruskal–Wallis test was used for nonnormal data. Categorical data were compared using the chi-square test [16, 17].

The initiation of EN was considered as survival data and a multivariable Cox proportional hazard regression model was performed to investigate factors associated with the initiation of EN [17]. Variables included in the model were the source of admission, age, SOFA and APACHE II scores, gender, AGI grading and GCS score. The selection of variables was based on the literature finding that severity of illness was associated with EN delivery, and demographics of age and gender were also included.

The factors associated with the proportion of EN on the study day were investigated in a multivariable linear regression model [18]. Variables were selected using a stepwise backward elimination method and the best model was judged by the Akaike information criterion (AIC) [19]. The final model included variables such as the source of admission, age, APACHE II score, SOFA score, gender, AGI grading and GCS score. The cross-sectional study was subject to length bias, in which patients who had a longer length of stay in the ICU were more likely to be included in the study. Thus, the patients included in the study were not an unbiased sample of the target population. To adjust for this length bias, we assigned higher weights to patients with shorter LOS in the ICU [20].

All statistical analyses were performed using R (version 3.3.2) [21]. Two-tailed p < 0.05 was considered to indicate statistical significance.

Characteristics of enrolled patients

The baseline characteristics and variables for the study day are presented in Table 1. A total of 1953 patients were included for analysis, including 1483 survivors and 312 nonsurvivors. The median study day was day 7 (IQR 2–19 days) after ICU entry, and there was no difference between survivors and nonsurvivors. Nonsurvivors were significantly older than survivors (median 70 vs 66 years; p = 0.005 for two-group comparison). There were more males among both survivors (68%) and nonsurvivors (63%). Most of the study patients came from emergency departments (35%), followed by surgery (29%) and internal medicine (19%). The distribution of admission wards was similar for survivors and nonsurvivors. Past and temporal weight (median 65 vs 65 kg, 65 vs 63 kg; p = 0.369, 0.198 for two-group comparison) and EN start day (median 1 (0, 2.5) vs 1 (0, 2) days; p = 0.912 for two-group comparison) showed no difference between the two groups.

For laboratory variables, CRP (median 40.1 vs 74.1 mg/dl; p < 0.001), maximum BG (median 10 vs 11.1 mmol/l; p < 0.001), minimum BG (median 6.3 vs 6.8 mmol/l; p < 0.001), lactate (median 1.3 vs 1.7 mmol/l; p < 0.001) and total bilirubin (median 12 vs 14 mmol/l; p < 0.001) were significantly lower in survivors than those in nonsurvivors, while the GCS was significantly higher (median 13 vs 9; p < 0.001). As expected, the SOFA score (median 7 vs 5; p < 0.001) and APACHE II score (median 20 vs 16; p < 0.001) were higher in the nonsurvivors than that in the survivors.

EN delivery

There were 1440 patients in AGI 1 (86%), 210 patients in AGI 2–3 (12%) and 33 patients in AGI 4 (2%). The severity of AGI was associated with mortality, as there were significantly more patients with AGI 1 in the survivors (88% vs 78%; p < 0.001) than in the nonsurvivors, and more patients had AGI 2–3 (18% vs 11%; p = 0.001) and AGI 4 (4% vs 1%) in the nonsurvivors than in the survivors.

There were 1270 patients using EN (66%) on the study day and there was no difference between survivors and nonsurvivors (67% vs 66%). The mean total volume of EN was 1000 ml in the two groups. The energy density of EN was 1 or 1.5 kcal/ml for different formulas. However, the proportion of starting EN within 24, 48 and 72 h after ICU entry was 23.9% (84/352), 32.7% (150/459) and 37.0% (200/541), respectively. The proportion of subjects receiving > 80% estimated energy target within 24, 48, 72 h and 7 days after ICU entry was 10.5% (37/352), 10.9% (50/459), 11.8% (64/541) and 17.8% (162/910), respectively. Data on the proportion of EN were only available for the study day. Thus, the calculation of EN proportion within 24, 48, 72 and 7 days only included patients for whom the study days were within the relevant time window (a subset of all patients) and not from all recruited patients. Figure 1 plots the number of patients with different EN proportions against the study days (e.g., with reference to the day of ICU entry), showing that the use of EN increased over time.

Number of patients with different EN proportions stratified by different study days. Nearly two thirds of patients did not use EN within 2 days after ICU entry. Both the proportion of patients using EN and EN proportions of total energy target increased over time, but less than half of patients received > 60% total target after first week in ICU. d days, EN enteral nutrition, ICU intensive care unit

Full size image

Gastric feeding was used in 85% patients, followed by postpyloric feeding (11%), PEG/J (2%) and jejunostomy (2%). Continuous pump (88%) was the predominant EN infusion style, followed by intermittent pump (7%) and infusion by gravity (3%). The majority of patients (83%) had their bed head elevated more than 30°. Warming was used in 47% of patients and EN was not diluted for the majority of cases (96%). The frequency and the total volume of GRV measurement were 2 (IQR 0–4) times and 0 (IQR 0–50) ml. The gastrointestinal depression volume and abdominal pressure were 100 (IQR 40–250) ml and about 14 (IQR 11–18) mmHg. EN was not discontinued for the majority of patients in the study day; the reasons for EN discontinuation were intolerance (52 patients, 5%), examination (98 patients, 9%) and others (52 patients, 5%). However, none of these EN-related parameters was significantly associated with mortality outcome, except for the proportion of gastrointestinal decompression (16% vs 25%; p < 0.001).

Factors associated with EN initiation

In the Cox proportional hazard model (Table 2), we found that the AGI and GCS were associated with initiation of EN. Using AGI 1 as the reference, patients with AGI 2–3 were associated with reduced likelihood of EN initiation (HR 0.46, 95% CI 0.353–0.599; p < 0.001). AGI 4 was significantly associated with lower hazard of EN administration (HR 0.056, 95% CI 0.008–0.398; p = 0.004). Although a greater value of GCS was statistically associated with reduced likelihood of EN initiation, the clinical significance was marginal (HR 0.945, 95% CI 0.921–0.969; p < 0.001).

Factors associated with proportion of EN delivery

In the multivariable linear regression model investigating factors associated with the EN proportion (Table 3), greater SOFA scores (coefficient – 0.002, 95% CI – 0.008 to − 0.001; p = 0.024) and male gender (coefficient – 0.144, 95% CI – 0.203 to − 0.085; p < 0.001) were found to be associated with lower EN proportion. As compared with AGI 1, AGI 2–3 was associated with lower EN proportion (coefficient – 0.206, 95% CI – 0.273 to − 0.139; p < 0.001).

This was the largest cross-sectional study in mainland China covering 118 ICUs in 116 hospitals. The study enrolled 1953 patients, including 1483 survivors and 312 nonsurvivors. The results showed that EN delivery was suboptimal in China, evidenced by only 32.7% of critically ill patients received EN feeding within 48 h after ICU admission. Only a minority of patients (17.8%) received 80% of the estimated energy target by EN within 7 days after ICU entry. In the multivariable analysis, AGI was independently associated with the initiation of EN, indicating that gastrointestinal dysfunction was the major obstacle prohibiting EN delivery. The severity of illness as reflected by the SOFA score was negatively associated with the EN proportion.

In a worldwide study involving 46 countries and 880 units, enteral feeding was prescribed to only 10% of patients on the first day but this number increased to more than 40% of patients after 5 days [22]. In our study, 23.9% patients received EN on the first day and this percentage increased to 37% within 3 days. In another study conducted in Latin America, 59.7% of patients received > 90% estimated daily target within 24 h after ICU admission [23]. The study included parenteral nutrition, enteral nutrition or a combination, and thus the proportion was much higher than that in our study. The time to achieve > 80% estimated energy target was longer in the present study as compared with other reports. For example, Yip et al. [24] reported that the mean time to achieve > 80% target was 1.8 days (SD 1.5 days). However, the proportion of patients receiving > 80% energy target in the first 2 days was only 10.9% (50/459) in our study. Meanwhile, the median volume of EN in our study was 1000 exactly, suggesting that a more standard 1 L per 24 h is used and not related to target or body weight or BMI. This seems more cultural, nonprotocolized. According to the recent clinical practice guidelines [25], EN should be initiated within 24–48 h for critically ill patients who are unable to maintain volitional intake. However, this target is far from being reached in most Chinese hospitals. A cohort study from China reported that the EN feeding protocol was able to increase the proportion of EN on day 2 (41.8 ± 22.3% vs 50.0 ± 28.3%; p = 0.006) [13]. Thus, every effort should be made to standardize the EN feeding protocol to prompt EN delivery for critically ill patients.

Factors associated with EN initiation were investigated in the study. The results showed that patients with advanced AGI stage were less likely to receive EN, which was consistent with other reports. In a multicenter study, AGI grading was found to be associated with EN intolerance and the discontinuation of EN feeding [26]. Although our study is cross-sectional in design and cannot deduce the causal relationship between AGI grade and EN discontinuation, we could propose that EN discontinuation might be attributable to the compromised gastrointestinal function. However, this practice may not be supported by empirical evidence, because there is evidence showing that advanced AGI grade should not be a contraindication of EN feeding. On the contrary, Jin et al. [27] showed that the AGI grade was improved in the early EN group (F-statistic = 4.468; p < 0.05), indicating that EN delivery for patients with moderate gastrointestinal impairment may help to accelerate the recovery of gastrointestinal function. Actually, any association between higher AGI and/or higher SOFA score and lower EN given indicates that increasing severity of disease is associated with decreasing EN tolerance as well as poorer outcome. It was counterintuitive to see that a greater value of GCS was statistically associated with reduced likelihood of EN initiation, but the clinical significance was marginal (HR 0.945, 95% CI 0.921–0.969; p < 0.001). This could be explained by the fact that for patients with low GCS caused by trauma or stroke, we usually insert a GI tube on admission and start EN very early. After all, gastrointestinal dysfunction or shock, one of the major factors limiting EN use, is less common in stroke or traumatic brain injury patients.

In a linear regression model, we investigated factors associated with the proportion of the estimated energy requirement delivered via EN. The results showed that the SOFA score, gender and AGI were independently associated with the proportion. No wonder that a higher SOFA score is an obstacle of EN increments, as it is a measurement of organ dysfunction for critically ill patients and circulatory failure is an important component of the total score [28]. The guidelines also suggest a delay of EN delivery in patients with circulatory failure [25]. A novel finding in the study was that male patients were less likely to reach the estimated energy target as compared to females. Some may argue that this can be a false positive finding due to multiple testing. However, the effect size was large with male patients 26% lower in proportion than female patients (p < 0.001). Such a large effect size cannot be simply explained by random errors [29]. Probably, the physicians were more likely to give a fixed dose to patients without considering their body weight. Since men’s body weight is usually greater than women’s, they might receive a lower proportion of EN than female patients. In the linear model, AGI appeared to be an independent predictor of EN proportion. It is reasonable that physicians reduce the amount of EN in the presence of signs and symptoms that exaggerate the AGI grading [30]. According to current guidelines, AGI was not a factor prohibiting the initiation of EN. For a patient with greater AGI, physicians can initiate but maybe not achieve the target so easily. Initiation has to do with protocols, awareness of the importance, preserving gut function and so on. Furthermore, the study showed that GI intolerance was a very uncommon reason to interrupt EN. Thus, the problem is not AGI but the culture of not starting EN.

Several limitations must be acknowledged. First, the study was a cross-sectional study that causal inference cannot be confirmed. For example, we could only conclude that there was an association between AGI and EN delivery, the causal relation cannot be well defined. While a high AGI score may prohibit physicians to add EN, it is also possible that EN delivery exacerbates the already compromised gastrointestinal function. Second, the study only recorded data on the study day and the EN delivery and gastrointestinal function before and after the study day remain unknown, which precluded the analysis of a temporal trend of the EN delivery. However, since we have a large number of patients, and the study day involved all consecutive days after ICU admission, this gives us an opportunity to investigate the current trend of EN delivery.

The study showed that EN delivery was suboptimal in Chinese ICUs, and that EN use in the early days after ICU admission should be paid more attention. One reason for this is that most Chinese hospitals lack a standardized EN feeding protocol for critically ill patients.

AGI:

Acute gastrointestinal injury

APACHE II:

Acute Physiology and Chronic Health Evaluation II

BG:

Blood glucose

CI:

Confidence interval

CRP:

C-reactive protein

EN:

Enteral nutrition

GCS:

Glasgow Coma Score

GRV:

Gastric residual volume

HR:

Hazard ratio

ICU:

Intensive care unit

IQR:

Interquartile range

PEG/J:

Percutaneous endoscopic gastrostomy/jejunostomy

SOFA:

Sequential Organ Failure Assessment

The authors thank Unimed Scientific Inc. for their efforts and dedication in the development of the website and database. The authors thank Jingyi Wu and Xuemei Gong for their help in data gathering.

Funding

This study was funded by The Standard Treatment to Key Diseases of Jiangsu Social Development Programme (BE2015685).

Availability of data and materials

The datasets generated and analyzed in this article are not publicly available due to health privacy concerns. However, they are available from the corresponding author and will be obtainable by the public when the database construction is complete.

Affiliations

1. Nanjing General Hospital of Nanjing Military Command, No.305 Zhongshan East Road, Nanjing, 210002, China

   • Juan Xing
   • , Lu Ke
   • , Jing Zhou
   •  & Weiqin Li

2. Department of emergency medicine, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou, China

   • Zhongheng Zhang

3. Henan Provincial People’s Hospital, Zhengzhou, China

   • Bingyu Qin

4. Zhongshan People’s Hospital, Zhongshan, China

   • Hongkai Liang

5. Shandong University Qilu Hospital, Jinan, China

   • Xiaomei Chen

6. Changzhou No.2 People’s Hospital affiliated to Nanjing Medical University, Nanjing, China

   • Wenming Liu

7. Jilin University First Hospital, Changchun, China

   • Zhongmin Liu

8. Xiangya Hospital Central South University, Changsha, China

   • Yuhang Ai

9. Guizhou Medical University affiliated hospital, Guiyang, China

   • Difeng Wang

10. Wuxi People’s Hospital, Wuxi, China

    • Qiuhui Wang

11. Hubei Provincial People’s Hospital, Wuhan, China

    • Qingshan Zhou

12. Tai’an Central Hospital, Tai’an, China

    • Fusen Zhang

13. First Affiliated Hospital of Nanchang University, Nanchang, China

    • Kejian Qian

14. Third Military Medical University Daping Hospital, Chongqing, China

    • Dongpo Jiang

15. China Medical University Second Affiliated Hospital, Shenyang, China

    • Bin Zang

16. Guangzhou Medical University First Affiliated Hospital, Guangzhou, China

    • Yimin Li

17. Sichuan Provincial People’s Hospital, Chengdu, China

    • Xiaobo Huang

18. Qingdao Municipal Hospital Group, Qingdao, China

    • Yan Qu

19. Jining First People’s Hospital, Jining, China

    • Yinguang Xie

20. Guangzhou First Municipal People’s Hospital, Guangzhou, China

    • Donglin Xu

21. Xiehe Affiliated Hospital of Fujian Medical University, Fuzhou, China

    • Zhiqiang Zou

22. Dazhou Central Hospital, Dazhou, China

    • Xiangde Zheng

23. Inner Mongolia People’s Hospital, Huhehaote, China

    • Jianbo Liu

24. Zhejiang University School of Medicine Sir Run Run Shaw Hospital, Hangzhou, China

    • Feng Guo

25. Qindao University Medical College Affiliated Yantai Yuhuangding Hospital, Qingdao, China

    • Yafeng Liang

26. Tianjing People’s Hospital, Tianjin, China

    • Qiang Sun

27. Tianjing First Central Hospital, Tianjin, China

    • Hongmei Gao

28. Tangshan Gongren Hospital, Tangshan, China

    • Yang Liu

29. Southern Medical University Zhujiang Hospital, Guangzhou, China

    • Ping Chang

30. North China University of Science and Technology Affiliated Hospital, Tangshan, China

    • Aibin Ceng

31. Dalian Central Hospital, Dalian, China

    • Rongli Yang

32. Peking University Third Hospital, Beijing, China

    • Gaiqi Yao

33. Anhui Medical University Second Affiliated Hospital, Hefei, China

    • Yun Sun

34. Wenzhou Medical University First Affiliated Hospital, Wenzhou, China

    • Xiaorong Wang

35. Shanxi Provincial People’s Hospital, Taiyuan, China

    • Yi Zhang
    •  & Dongwu Shi

36. Guangzhou Medical University Second Affiliated Hospital, Guangzhou, China

    • Yichao Wen

37. Second Affiliated Hospital of Dalian Medical University, Dalian, China

    • Jian Yu

38. Zhengzhou University First Affiliated Hospital, Zhengzhou, China

    • Rongqing Sun

39. First People’s Hospital of Kunming, Kunming, China

    • Zhiwei Li

40. Union Hospital Affiliated to Tongji Medical College of Huanzhong University of Science and Technology, Wuhan, China

    • Shiying Yuan

41. Xinjiang Medical University Affiliated First Hospital, Wulumuqi, China

    • Yunlin Song

42. Chengdu University of Traditional Chinese Medicine Affiliated Hospital, Chengdu, China

    • Peiyang Gao

43. First Affiliated Hospital of Anhui Medical University, Hefei, China

    • Haiyan Liu

44. Yichang Central People’s Hospital, Yichang, China

    • Zhaohui Zhang

45. Suzhou Municipal Hospital, Suzhou, China

    • Yunfu Wu

46. Jining Medical College Affiliated Hospital, Jining, China

    • Biao Ma

47. First Affiliated Hospital of Soochow University, Suzhou, China

    • Qiang Guo

48. Qindao University Medical College Affiliated Hospital, Qindao, China

    • Feng Shan

49. Central South University Third Xiangya Hospital, Changsha, China

    • Mingshi Yang

50. 401 Military Hospital of China, Qindao, China

    • Hailing Li

51. Changsha Central Hospital, Changsha, China

    • Yuanfei Li

52. Yijishan Hospital of Wannan Medical College, Wuhu, China

    • Weihua Lu

53. Shanxi Medical University First Affiliated Hospital, Taiyuan, China

    • Lei Wang

54. Kuming Medical University First Affiliated Hospital, Kuming, China

    • Chuangyun Qian

55. Hebei Medical University Third Affiliated Hospital, Shijiazhuang, China

    • Zhiyong Wang

56. First Affiliated Hospital of Fujian Medical University, Fuzhou, China

    • Jiandong Lin

57. Zibo Central Hospital, Zibo, China

    • Rumin Zhang

58. First People’s Hospital of Yichang, Yichang, China

    • Peng Wan

59. Wuhan University Zhongnan Hospital, Wuhan, China

    • Zhiyong Peng

60. Wenzhou Medical University Second Affiliated Hospital, Wenzhou, China

    • Yuqiang Gong

61. Shantou University Medical College First Affiliated Hospital, Shantou, China

    • Linxi Huang

62. Zhongnan University Xiangya Second Hospital, Changsha, China

    • Guobao Wu

63. Yunnan Second People’s Hospital, Kunming, China

    • Jie Sun

64. Yancheng First People’s Hospital, Yancheng, China

    • Yijun Deng

65. First People’s Hospital of Foshan, Foshan, China

    • Lixin Zhou

66. Chongqing Medical University First Affiliated Hospital, Chongqing, China

    • Fachun Zhou

67. Xi’an Jiao Tong University First Affiliated Hospital, Xi’an, China

    • Qindong Shi

68. Beijing Wujing General Hospital, Beijing, China

    • Xiaodong Guo

69. Shenzhen People’s Hospital, Shenzhen, China

    • Xueyan Liu

70. Shanxi Dayi Hospital of Shanxi Academy of Medical Science, Taiyuan, China

    • Weidong Wu

71. 254 Military Hospital of China, Tianjin, China

    • Xiangzhong Meng

72. Qingdao University Affiliated Hospital, Qingdao, China

    • Liandi Li

73. Linhai First People’s Hospital, Linhai, China

    • Weiwei Chen

74. Tongji Hospital of Tongji Medical College of Huazhong University of Science and Technology, Wuhan, China

    • Shusheng Li

75. First Affiliated Hospital of Dalian Medical University, Dalian, China

    • Xianyao Wan

76. Beijing Chaoyang Hospital, Beijing, China

    • Zhixin Chao

77. Chongqing Medical University Second Affiliated Hospital, Chongqing, China

    • An Zhang

78. People’s Hospital of Yuxi City, Yuxi, China

    • Liming Gu

79. Shijitan Hospital of Capital Medical University, Beijing, China

    • Wei Chen

80. Shenzhen Nanshan People’s Hospital, Shenzhen, China

    • Jinglan Wu

81. Affiliated Hospital of Inner Mongolia Medical College, Huhehaote, China

    • Lihua Zhou

82. Benq Hospital, Nanjing, China

    • Zhenhuan Zhang

83. Luhe Hospital, Beijing, China

    • Yibing Weng

84. Beijing Jingmei Group General Hospital, Beijing, China

    • Yongshun Feng

85. Jiangxi Provincial People’s Hospital, Nanchang, China

    • Chunli Yang

86. Jinan University First Affiliated Hospital, Jinan, China

    • Yongjian Feng

87. General Hospital of Rocket Army, Beijing, China

    • Sumin Zhao

88. Hebei Medical University Second Affiliated Hospital, Shijiazhuang, China

    • Fei Tong

89. Binzhou Medical College Affiliated Hospital, Binzhou, China

    • Dong Hao

90. Chinese PLA General Hospital, Beijing, China

    • Hui Han

91. Yantai Mountain Hospital, Yantai, China

    • Baocai Fu

92. Tianjing Hospital of ITCWM Nankai Hospital, Tianjing, China

    • Chuanyong Gong

93. Hunan Provincial People’s Hospital, Changsha, China

    • Zhiping Li

94. People’s Hospital of Guangxi Zhuang Autonomous Region, Nanning, China

    • Kunlin Hu

95. Sun Yat-sen University Sixth Affiliated Hospital, Guangzhou, China

    • Qiuye Kou

96. China Academy of Chinese Medical Sciences Xiyuan Hospital, Beijing, China

    • Han Zhang

97. Wuhan General Hospital of Guangzhou Military Region, Wuhan, China

    • Jie Liu

98. First People’s Hospital of Yunnan, Kunming, China

    • Chuming Fan

99. Xinjiang Military General Hospital, Wulumuqi, China

    • Xin Zhou

100. Miyun Hospital of Beijing, Beijing, China

     • Xiumei Chen

101. Luoyang Central Hospital, Luoyang, China

     • Junli Sun

102. Huairou First Hospital of Beijing, Beijing, China

     • Xuejun Zhou

103. Military General Hospital of Beijing PLA, Beijing, China

     • Bin Song

104. Guangdong Provincial People’s Hospital, Guangzhou, China

     • Cheng Sun

105. Guangdong Second TCM Hospital, Guangzhou, China

     • Liyun Zhao

106. Beijing Hospital of TCM, Beijing, China

     • Xinglu Dong

107. Anhui Provincial Hospital, Hefei, China

     • Linlin Zhang

108. Shenyang First People’s Hospital, Shenyang, China

     • Dafei Tong

109. Guangzhou Military General Hospital, Guangzhou, China

     • Zhiguo Pan

110. Sun Yat-sen University First Affiliated Hospital, Guangzhou, China

     • Chuangjie Cai

111. Tianjing Cancer Hospital, Tianjing, China

     • Donghao Wang

112. Shanxi Cancer Hospital, Taiyuan, China

     • Yingjun Dong

113. Nanchang University Second Affiliated Hospital, Nanchang, China

     • Yuanqi Gong

114. Beijing University of Chinese Medicine Affiliated Dongfang Hospital, Beijing, China

     • Zhisong Wu

115. Shenzhen Second People’s Hospital, Shenzhen, China

     • Xinke Meng

116. Chendu Fifth People’s Hospital, Chendu, China

     • Ping Wang

Contributions

JX and ZZ participated in the study design and statistical analysis, and wrote the manuscript. WqL contributed to the study design, statistical analysis plan and critical revision of the manuscript. LK and JZ helped in design and editing of the manuscript. JX and the remaining authors not already specified participated in the data collection. JX helped with quality checking and authors’ connection. All authors approved the final manuscript.

Corresponding author

Correspondence to Weiqin Li.

Ethics approval and consent to participate

This study was approved by the medical ethics committee of Nanjing General Hospital of Nanjing Military Command (2017NZKY-010-01).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Reprints and Permissions