Skip to main content

Influences of changes in obesity-related parameters on progression of arteriosclerosis in hemodialysis patients: a prospective cohort study

Abstract

Background

Hemodialysis patients have a greater risk of cardiovascular disease, their major cause of death. Obesity is a risk factor for arteriosclerosis, cardiovascular disease, and mortality in the general population. Conversely, an increased body mass index (BMI) is known to be associated with better survival of hemodialysis patients. However, the influences of changes in parameters associated with obesity on arteriosclerosis in these patients have not been evaluated. Here, we evaluated the relationship between changes in parameters associated with obesity and brachial-ankle pulse wave velocity (baPWV), as an index of arteriosclerosis, in hemodialysis patients.

Methods

Two hundred and fifty-eight subjects were selected from patients undergoing stable hemodialysis. BMI, abdominal circumference (AC), visceral fat area (VFA) and subcutaneous fat areas (SFA), and baPWV were assessed at baseline and 4 years later. Ankle-brachial index and baPWV were measured on a non-dialysis day using volume-plethysmographic apparatus. Abdominal circumference was measured at the level of umbilicus by computed tomography (CT). Body fat distribution was determined using CT. All CT scans were performed on a non-dialysis day with the subject in a supine position.

Results

Sixty-one patients died during the 4-year period. Abdominal circumference, BMI, VFA, and SFA were higher, and baPWV was lower in patients who survived than in those who died. During the study period, baPWV increased and was positively correlated with the change in visceral plus subcutaneous fat area in patients whose initial VFA was > 100 cm2. There were negative correlations between the change in baPWV and changes in AC, and SFA in individuals whose VFA was < 100 cm2.

Conclusions

These data suggest that patients with a large VFA should adopt strategies aimed at reducing body weight, while patients with a small VFA should avoid undernutrition.

Background

Patients undergoing hemodialysis (HD) have a greater risk of arteriosclerotic complications and a poor prognosis due to the greater prevalence of cardiovascular (CV) diseases. The risk factors include hypertension, diabetes mellitus, vascular calcification, and inflammation. Obesity has been recognized as a risk factor for arteriosclerosis and higher mortality and morbidity from CV disease in the general population. In addition, a positive relationship between visceral fat accumulation and arteriosclerosis has been reported in diabetic HD patients [1]. In contrast, an inverse association has been reported in female HD patients [2]. It has been proposed that the relationships between obesity-associated parameters and arteriosclerosis are complicated by the existence of malnutrition-inflammation-atherosclerosis (MIA) syndrome in these patients [3], and the influences of changes in obesity-related parameters on arteriosclerosis in patients undergoing long-term HD have not been evaluated.

Brachial-ankle pulse wave velocity (baPWV), which is an easily performed, non-invasive measurement, is a marker of arteriosclerosis of both central elastic and peripheral muscular components of the circulation [4]. In the general population, high baPWV is an independent predictor of mortality [5], and we have shown that baPWV is a predictor of CV outcome in patients undergoing HD [6], indicating that baPWV is a useful index of arteriosclerosis in HD patients.

Therefore, in this study, we aimed to determine the relationship between changes in parameters associated with obesity and baPWV in HD patients.

Methods

Patient recruitment

This study was approved by the ethics committee of Moriguchi Keijinkai Hospital (approval no.: 3003). All participants were enrolled after providing written informed consent. The participants were outpatients undergoing maintenance dialysis at Kadoma Keijinkai Clinic in Kadoma City, Neyagawa Keijinkai Clinic in Neyagawa City, or Moriguchi Keijinkai Clinic in Moriguchi City, Osaka, Japan. Patients were selected from those who were undergoing regular dialysis treatment and were hemodynamically stable during dialysis. Patients with malignancy and amputation of both legs were excluded in this study. The enrollment period was from January 2013 to December 2013.

Patient characteristics

At enrollment, information was collected regarding the age, gender, BMI, duration of HD, presence of diabetes, and use of anti-hyperlipidemic drugs. Diabetes was diagnosed when the patients were already under treatment with an oral hypoglycemic agent or insulin, or met the World Health Organization criteria [7] in the past or present. Pre-dialysis systolic and diastolic blood pressures (BPs) and pre-dialysis values of hemoglobin, high-density lipoprotein cholesterol (HDL-c), low-density lipoprotein cholesterol (LDL-c), triglyceride, uric acid, creatinine, c-reactive protein (CRP), calcium, inorganic phosphorus, intact parathyroid hormone, albumin, and hemoglobin A1c were measured on the first dialysis day of the week. Post-dialysis value of human atrial natriuretic peptide (hANP), brain natriuretic peptide (BNP), and cardiothoracic ratio (CTR) were measured on the first dialysis day of the week.

Ankle-brachial index and brachial-ankle pulse wave velocity

Ankle-brachial index (ABI) and baPWV were measured on a non-dialysis day using a volume-plethysmographic apparatus (Omron Healthcare Co., Ltd., Kyoto, Japan), as previously described [8]. The higher of the values for each leg was accepted as the baPWV value. Because baPWV cannot be estimated accurately when ABI is < 0.9, due to a reduction induced by arterial occlusion [9], lower limbs with an ABI of < 0.9 were excluded from the analysis of baPWV.

Computed tomography

Abdominal circumference (AC) was measured at the level of the umbilicus by computed tomography (CT) using Ziostation 2 software (Ziosoft, Tokyo, Japan). Body fat distribution was determined using CT. All CT scans were performed on a non-dialysis day with the subject in a supine position, using an 80-slice multidetector CT scanner (Aquilion PRIME, Canon medical systems Corporation, Tochigi, Japan). Three-millimeter-thick slices were acquired. The subcutaneous fat layer was defined as the layer lying between the skin and abdominal muscle that showed attenuation from − 60 to − 160 HU. The component of the image with the same density as the subcutaneous fat layer in the abdominal cavity was defined as the visceral fat. The subcutaneous fat area (SFA) and VFA were measured in a slice at the level of the umbilicus.

The CT attenuation of adipose tissue is thought to be approximately − 60 to − 160 HU. However, when adipose tissue is edematous, this value is higher because of the presence of water (attenuation = 0). In such cases, accurate measurement of fat area by CT is impossible. Therefore, we did not include patients with gross edema in our analysis. In our study, the presence of gross edema was defined by the area with a radiodensity of − 60 to − 160 HU being smaller than that with a radiodensity of − 30 to − 50 HU at the level of the umbilicus. Consequently, 12 patients were not included in our analysis.

Study protocol

Patient characteristics were recorded at baseline, and the BMI, AC, VFA, SFA, and baPWV were measured at baseline and 4 years later. Comparison of patients’ characteristics between those who survived and those who died was done. Of the patients who survived, those with complete data at baseline and 4 years later were analyzed for the influences of changes in obesity-related factors on arteriosclerosis over a 4-year period.

Statistical analysis

Data are expressed as mean ± SD. Statistical comparisons between patients who survived and those who died were performed using Student’s t test or Fisher’s exact test. Paired t tests were used to assess changes over the 4-year period in patients who survived and had complete data at baseline and 4 years later. Pearson’s correlation coefficient was used to evaluate the relationships between the change in baPWV and changes in other obesity-related parameters over the period of the study. P < 0.05 was considered to represent statistical significance. All analyses were performed using EZR (Saitama Medical Center, Jichi Medical University, Saitama, Japan) [10].

Results

A total of 258 patients were enrolled in the study. Table 1 shows details of the characteristics of the study subjects. Out of 258 patients, 61 patients died within 4 years. The patients’ flowchart is presented in Fig. 1.

Table 1 Characteristics of the study subjects
Fig. 1
figure 1

Patients’ flowchart

Comparison of patients’ characteristics between those who survived and those who died is shown in Table 2. Age, presence of diabetes, hANP, BNP, CTR, and baPWV were significantly higher in patients who died than in those who survived. BMI, serum creatinine, serum albumin, AC, VFA, and SFA were significantly lower in patients who died than in those who survived.

Table 2 Comparison of characteristics of the patients who survived and those who died

Of the patients who survived, 102 patients with complete data were analyzed for the influences of changes in obesity-related factors on arteriosclerosis over a 4-year period. The characteristics of these 102 patients at baseline are shown in Table 3. The BMI (initial value 22.4 ± 3.0 kg/m2, final value 22.2 ± 3.3 kg/m2), AC (initial value 83.4 ± 8.6 cm, final value 82.6 ± 9.4 cm), VFA (initial value 86.2 ± 54.4 cm2, final value 85.7 ± 54.0 cm2), and SFA (initial value 120.1 ± 62.7 cm2, final value 117.9 ± 60.2 cm2) did not significantly change during the study period. baPWV significantly increased (initial value 1913 ± 446 cm/s, final value 2029 ± 546 cm/s; P < 0.01), while BPs did not significantly change (sBP initial value 145 ± 22 mmHg, final value 145 ± 21 mmHg; dBP initial value 79 ± 12 mmHg, final value 77 ± 13 mmHg).

Table 3 Characteristics of the surviving patients who were included in the final analysis

Of these 102 patients, 36 had an initial VFA of ≥ 100 cm2, which meets the Japanese obesity criteria [11]. In these patients, the change in baPWV over 4 years was significantly positively correlated with the changes of VFA + SFA (r = 0.33, P < 0.05) (Table 4). Sixty-six patients had an initial VFA of < 100 cm2. In these patients, the change in baPWV over 4 years was significantly negatively correlated with the changes in AC (r = − 0.333, P < 0.01), SFA (r = − 0.367, P < 0.01), and VFA + SFA (r = − 0.297, P < 0.05) (Table 4).

Table 4 Correlation with ΔbaPWV

Discussion

There were three main findings in this study. First, baPWV was higher and BMI, AC, VFA, and SFA were lower in the patients who died than in those who survived. Second, baPWV increased over the 4-year period of the study. Third, the change in baPWV was positively correlated with the changes in obesity-related parameters in patients whose initial VFA was ≥ 100 cm2 and was negatively correlated with the change in obesity-related parameters in those whose initial VFA was < 100 cm2. These findings suggest that increasing obesity in individuals with a large visceral fat area and decreasing adiposity in individuals with a low visceral fat area are both associated with the progression of arteriosclerosis.

Characteristics of patients who survived and those who died

There were significant differences in some of the characteristics of the patients who died and those who survived the 4-year study period. BMI, serum albumin, and adiposity were significantly lower among those who died. Thus, in this study, good nutritional status and slightly greater fat accumulation were beneficial for HD patients, trends which are consistent with previous reports [12, 13]. Higher AC was also beneficial in our study. The hANP and BNP are known to be independent risk factors for CV and overall mortality in HD patients [14, 15], and in our study, they were significantly higher in the patients who died. High baPWV is thought to be a risk factor for CV and overall mortality [5], and in terms of the mortality of HD patients, the same trends of high CTR [16] and concomitant diabetes [13] were identified in this study as in previous reports. In addition, the mortality rate during 4 years was 24.6% in this study, which is almost equivalent to that of Japanese dialysis patients whose crude death rate is about 10%/year [17]. On the basis of these consistent findings, the patients in our study appear to be representative HD patients.

Time-dependent change in brachial-ankle pulse wave velocity in hemodialysis patients

Patients undergoing HD have a greater risk of arteriosclerotic complications and a poor prognosis due to the greater prevalence of CV diseases [18,19,20]. Most HD patients are thought to have developed arteriosclerosis by the time they commence HD, and afterwards, arteriosclerosis is thought to be accelerated [21]. However, the natural course of PWV has not been comprehensively assessed in longitudinal studies. Avramovski et al. showed in a comparison of 80 HD patients and 60 controls that carotid-femoral PWV was higher in HD patients at baseline, but increased more in these patients over a 36-month period [22]. In the present study, we have shown that baPWV increased during a 4-year period in 102 HD patients. PWV has been reported to reflect not only arterial stiffness but also the BPs at the time of PW measurement [23]. In this study, however, the BPs did not show significant changes. Therefore, high baPWV is considered to reflect progression of arteriosclerosis, supporting the notion that the natural course of arteriosclerosis is accelerated in patients undergoing HD.

Influences of changes in obesity-related parameters on arteriosclerosis in hemodialysis patients

Accumulation of visceral fat is a risk factor for arteriosclerosis. Abdominal visceral fat mass is positively associated with carotid-femoral PWV in healthy subjects [24]. There is also a positive relationship between abdominal VFA and baPWV in diabetic patients [25] and diabetic HD patients [1]. However, an inverse association was found between abdominal visceral fatness and baPWV in female HD patients [2]. Thus, it is unclear whether the accumulation of visceral fat is a risk factor for arteriosclerosis in HD patients. In addition, overweight, hypertension, and hyperlipidemia, which are cardiovascular risk factors in the general population, have been reported to correlate with lower mortality in HD patients, a phenomenon that is known as “reverse epidemiology,” the “risk factor paradox,” and the “obesity paradox.”

In this study, increasing obesity in patients with high VFA and decreasing adiposity in patients with low VFA were associated with the development of arteriosclerosis. There are likely to be two different etiologies that may explain this phenomenon. One is the metabolic syndrome, and the other is the MIA syndrome. Obesity and insulin resistance are believed to be at the core of most cases of metabolic syndrome, and patients may have a greater risk of diseases such as cardiovascular disease [26]. The MIA syndrome features a strong relationship between malnutrition, inflammation, and atherosclerosis in dialysis patients, associated with a high mortality rate [3]. In this study, increasing obesity in patients with high VFA was associated with the development of arteriosclerosis, likely due to metabolic syndrome, whereas decreasing adiposity in patients with low VFA was also associated with the development of arteriosclerosis, but this was likely due to the MIA syndrome. Thus, the present study suggests that patients with a large visceral fat area are advised to adopt strategies aimed at reducing body weight, while patients with a small visceral fat area should avoid undernutrition. Patients failing to follow these instructions may have a high risk of arteriosclerosis and should be followed carefully.

Limitations

Several limitations of the present study should be mentioned. First, a small number of patients were studied. Second, because many patients were lost to follow-up because they died or because of missing data, no more than 40% of enrolled patients were analyzed for the influences of changes in obesity-related factors on arteriosclerosis. Third, we did not compare changes in baPWV with general population. Fourth, we do not have data investigating if reducing VFA and SFA can improve the prognosis of patients with large VFA. Moreover, associations between obesity-related parameters, such as adipocytokines and arteriosclerosis, were not evaluated in this study. Larger-scale longitudinal studies investigating the mechanism by which a change in VFA affects the progression of arteriosclerosis in HD patients are required.

Conclusions

In this study, we found that over the 4-year period, baPWV increased, and the change in baPWV was higher in patients who initially had a high VFA and gained additional adipose tissue mass, and was also higher in patients who initially had a low VFA and lost adipose tissue mass. Physicians should advise patients with a large visceral fat area not to gain weight, and conversely, those with small visceral fat area not to risk undernutrition.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ABI:

Ankle-brachial index

AC:

Abdominal circumference

baPWV:

Brachial-ankle pulse wave velocity

BMI:

Body mass index

BNP:

Brain natriuretic peptide

BPs:

Blood pressures

CRP:

C-reactive protein

CT:

Computed tomography

CTR:

Cardiothoracic ratio

CV:

Cardiovascular

HD:

Hemodialysis

hANP:

Human atrial natriuretic peptide

HDL-c:

High-density lipoprotein cholesterol

LDL-c:

Low-density lipoprotein cholesterol

MIA:

Malnutrition-inflammation-atherosclerosis

SFA:

Subcutaneous fat areas

VFA:

Visceral fat area

References

  1. Yurugi T, Morimoto S, Okamoto T, Amari Y, Kasuno Y, Fukui M, et al. Accumulation of visceral fat in maintenance hemodialysis patients. Clin Exp Nephrol. 2012;16(1):156–63.

    Article  CAS  Google Scholar 

  2. Kato A, Takita T, Kumagai H. Relationship between arterial stiffening and skeletal muscle atrophy in hemodialysis patients: a gender comparative study. J Cachexia Sarcopenia Muscle. 2014;5(3):247–9.

    Article  Google Scholar 

  3. Jeznach-Steinhagen A, Slotwinski R, Szczygiel B. Malnutrition, inflammation, atherosclerosis in hemodialysis patients. Rocz Panstw Zakl Hig. 2007;58(1):83–8.

    PubMed  Google Scholar 

  4. Ghiadoni L, Taddei S, Virdis A, Sudano I, Di Legge V, Meola M, et al. Endothelial function and common carotid artery wall thickening in patients with essential hypertension. Hypertension. 1998;32(1):25–32.

    Article  CAS  Google Scholar 

  5. Turin TC, Kita Y, Rumana N, Takashima N, Kadota A, Matsui K, et al. Brachial-ankle pulse wave velocity predicts all-cause mortality in the general population: findings from the Takashima study. Japan Hypertens Res. 2010;33(9):922–5.

    Article  Google Scholar 

  6. Morimoto S, Yurugi T, Aota Y, Sakuma T, Jo F, Nishikawa M, et al. Prognostic significance of ankle-brachial index, brachial-ankle pulse wave velocity, flow-mediated dilation, and nitroglycerin-mediated dilation in end-stage renal disease. Am J Nephrol. 2009;30(1):55–63.

    Article  Google Scholar 

  7. American DA. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2006;29(Suppl 1):S43–8.

    Google Scholar 

  8. Yamashina A, Tomiyama H, Takeda K, Tsuda H, Arai T, Hirose K, et al. Validity, reproducibility, and clinical significance of noninvasive brachial-ankle pulse wave velocity measurement. Hypertens Res. 2002;25(3):359–64.

    Article  Google Scholar 

  9. Yokoyama H, Shoji T, Kimoto E, Shinohara K, Tanaka S, Koyama H, et al. Pulse wave velocity in lower-limb arteries among diabetic patients with peripheral arterial disease. J Atheroscler Thromb. 2003;10(4):253–8.

    Article  Google Scholar 

  10. Kanda Y. Investigation of the freely available easy-to-use software ‘EZR’ for medical statistics. Bone Marrow Transplant. 2012;48:452.

    Article  Google Scholar 

  11. Examination Committee of Criteria for 'Obesity Disease' in J, Japan Society for the Study of O. New criteria for ‘obesity disease’ in Japan. Circ J. 2002;66(11):987–92.

  12. Kalantar-Zadeh K, Block G, Humphreys MH, Kopple JD. Reverse epidemiology of cardiovascular risk factors in maintenance dialysis patients. Kidney Int. 2003;63(3):793–808.

    Article  Google Scholar 

  13. Lowrie EG, Lew NL. Death risk in hemodialysis patients: the predictive value of commonly measured variables and an evaluation of death rate differences between facilities. Am J Kidney Dis. 1990;15(5):458–82.

    Article  CAS  Google Scholar 

  14. Cataliotti A, Malatino LS, Jougasaki M, Zoccali C, Castellino P, Giacone G, et al. Circulating natriuretic peptide concentrations in patients with end-stage renal disease: role of brain natriuretic peptide as a biomarker for ventricular remodeling. Mayo Clin Proc. 2001;76(11):1111–9.

    Article  CAS  Google Scholar 

  15. Zoccali C, Mallamaci F, Benedetto FA, Tripepi G, Parlongo S, Cataliotti A, et al. Cardiac natriuretic peptides are related to left ventricular mass and function and predict mortality in dialysis patients. J Am Soc Nephrol. 2001;12(7):1508–15.

    CAS  PubMed  Google Scholar 

  16. Ogata H, Kumasawa J, Fukuma S, Mizobuchi M, Kinugasa E, Fukagawa M, et al. The cardiothoracic ratio and all-cause and cardiovascular disease mortality in patients undergoing maintenance hemodialysis: results of the MBD-5D study. Clin Exp Nephrol. 2017;21(5):797–806.

    Article  CAS  Google Scholar 

  17. Nitta K, Masakane I, Hanafusa N, Taniguhi M, Hasegawa T, et al. Annual dialysis data report 2017, JSDT Renal Data Registry. Japanese Soc Dial Ther. 2018;51:699–766.

    Article  Google Scholar 

  18. Levey AS, Beto JA, Coronado BE, Eknoyan G, Foley RN, Kasiske BL, et al. Controlling the epidemic of cardiovascular disease in chronic renal disease: what do we know? What do we need to learn? Where do we go from here? National Kidney Foundation Task Force on Cardiovascular Disease. Am J Kidney Dis. 1998;32(5):853–906.

    Article  CAS  Google Scholar 

  19. Tong J, Liu M, Li H, Luo Z, Zhong X, Huang J, et al. Mortality and associated risk factors in dialysis patients with cardiovascular disease. Kidney Blood Press Res. 2016;41(4):479–87.

    Article  CAS  Google Scholar 

  20. Mark PB. Strategies to manage cardiovascular risk in chronic kidney disease. Nephrol Dial Transplant. 2018;33(1):23–5.

    Article  CAS  Google Scholar 

  21. Fischer EC, Zocalo Y, Galli C, Wray S, Bia D. Arterial stiffness and renal replacement therapy: a controversial topic. Int J Nephrol. 2015. https://doi.org/10.1155/2015/729609.

  22. Avramovski P, Janakievska P, Sotiroski K, Sikole A. Accelerated progression of arterial stiffness in dialysis patients compared with the general population. Korean J Intern Med. 2013;28(4):464–74.

    Article  Google Scholar 

  23. Tomiyama H, Matsumoto C, Shiina K, Yamashina A. Brachial-ankle PWV: current status and future directions as a useful marker in the management of cardiovascular disease and/or cardiovascular risk factors. J Atheroscler Thromb. 2016;23(2):128–46.

    Article  CAS  Google Scholar 

  24. Lee Y, Shin H, Vassy JL, Kim JT, Cho SI, Kang SM, et al. Comparison of regional body composition and its relation with cardiometabolic risk between BMI-matched young and old subjects. Atherosclerosis. 2012;224(1):258–65.

    Article  CAS  Google Scholar 

  25. Anan F, Masaki T, Umeno Y, Iwao T, Yonemochi H, Eshima N, et al. Correlations of visceral fat accumulation and atherosclerosis in Japanese patients with type 2 diabetes mellitus. Metabolism. 2008;57(2):280–4.

    Article  CAS  Google Scholar 

  26. Samson SL, Garber AJ. Metabolic syndrome. Endocrinol Metab Clin N Am. 2014;43(1):1–23.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Mark Cleasby, PhD, from Edanz Group (www.edanzediting.com/ac) for editing a draft of this manuscript.

Funding

This work was supported in part by grants from the Ministry of Education, Culture, Sports, Science and Technology of Japan (19K08497 to SM).

Author information

Authors and Affiliations

Authors

Contributions

All author contributed to this manuscript. Conception and design, analysis and interpretation of data: TI, SM, and AI. Drafting and revising the manuscript: TI and SM. Providing intellectual content of critical importance to the work described: TI and SM. Collection and assembly of data: TI, HO, YA, TY, and FN. Final approval of the version to be published: AI. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Satoshi Morimoto.

Ethics declarations

Ethics approval and consent to participate

This prospective study was conducted in accordance with the Declaration of Helsinki and good clinical practice guidelines and was approved by the ethics committee of Moriguchi Keijinkai Hospital (approval no.: 3003). All participants were enrolled after providing written informed consent.

Consent for publication

All participants agreed for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Iida, T., Morimoto, S., Okuda, H. et al. Influences of changes in obesity-related parameters on progression of arteriosclerosis in hemodialysis patients: a prospective cohort study. Ren Replace Ther 6, 54 (2020). https://doi.org/10.1186/s41100-020-00304-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s41100-020-00304-w

Keywords