Patients with end-stage renal disease often have extensive CVD, with estimates of annual mortality as high as 20 %(1). This ongoing high mortality rate serves as an incentive to seek and establish more preventive strategies for dialysis-dependent patients.
Maintenance haemodialysis (MHD) is the use of a machine to filter wastes, salts and fluid from blood for at least 3 months to prolong the life of patients with advanced kidney failure. Several previous studies have found that MHD patients tend to have low dietary energy intake (DEI) and to be underweight(Reference Dwyer, Cunniff and Maroni2–Reference Therrien, Byham-Gray and Beto5). Inadequate energy intake is a frequent and important cause of protein-energy wasting, which closely related to major adverse clinical outcomes in MHD patients(Reference Ikizler, Cano and Franch6,Reference Chen, Qin and Li7) . Hence, ensuring the optimal amount of energy intake is very important to improve the clinical outcomes of MHD patients. A metabolic balance study on six MHD patients suggested that a mean energy intake of 146·4 kJ/kg per d was necessary to maintain both a neutral nitrogen balance and an unchanged body composition(Reference Slomowitz, Monteon and Grosvenor8). Mainly based on this study, the European Best Practice Guidelines on nutrition recommend a DEI of 125·5–167·4 kJ/kg of ideal body weight (IBW) per d in clinically stable MHD patients(Reference Fouque, Vennegoor and ter Wee9). The guidelines from the Kidney Disease Outcome Quality Initiative recommend a daily DEI of 146·4 kJ/kg body weight per d for MHD patients who are <60 years of age and a DEI of 125·5–146·4 kJ/kg body weight per d for those 60 years of age or older(4). And, the International Society of Renal Nutrition and Metabolism consensus indicates that the minimum energy requirements for MHD patients are 125·5–146·4 kJ/kg of IBW per d based on the physical activity level(Reference Ikizler, Cano and Franch6).
However, few epidemiological studies have prospectively examined the association between DEI and risk of mortality in MHD patients, and the optimal energy intake for prevention of mortality in such patients is still unknown. Estimation of DEI can only be done by monitoring intake. The most common methods for estimating dietary intake are FFQ, food intake records and 24-h dietary recalls(Reference Bross, Noori and Kovesdy10,Reference Hebert, Hurley and Chiriboga11) . Given evidence from biomarker-based validation, the 24-h dietary recalls capture dietary intake with less bias than do FFQ(Reference Subar, Kipnis and Troiano12,Reference Prentice, Mossavar-Rahmani and Huang13) . Moreover, the US Department of Agriculture five-step multiple-pass method for dietary recall has been validated by a lot of previous studies(Reference Burrows, Martin and Collins14–Reference Conway, Ingwersen and Moshfegh17). Therefore, we make the assumption that the 24-h recalls used in the present study reflected the DEI because of earlier validations. Our purpose was to explore the relationship of DEI with all-cause and CVD mortality and to examine any possible effect that modifies among MHD patients.
Methods
Study design and participants
The present study was a multi-centre, prospective cohort study conducted from January 2014 to December 2015, in eight outpatient dialysis centres (including Nanfang Hospital, the First People’s Hospital of Foshan, Huadu District People’s Hospital of Guangzhou, Guangzhou Red Cross Hospital, Guangzhou Overseas Chinese Hospital, the Third Affiliated Hospital of Southern Medical University, Nanhai District People’s Hospital of Foshan and the Second People’s Hospital of Shenzhen) in Guangdong province in China. Eligible participants were males and females aged more than 18 years who had received maintenance dialysis for at least 3 months and had a normal oral intake. We excluded those who had hyperthyroidism, acute infection, liver cirrhosis, multiple organ failure, serious gastrointestinal disease, cognitive disorder and advanced malignant tumour.
Participants were followed up at each routine dialysis visit where vital signs and possible endpoint events were documented by trained research staff and physicians.
The study was approved by the Medical Ethics Committee of Nanfang Hospital and all the participants signed an informed consent.
Data collection and measurements
Baseline data were collected by trained research staff according to standard operating procedures. Interviews were conducted on each participant using a standardised questionnaire designed specifically for the present study.
Diabetes mellitus was defined as having a medical history of diabetes, or undergoing glucose-lowing therapy, including insulin or oral hypoglycaemic agents. Patients who reported antihypertensive drugs use or had a medical history of hypertension were defined as having hypertension. The history of CVD was defined as history of angina, congestive heart failure, transient ischemic attack, history of myocardial infarction or cerebrovascular accident and peripheral arterial disease. The question about physical activity was phrased as follows, ‘How do you describe your daily physical activity level?’ and a choice of three responses: low, moderate and high was provided(Reference Wang, Li and Qin18,Reference Qin, Li and Zhang19) .
Anthropometric measures including weight and height were obtained during the physical examination. All measurements were conducted post-dialysis when the patient was at dry weight. BMI was calculated as weight/height squared (kg/m2).
Dietary intake was assessed by trained interviewers using the automated multiple-pass method where 24-h diet recalls were administered on 3 d (including one dialysis day and two non-dialysis days) within 1 week. The 24-h diet recalls are relatively quick assessment modalities to obtain the most recent information about food intake and are widely used in epidemiological studies(Reference Bross, Noori and Kovesdy10). Briefly, the automated multiple-pass method was developed by the US Department of Agriculture and utilises a five-step, multiple-pass, interviewer-administered, computerised process for obtaining diet recalls. The automated multiple-pass method consists of five following steps: (a) the quick list, which is an uninterrupted listing by the subject of foods and beverages consumed; (b) the forgotten foods list, which queries the subject on categories of foods that have been documented as frequently forgotten; (c) a time and occasion at which foods were consumed; (d) the detail cycle, which elicits descriptions of foods and amounts eaten aided by the use of the Food Model Booklet; and finally (e) the final probe review(Reference Subar, Kipnis and Troiano12,Reference Conway, Ingwersen and Moshfegh17) . The energy intake of all food and drink items was sourced using a computer-aided dietary software (China software version 2.0), in which nutrient models were based on the Chinese Food Composition Table developed by the Chinese Center for Disease Control and Prevention, 2009.
DEI was expressed in three units of measure: actual energy intake in kJ/d, energy intake normalised to actual body weight (ABW) in kJ/kg ABW per d and energy intake normalised to IBW in kJ/kg IBW per d. As guidelines recommended, we used IBW instead of ABW, since normalising energy intake to ABW would result in erroneously high- or low-energy requirements in overweight or underweight patients(Reference Ikizler, Cano and Franch6,Reference Fouque, Vennegoor and ter Wee9) . In the main analysis of the present study, IBW was calculated using improved Broca formula adapted for Asian populations: IBW = (height (cm) – 100) × 0·9 (kg) for males and (height (cm) – 100) × 0·9–2·5 (kg) for females(Reference Shah, Sucher and Hollenbeck20,Reference Rao and Balakrishna21) . Moreover, IBW was also calculated by multiplying an ideal BMI of 22 kg/m2 with a person’s actual height (m) squared(Reference Shah, Sucher and Hollenbeck20,Reference Watanabe, Machida and Matsumoto22) in the sensitivity analysis.
Blood samples of all participants were obtained from each subject prior to the haemodialysis session at baseline. Serum lipids, albumin, Ca, phosphate levels and total CO2 were measured using automatic clinical analysers following the same standard protocol at each local dialysis centre.
Study outcomes
All-cause mortality was the primary outcome in this analysis and included death due to any reason. The secondary outcome was death from CVD, which included sudden cardiac death, myocardial infarction, heart failure, stroke, cardiovascular haemorrhage and death due to other known vascular causes. Evidence for death included death certificates from hospitals or reports from investigator visits.
Statistical analysis
We assumed that the annual mortality rate of MHD patients is 20 %(1), the hazard ratio (HR) of the DEI between <104·6 and ≥104·6 kJ/kg IBW per d was about 1·67(Reference Araújo, Kamimura and Draibe23), with an α of 0·05, and 124 MHD patients followed up for 2 years would have a power of 90 %.
Baseline characteristics are presented as means and standard deviations for normally distributed continuous variables, medians (Q1–Q3) for skewed continuous variables and frequencies (percentages) for categorical variables. Differences in baseline characteristics according to DEI categories (<125·5, 125·5–146·4, 146·4–167·4 and ≥167·4 kJ/kg IBW per d) were compared using χ 2 tests or ANOVA tests as appropriate.
Cox proportional hazards models were used to estimate the HR and 95 % CI for the risk of all-cause mortality and CVD mortality associated with DEI, without and with adjustments for dialysis centre, age, sex, smoking, alcohol use, physical activity, BMI, albumin, TAG, total cholesterol (TC), C-reactive protein (CRP), Ca, phosphate, Hb, dialysis vintage, Kt/V, diabetes, hypertension and history of CVD. We applied two piecewise regression models to examine the threshold effect of DEI on the risk of all-cause mortality or CVD mortality using a smoothing function. The threshold level (turning point) was determined using likelihood-ratio tests and bootstrap re-sampling methods.
In a stratified analysis, possible modifications of the association between DEI and all-cause mortality were assessed for the variables, age (<60 v. ≥60 years), sex, diabetes (yes v. no), history of CVD (yes v. no), CRP (<8 v. ≥8 mg/l)(Reference Zimmermann, Herrlinger and Pruy24), TC (<5·2 v. ≥5·2 mmol/l)(Reference Sarnak, Bloom and Muntner25) and BMI (<23 v. ≥23 kg/m2)(26) levels at baseline.
A two-tailed P < 0·05 was considered statistically significant in all analyses. All analyses were performed using Empower (R) (www.empowerstats.com, X&Y Solutions, Inc.) and R (http://www.R-project.org).
Results
Study participants and baseline characteristics
As illustrated in the flow chart (online Supplementary Fig. S1), 1039 MHD patients were included in the final analyses.
At baseline, 577 (55·5 %), 222 (21·4 %), 154 (14·8 %) and 86 (8·3 %) participants had DEI <125·5, 125·5 to <146·4, 146·4 to <167·4 and ≥167·4 kJ/kg IBW per d (<30, 30–<35, 35–<40 and ≥40 kcal/kg IBW per d), respectively. Mean values of DEI in total population were 6695·7 (sd 1860·2) kJ per d and 122·6 (sd 33·1) kJ/kg IBW per d, respectively.
Participant characteristics by DEI categories are listed in Table 1. Patients with lower DEI tended to be older, had lower TC, serum Ca, serum phosphate, CRP, dietary protein intake, and were more likely to be males.
Table 1. Baseline characteristics of the maintenance haemodialysis patients*
(Means and standard deviations; medians and interquartile ranges (IQR); numbers and percentages)
DEI, dietary energy intake; IBW, ideal body weight; TC, total cholesterol; CRP, C-reactive protein; ABW, actual body weight; DPI, dietary protein intake.
* Continuous variables are expressed as means and standard deviations or medians and 25th–75th percentiles; and categorical variables are given as n and %.
† To convert kcal to kJ, multiply by 4·184.
‡ Kt/V, Kt showed effective urea clearance and duration of dialysis, and V represents the volume of distribution of urea in the body, calculated as Kt/V = −ln (post-BUN/pre-BUN – 0·008 × t) + (4 – 3·5 × post-BUN/pre-BUN) × UF/post-weight, where t is the effective dialysis time, BUN is serum blood urea nitrogen and UF is ultrafiltration.
§ Ideal body weight was defined as a BMI of 22 kg/m2.
Relationship of dietary energy intake with all-cause and CVD mortality
Over a median follow-up of 28 months, eighty-four patients underwent kidney transplantation, five patients switched to peritoneal dialysis and fifty-five patients were lost to follow-up. These patients were censored at the last date of follow-up. All-cause and CVD mortality occurred in 230 (22·1 %) and 140 (13·5 %) participants, respectively.
Overall, a U-shaped relationship was observed between DEI and all-cause mortality (Fig. 1(a)). The two piecewise regression models showed that a DEI value of 167·4 kJ/kg IBW per d (40 kcal/kg IBW per d) was a turning point for all-cause mortality after multi-variate adjustment. That is, the risk of all-cause mortality decreased significantly with the increase of DEI in participants with DEI <167·4 kJ/kg IBW per d (40 kcal/kg IBW per d) (HR 0·98; 95 % CI 0·96, 1·00; P = 0·027), but increased significantly with the increase of DEI in participants with DEI ≥167·4 kJ/kg IBW per d (40 kcal/kg IBW per d) (HR 1·12; 95 % CI 1·04, 1·20; P = 0·002) (Table 2). Consistently, when compared with participants with DEI of 146·4 to <167·4 kJ/kg IBW per d (35 to <40 kcal/kg IBW per d), a significantly higher risk of all-cause mortality was observed in those with DEI <125·5 kJ/kg IBW per d (30 kcal/kg IBW per d) (HR 1·85; 95 % CI 1·20, 2·85), and a higher but insignificant risk of all-cause mortality was also found in participants with DEI ≥167·4 kJ/kg IBW per d (40 kcal/kg IBW per d) (HR 1·83; 95 % CI 0·99, 3·40) (Table 3).
Fig. 1. Association between dietary energy intake (DEI) normalised to ideal body weight (IBW) and risk for all-cause mortality (a) and CVD-related mortality (b)*. HR, hazard ratio. * Adjusted for dialysis centre, age, sex, smoking, alcohol use, physical activity, BMI, albumin, TAG, total cholesterol, C-reactive protein, calcium, phosphate, Hb, dialysis vintage, Kt/V, systolic blood pressure, diabetes and history of CVD. † To convert kcal to kJ, multiply by 4·184.
Table 2. Threshold effect analyses of dietary energy intake (DEI) levels on the risk of all-cause mortality and CVD mortality using two piecewise regression models
(Numbers and percentages; hazard ratios (HR) and 95 % confidence intervals)
IBW, ideal body weight.
* To convert kcal to kJ, multiply by 4·184.
† Adjusted for dialysis centre, age, sex, smoking, alcohol use, physical activity, BMI, albumin, TAG, total cholesterol, C-reactive protein, Ca, phosphate, Hb, dialysis vintage, Kt/V, systolic blood pressure, diabetes and history of CVD.
Table 3. Univariate and multi-variate Cox regression models of all-cause mortality and CVD-related mortality for dietary energy intake (DEI) groups
(Numbers and percentages; hazard ratios (HR) and 95 % confidence intervals)
IBW, ideal body weight; Ref, reference.
* To convert kcal to kJ, multiply by 4·184.
† Adjusted for dialysis centre, age, sex, smoking, alcohol use, physical activity, BMI, albumin, TAG, total cholesterol, C-reactive protein, Ca, phosphate, Hb, dialysis vintage, Kt/V, systolic blood pressure, diabetes and history of CVD.
Similarly, a U-shaped association was also found between DEI and CVD mortality (Fig. 1(b)). The two piecewise regression models showed that a DEI value of 152·7 kJ/kg IBW per d (36·5 kcal/kg IBW per d) was a turning point for CVD mortality after multi-variate adjustment. The risk of CVD mortality decreased significantly with the increase of DEI in participants with DEI <152·7 kJ/kg IBW per d (36·5 kcal/kg IBW per d) (HR 0·96; 95 % CI 0·93, 0·99; P = 0·006) but increased significantly with the increase of DEI in participants with DEI ≥ 152·7 kJ/kg IBW per d (36·5 kcal/kg IBW per d) (HR 1·11; 95 % CI 1·04, 1·18; P = 0·002) (Table 2). Consistently, when compared with participants with DEI of 146·4 to <167·4 kJ/kg IBW per d (35 to <40 kcal/kg IBW per d), a significantly higher risk of CVD mortality was observed in those with DEI <125·5 kJ/kg IBW per d (30 kcal/kg IBW per d) (HR 1·94; 95 % CI 1·12, 3·36) or ≥167·4 kJ/kg IBW per d (40 kcal/kg IBW per d) (HR 2·39; 95 % CI 1·14, 5·03) (Table 3).
There was a similar trend in the relationship between DEI by ABW and all-cause or CVD mortality (online Supplementary Fig. S2). Moreover, DEI normalised to IBW calculated from BMI showed the similar relationship with mortality (online Supplementary Fig. S3). Further adjustment for DPI did not substantially change the results (online Supplementary Table S1).
Stratified analyses by potential effect modifiers
Stratified analyses were performed by age, sex, diabetes, history of CVD, CRP, TC and BMI. The lowest risk of mortality was found in patients with DEI at 125·5 to <167·4 kJ/kg IBW per d (30 to <40 kcal/kg IBW per d) in all the sub-groups. There were no significant interactions in any of the sub-groups (P for interaction >0·05 for all comparisons; Fig. 2).
Fig. 2. Sub-group analysis of the relationship of dietary energy intake (DEI) ((a) <30 v. 30 to <40 kcal/kg IBW per d and (b) ≥40 v. 30 to <40 kcal/kg IBW per d) with all-cause mortality*. DM, diabetes mellitus; HR, hazard ratio; IBW, ideal body weight; TC, total cholesterol. * Adjusted for dialysis centre, age, sex, smoking, alcohol use, physical activity, BMI, albumin, TAG, total cholesterol, C-reactive protein, calcium, phosphate, Hb, dialysis vintage, Kt/V, systolic blood pressure, diabetes and history of CVD. To convert kcal to kJ, multiply by 4·184. To convert cholesterol in mg/dl to mmol/l, multiply by 0·0259.
Discussion
Our study is the first to demonstrate a U-shaped association between DEI and all-cause or CVD mortality in MHD patients. Participants with a DEI at 146·4 to <167·4 kJ/kg IBW per d (35 to <40 kcal/kg IBW per d) had the lowest risk of all-cause and CVD mortality.
Haemodialysis patients are at increased risk of energy deficiency for several reasons, including inadequate energy intake, inflammation and multiple co-morbidities, metabolic and hormonal derangements, haemodialysis-associated catabolism and increased resting energy expenditure(Reference Ikizler, Cano and Franch6,Reference Ikizler, Wingard and Sun27–Reference Pupim and Ikizler30) . A previous retrospective study including 344 incident haemodialysis patients found that energy intake lower than 104·6 kJ/kg per d was related to worse survival(Reference Araújo, Kamimura and Draibe23). Another retrospective observational study also reported that lower energy intake (<104·6 kJ/kg per d) was associated with all-cause mortality among 144 patients on MHD(Reference Kang, Chang and Park31). Our findings are in agreement with these last two retrospective studies with small sample size, for we also found that DEI levels were inversely associated with all-cause mortality in participants with DEI <167·4 kJ/kg IBW per d (40 kcal/kg IBW per d). The mechanisms linking lower energy intake and the risk of mortality are still unclear. When DEI is deficient, body fat and muscle mass are depleted to maintain BMR and protein stores cannot be able to well preserved(Reference Ikizler, Cano and Franch6,Reference Fouque, Kalantar-Zadeh and Kopple32) . In addition, pro-inflammatory cytokines, such as IL-6 or TNF-α, are higher in MHD patients with reduced nutrient intake(Reference Qureshi, Alvestrand and Divino-Filho33–Reference Kalantar-Zadeh, Block and McAllister36). Loss of muscle and fat stores and inflammation may increase the risk of all-cause mortality. Moreover, patients with low energy intake might also have a higher prevalence of latent chronic disease(Reference Leosdottir, Nilsson and Nilsson37).
Our study had a relatively higher DEI compared with some studies. In our study, the mean and standard deviation DEI were 121·8 (sd 34·7) kJ/kg of ABW per d and 122·6 (sd 33·1) kJ/kg IBW per d. In comparison, the mean DEI by ABW in haemodialysis patients was 95·8 (sd 35·1) kJ/kg per d in the Haemodialysis (HEMO) Study cohort(Reference Rocco, Paranandi and Burrowes38), 107·9 kJ/kg per d in a Korean cohort(Reference Kang, Chang and Park31) and 109·6 kJ/kg per d in a Brazilian cohort(Reference Araújo, Kamimura and Draibe23). The differences in mean DEI across these different studies may be due to several reasons. First, dietary habits vary greatly across different countries and cultures. Second, the nutritional status of a population is known to improve with its economic development; the living standards for these populations may have improved considerably since each study was conducted. The HEMO study(Reference Rocco, Paranandi and Burrowes38) was conducted between May 1995 and March 1997, the study from Korea(Reference Kang, Chang and Park31) was conducted during April 2006 and the study from Brazil(Reference Araújo, Kamimura and Draibe23) was conducted from January 1992 to December 2002. By contrast, our study was conducted more recently, from January 2014 to December 2015. Third, each study used a different method to assess DEI thereby potentially over or under estimating the true energy intake. DEI in the HEMO study was assessed using a 2-d dietary intake diary on one dialysis day and one non-dialysis day and was recorded by patients. However, in our study, DEI was assessed through face-to-face or telephone interview by trained staff on two non-dialysis days and one dialysis day. Fourth, at each follow-up visit, all patients in the present study received information on how to improve their protein and energy intake and may be a major contributor to the relatively high DEI observed in our study. Still, some studies observed similar results compared with our study. In a French study, Chauveau et al. found results with a DEI of 124·7 kJ/kg per d(Reference Chauveau, Grigaut and Kolko39). On the other hand, Arslan and Kiziltan assessed a higher DEI of 143·1 kJ/kg per d in ninety-three MHD patients(Reference Arslan and Kiziltan40).
The relatively high DEI level and large sample size of the present study offer an exceptional opportunity to examine the relationship of higher DEI levels with all-cause and CVD mortality in MHD patients. Our study further indicated that all-cause mortality risk increased significantly with the increase of DEI in participants with DEI ≥167·4 kJ/kg IBW per d (40 kcal/kg IBW per d). Consistently, previous studies reported that high energy intake was associated with increased risks of cancer, incident diabetes, all-cause mortality and CVD mortality in general populations(Reference Silvera, Jain and Howe41–Reference Hingle, Wertheim and Neuhouser43). Iff et al.(Reference Iff, Wong and Webster44) found that higher relative energy intake increases mortality in patients with estimated glomerular filtration rate <60 ml/min per 1·73 m2. The biological mechanisms linking higher DEI to mortality risk are not well known. It has been suggested that excessive energy exposure may increase oxidative stress, leading to diminished mitochondrial function, decreased ATP productive capacity, cell death and multiple organ failure(Reference Cano, Fouque and Roth45). However, more studies are needed to further investigate our findings and examine the underlying mechanisms.
Due to the exclusion of patients with acute infection, CRP was at a relatively low level in the present study. Furthermore, our study showed that patients with higher DEI tended to had higher CRP levels in the present study. At same time, participants with higher DEI also had higher TC levels and there was a positive association between CRP and TC (P = 0·038). It has been reported that cholesterol accumulation in macrophages and other immune cells may promote inflammatory responses, including augmentation of Toll-like receptor signalling, inflammasome activation and the production of monocytes and neutrophils(Reference Tall and Yvan-Charvet46). Although the underlying mechanisms still need to be further investigated in future studies, we speculated that elevated TC associated with higher DEI may partly explain the increased CRP levels in this population. More importantly, the stratified analysis indicated that both TC and CRP levels did not significantly modify the association between DEI and the mortality risk.
Our findings should be interpreted in the context of several limitations. First, the DEI was assessed using the 24-h diet recalls, which was dependent on reporting by patients(Reference Bross, Noori and Kovesdy10). However, it is one of the most common methods for obtaining dietary intake data as there are no alternatives to calculate accurate nutrient intake(4). To minimise the potential errors, we used extensive quality control procedures and three 24-h diet recalls including one dialysis day and two non-dialysis days to yield more accurate data. More importantly, the automated multiple-pass method approach to 24-h diet recalls that was used in the present study has been validated by a lot of previous studies(Reference Burrows, Martin and Collins14,Reference Bingham, Gill and Welch16,Reference Conway, Ingwersen and Moshfegh17,Reference Johnson, Driscoll and Goran47,Reference Sawaya, Tucker and Tsay48) . Although the determination of DEI may possible be unreliable in both obese and very thin(Reference Poslusna, Ruprich and de Vries49), the exclusion of patients with BMI < 18·5 or ≥30 kg/m2 did not substantially change our results (online Supplementary Table S2). Second, dietary intake was assessed at only a single time which did not take into consideration the changing pattern of eating and physical behaviour over the follow-up period. Third, this was an observational study. Despite extensive adjustment for known confounders in our analysis, we cannot exclude the possibility that unrecorded factors may explain some of our findings. Fourth, the study was carried out in China. Whether the results can be extrapolated to other populations requires further verification. Therefore, further confirmation of our findings in an independent trial is greatly needed.
In conclusion, our study suggests a U-shaped association between DEI and all-cause or CVD mortality, with a turning point at about 167·4 and 152·7 kJ/kg IBW per d (40 and 36·5 kcal/kg IBW per d), respectively, in MHD patients. Our findings provide new insights with respect to the benefit:risk ratio of energy supplementation in MHD patients.
Acknowledgements
The authors thank the participants, investigators and staff for their contribution to the present study.
The present study was supported by the National Key Technology Support Program of China (grant no. 2015BAI12B00) and the High-level Matching Funds of Nanfang Hospital (2014070).
Dr X. Q. reports grants from the National Natural Science Foundation of China (81730019) and Outstanding Youths Development Scheme of Nanfang Hospital, Southern Medical University (2017J009). Dr F. H. reports grants from the Science and Technology Planning Project of Guangzhou and the National Key Research and Development Program.
Research idea and study design: M. L., S. L. Y., Y. M. L. and Y. Y. Y.; data acquisition: S. L. Y., Y. M. L., Y. Y. Y., Y. L., Z. H. L., Y. Z. K., Y. L., Q. W., F. N. L., S. H., A. Q. L. and Q. J. W.; data analysis/interpretation: Y. M. L, Y. Y. Y., Y. Z. K., Y. L., Q. W., F. N. L., S. H., A. Q. L. and Q. J. W.; statistical analysis: Y. Y. Y. and X. H. Q.; supervision or mentorship: M. L. and X. H. Q. Each author contributed important intellectual content during manuscript drafting or revision and accepts accountability for the overall work by ensuring that questions pertaining to the accuracy or integrity of any portion of the work are appropriately investigated and resolved.
There were no conflicts of interest.
Supplementary material
For supplementary material referred to in this article, please visit https://doi.org/10.1017/S0007114519002861
