- 1Department of Cardiology, Toho University Sakura Medical Center, Chiba, Japan
- 2Mihama Hospital, Chiba, Japan
During hemodialysis (HD), blood pressure (BP) changes are frequently observed. However, the mechanism of BP changes during HD has not been fully clarified. The cardio–ankle vascular index (CAVI) reflects the arterial stiffness of the arterial tree from the origin of the aorta to the ankle independent from BP during measurement. Additionally, CAVI reflects functional stiffness in addition to structural stiffness. We aimed to clarify the role of CAVI in regulating the BP system during HD. We included 10 patients undergoing 4-hour HD (total 57 HD sessions). Changes in the CAVI and various hemodynamic parameters were evaluated during each session. During HD, BP decreased and CAVI significantly increased (CAVI, median [interquartile range]; 9.1 [8.4–9.8] [0 min] to 9.6 [9.2–10.2] [240 min], p < 0.05). Changes in CAVI from 0 min to 240 min were significantly correlated with water removal rate (WRR) (r = −0.42, p = 0.002). Changes in CAVI at each measurement point were negatively correlated with ΔBP (Δsystolic BPeach MP, r = −0.23, p < 0.0001; Δdiastolic BPeach MP, r = −0.12, p = 0.029). Whereas one patient exhibited a simultaneous decrease in BP and CAVI during the initial 60 min of HD. Arterial stiffness monitored with CAVI generally increased during HD. CAVI elevation is associated with decreased WWR and BP. An increase in CAVI during HD may reflect the contraction of smooth muscle cells and play an important role in BP maintenance. Hence, measuring CAVI during HD may distinguish the cause of BP changes.
1 Introduction
Blood pressure (BP) changes can occur during hemodialysis (HD) and hypotension attack and is associated with severe clinical symptoms such as fainting, discomfort, and sweating. Intra-dialytic hypotension occurs in 10%–12% of HD cases and is related to poor prognosis (Gul et al., 2016; Sars et al., 2020). The changes in BP during HD is thought to be mainly caused by the removal of the volume of the water component. Other factors involved in BP changes are the autonomic nervous system and cardiac function (Sars et al., 2020). In addition, the contraction of arterial smooth muscle cells may be involved in controlling BP, but understanding the state of arterial smooth muscle contraction or relaxation is difficult.
Arterial stiffness consists of structural and functional stiffness. The former mainly reflects arteriosclerosis and vascular aging. The latter reflects arterial smooth muscle contraction (Lacolley et al., 2017). As an index of arterial stiffness, pulse wave velocity (PWV) has been used to mainly evaluate the degree of arteriosclerosis or vascular aging, and there were several reports that PWV played a predictive role in cardiovascular events (Blacher et al., 1999; Vlachopoulos et al., 2010). As for changes in PWV during HD, Lafta et al. (2022) demonstrated that there were no significant changes in PWV during HD. Öğünç et al. (2015) reported that BP and PWV decreased during HD. However, PWV is dependent on BP during measurement (Cecelja and Chowienczyk, 2009; Schillaci et al., 2015). Therefore, PWV is not an appropriate index for the evaluation of arterial stiffness if BP changes during HD.
Recently, the cardio–ankle vascular index (CAVI) was developed as a BP-independent arterial stiffness index by applying the stiffness parameter ß theory (Shirai et al., 2006). The CAVI reflects arterial stiffness from the origin of the aorta to the ankle, and arterial stiffness indicated by CAVI is composed of structural stiffness, which reflects the degree of arteriosclerosis (Shirai et al., 2011a), and functional stiffness (Shirai et al., 2011b). Functional stiffness reflects contraction or relaxation of arterial smooth muscle cells. Thus, CAVI may be a useful index for evaluating the role of arterial smooth muscle contraction in various pathological conditions. Nagasawa et al. (2022) reported that CAVI increased in hypovolemic states caused by acute blood removal in rabbits. This result strongly suggests that CAVI could reflect systemic vascular contraction in hypovolemic states. This hypovolemic state resembles the hemodynamic condition during HD. Therefore, measuring CAVI during HD may provide insight into understanding vascular conditions and may lead to clarify the pathophysiological mechanisms of the BP change during HD. CAVI could be a novel modality for assessing vascular physiology during HD.
In this study, we measured CAVI during HD to clarify the role of CAVI in maintaining systemic circulation.
2 Materials and methods
We enrolled 10 patients who underwent a 4-hour HD therapy at Mihama Hospital between March 2015 and May 2015 for a total of 57 HD sessions (6 HD sessions per patient, 1 patient had 3 HD sessions). Patients with persistent atrial fibrillation, peripheral arterial disease with ankle–brachial index <0.9, and severe aortic valvular heart disease were excluded from this study.
We analyzed the systolic BP (SBP), diastolic BP (DBP), heart rate (HR), heart–ankle PWV (haPWV), CAVI, and stroke volume (SV) during HD. The BP, HR, haPWV, and CAVI were measured using VaSela VS-1500 (Fukuda Denshi Co. Ltd., Tokyo, Japan) after 10 min of rest in the supine position. Simultaneously, SV was non-invasively measured using an impedance cardiac output meter (AESCULON; Osypka Medical Inc., Berlin, Germany). Systemic vascular resistance (SVR) was calculated as the ratio of mean BP to cardiac output (CO).
To evaluate the intra-dialytic changes in arterial stiffness and hemodynamics, each parameter was measured at the following points during each HD session: before HD (0 min), 30 min, 90 min, 150, 210, and 240 min, and after blood return. To examine the relationship between changes in CAVI or various hemodynamic parameters and water removal volume (WRV), the correlation between changes in CAVI or other parameters from 0 min to 240 min (ΔCAVI0–240 min, ΔSBP0–240 min, ΔDBP0–240 min, ΔHR0–240 min, ΔSV0–240 min, ΔCO0–240 min, and ΔSVR0–240 min) and water removal rate (WRR), which was defined as the percentage of WRV per body weight before HD, were analyzed. To examine the relationship between CAVI and other parameters at each measurement point (MP), the correlation between the changes in CAVI at each MP (ΔCAVIeach MP) and the changes in other parameters at the same time (ΔSBPeach MP, ΔDBPeach MP, ΔHReach MP, ΔSVeach MP, ΔCOeach MP, and ΔSVReach MP) were also analyzed.
Furthermore, CAVI was measured every 30 min in one patient whose BP decreased in the initial 60 min during HD.
2.1 Statistical analyses
Data normality was assessed using the Shapiro–Wilk test. For patient characteristics, continuous variables are expressed as mean ± standard deviation or median (interquartile range [IQR]), as appropriate. Categorical variables are expressed as percentages. Hemodynamic parameters, including CAVI, are expressed as median (IQR). Comparisons of hemodynamic parameters during HD were performed using the Friedman test and the post hoc Steel–Dwass test. We analyzed the association between changes in CAVI and other hemodynamic parameters using Spearman’s correlation coefficient. The statistical significance was set at a p < 0.05. Statistical analysis was performed using EZR on R commander statistical software version 1.55 (Saitama Medical Center, Jichi Medical University, Saitama, Japan) (Kanda, 2013) and JMP computer software version 14.2 (SAS Institute, Cary, NC, United States).
2.2 Ethics approval
Informed consent was obtained from all patients. This study was approved by the ethics committee of Toho University Sakura Medical Center (S18019) and Mihama Hospital (17-010) and was performed in accordance with the latest version of the tenets of the Declaration of Helsinki.
3 Results
3.1 Patient characteristics
A total of 10 patients undergoing HD were enrolled in this study. All patients received 4-hour HD for 3 times per week; 9 patients who had 6 HD sessions and 1 patient who had 3 HD sessions were enrolled in the study, thus providing a total of 57 HD sessions for analysis.
Table 1 summarizes patient characteristics. The study sample comprised 50% men, and the mean age was 69.3 ± 9.3 years. The median dialysis duration was 4.7 years. The percentage of patients with diabetes mellitus was 20%, and the percentage of those with a history of coronary artery disease was also 20%. Echocardiographic findings revealed a mean left ventricular ejection fraction of 64.0% ± 6.7%. The mean E/e' was 17.5 ± 4.3, thus indicating left ventricular diastolic dysfunction. The rate of use of oral antihypertensive agents and hypoglycemic agents was low (Ca channel blocker, 10%; angiotensin Ⅱ receptor blocker, 20%; ß-blocker, 30%; and hypoglycemic agent, 10%). Mean total WRV per single HD session was 2992.6 ± 564.2 mL, and mean WRR was 5.5% ± 1.3%.
3.2 Changes in blood pressure, cardio–ankle vascular index, and other hemodynamic parameters during hemodialysis
Figure 1A depicts the changes in BP, CAVI, and haPWV during HD. The SBP significantly decreased from 30 min to 240 min and after blood return. The DBP significantly decreased from 90 min to 240 min. After blood return, both SBP and DBP increased compared with their values at 240 min. In contrast to the BP, CAVI significantly increased from 150 min to 240 min. For haPWV, there were no significant changes during HD or after blood return.
FIGURE 1. Changes in BP, CAVI, and other hemodynamic parameters during hemodialysis (A) Changes in BP, CAVI, and haPWV during hemodialysis. (B) Changes in CO, HR, SV, and SVR during hemodialysis. Data are shown as median (interquartile range). White bars, median; boxes, interquartile range; whisker, ranges excluding statistical outliers; plots, statistical outliers. SBP, systolic blood pressure; DBP, diastolic blood pressure; CAVI, cardio–ankle vascular index; haPWV, heart–ankle pulse wave velocity; CO, cardiac output; HR, heart rate; SV, stroke volume; SVR, systemic vascular resistance. *, p < 0.05 vs. 0 min; **, p < 0.05 vs. 240 min.
Figure 1B depicts the changes in CO, HR, SV, and SVR during HD. CO significantly decreased from 30 min to 240 min. HR did not change significantly during HD but increased slightly from 210 min to 240 min. SV significantly decreased from 90 min to 240 min and after blood return. After blood return, SV tended to increase compared with that at 240 min but was not statistically significant. SVR significantly increased at 90 min, but there was no significant change at other MPs.
3.3 Correlation between changes in cardio–ankle vascular index or blood pressure and water removal rate during hemodialysis
Figure 2 shows the correlation between ΔCAVI0–240 min or ΔBP0–240 min and WRR. ΔCAVI0–240 min was significantly correlated with WRR (r = 0.42, p = 0.002; Figure 2A). ΔSBP0–240 min and ΔDBP0–240 min tended to correlate with WRR, but the correlation was not statistically significant (ΔSBP0–240 min: r = −0.16, p = 0.25; ΔDBP0–240 min: r = −0.15, p = 0.31 vs. WWR; Figures 2B,C). ΔHR0–240 min, ΔSV0–240 min, ΔCO0–240 min, and ΔSVR0–240 min were not significantly correlated with WRR (ΔHR0–240 min: r = 0.22, p = 0.12; ΔSV0–240 min: r = 0.10, p = 0.48; ΔCO0–240 min: r = 0.17, p = 0.22; ΔSVR0–240 min: r = −0.17, p = 0.21 vs. WWR).
FIGURE 2. Correlation between ΔCAVI0–240 min or ΔBP0–240 min, and WRR during hemodialysis (A) Correlation between ΔCAVI 0–240 min and WRR during hemodialysis. (B) Correlation between ΔSBP0–240 min and WRR during hemodialysis. (C) Correlation between ΔDBP0–240 min and WRR during hemodialysis. r, Spearman’s correlation coefficient; CAVI, cardio–ankle vascular index; BP, blood pressure; WRR, water removal rate; SBP, systolic blood pressure; DBP, diastolic blood pressure.
3.4 Correlation between changes in cardio–ankle vascular index and blood pressure at each measurement point during hemodialysis
Figure 3 shows the correlation between ΔCAVI and ΔSBP or ΔDBP at each MP during HD. ΔCAVIeach MP significantly correlated with ΔSBPeach MP and ΔDBPeach MP (ΔSBPeach MP: r = −0.23, p < 0.0001; ΔDBPeach MP: r = −0.12, p = 0.029; Figures 3A,B). In other parameters, ΔSVeach MP and ΔCOeach MP did not correlate with ΔCAVIeach MP (ΔSVeach MP: r = 0.08, p = 0.15; ΔCOeach MP: r = −0.008, p = 0.88). ΔCAVIeach MP was positively correlated with ΔHReach MP and was negatively correlated with ΔSVReach MP (ΔHReach MP: r = 0.21, p = 0.001; ΔSVReach MP: r = −0.14, p = 0.01).
FIGURE 3. Correlation between ΔCAVI and ΔBP at each MP during hemodialysis (A) Correlation between ΔCAVIeach MP and ΔSBPeach MP during hemodialysis. (B) Correlation between ΔCAVIeach MP and ΔDBPeach MP during hemodialysis. r, Spearman’s correlation coefficient; CAVI, cardio–ankle vascular index; MP, measurement point; BP, blood pressure; SBP, systolic blood pressure; DBP, diastolic blood pressure.
3.5 A case with decreased blood pressure and cardio–ankle vascular index
An 85-year-old man exhibited a decrease in BP at 60 min during HD, as shown in Figure 4. At that moment, CAVI also decreased accompanying with the decrease in BP.
FIGURE 4. A case in which CAVI decreases with a reduction in BP at the early phase of hemodialysis. Case. An 85-year-old man with diabetic nephropathy. CAVI, cardio–ankle vascular index; SBP, systolic blood pressure; DBP, diastolic blood pressure.
4 Discussion
In present study, we observed changes in hemodynamic parameters and arterial stiffness monitored with CAVI during HD, and the relationship between hemodynamic parameters and CAVI was evaluated. The SBP and DBP started to decrease at 30 min; this decrease continued until 240 min and then increased slightly after blood return, as shown in Figure 1A. A decrease in CO is also shown in Figure 1B. The observed decrease in these parameters was due to the decrease in circulating blood volume by water removal with HD. The WRR was negatively related to the changes in SBP and DBP but not significantly (SBP: r = −0.16, p = 0.25; DBP: r = 0.15, p = 0.31; Figures 2B,C).
Conversely, CAVI gradually increased from 30 min, and the increase was significant at 150 min to 240 min. The contradictory changes between BP and CAVI were also confirmed by significant correlations between ΔSBP and ΔDBP and ΔCAVI at each MP, as shown in Figure 3.
Reportedly, CAVI reflects functional stiffness in addition to structural stiffness of the arterial tree (Shirai et al., 2011b; Takahashi et al., 2012; Chiba et al., 2015; Shimizu et al., 2016; Sato et al., 2020). The adrenaline receptor α blocker doxazosin, which dilates arterial smooth muscles cells, is reported to decrease BP and decrease CAVI (Shirai et al., 2011b). Nitroglycerine, which also dilates smooth muscle cells, decreases CAVI (Chiba et al., 2015; Shimizu et al., 2016). Nagasawa et al. (2022) reported that blood removal provokes a decrease in BP and an increase in CAVI in rabbits. Moreover, in this experimental study, actual vasoconstriction was observed using intra-vascular ultrasound. This fact strongly suggests that the CAVI increase reflects vasoconstriction during blood removal. These results indicate that CAVI reflects the state of contraction and/or dilatation of smooth muscle cells in the arterial tree and that CAVI may be involved in the control system of BP as a part of resistance. From this viewpoint, an increase in CAVI during HD reflects smooth muscle contraction in the arterial tree from the origin of the aorta to the ankle and increases BP as a part of systemic resistance. The mechanism by which CAVI increases during blood removal may be via the stimulation of the sympathetic autonomic nervous system by hypotension. However, more studies are required to understand this further.
A weak negative relationship between WRR and BP (Figures 2B,C) was observed, but it was not statistically significant. The lack of significance in the relationship between WRR and BP may be owing to the compensatory reaction of CAVI; that is, an increase in CAVI by blood removal may have increased the peripheral arterial resistance and prevented a decrease in BP because the relationship between WRR and ΔCAVI was strongly positive (Figure 2A). From these results and speculations, we suggest that CAVI may play an important role in regulating systemic BP via peripheral artery resistance (Figure 5).
FIGURE 5. Role of CAVI in hemodynamic maintenance during hemodialysis. CAVI, cardio–ankle vascular index.
Furthermore, haPWV did not change during HD, as shown in Figure 1A. PWV essentially depends on BP upon measurement. Enhancing arterial stiffness brought about by water loss during HD should enhance PWV, but a decrease in BP reduced PWV, which depended on BP at measurement. Therefore, haPWV seemed to be unchanged. Furthermore, PWV is not a suitable index for reflecting proper arterial stiffness.
The above results suggested that an increase in CAVI during HD reflects the contraction of smooth muscle cells. This increase in CAVI might contribute to BP maintenance, which decreased during HD by decreasing the circulating blood brought about by water removal (Figure 5).
An 85-year-old man exhibited a decrease in BP at 60 min of HD, as shown in Figure 4. At this time, CAVI also decreased with the decrease in BP. This phenomenon demonstrates that smooth muscle cells dilate when BP decreases (Figure 5). The precise mechanism for smooth muscle dilation is not yet clear. However, those observations are important, because CAVI can differentiate the causes of a decrease in BP during HD. One general cause of a decrease in BP accompanying the elevation in CAVI is due to blood volume loss. The other cause of a decrease in BP accompanying a decrease in CAVI is due to smooth muscle cell dilatation or relaxation. Thus, the cause of a decrease in BP during HD might be identified by measuring CAVI.
In summary, arterial stiffness monitored with CAVI generally increased during HD. In addition, this CAVI increase was associated with an increase in WWR and a decrease in BP. An increase in CAVI during HD may reflect the contraction of smooth muscle cells and plays an important role in maintaining BP. On the basis of this hypothesis, measuring CAVI during HD may help to clarify the cause of BP change.
5 Study limitations
First, this was an observational study with a small sample size. We included only 10 HD patients that was not enough size for statistical analysis. Therefore, we analyzed total 57 HD sessions (6 sessions per patient except 1 patient) that were measured CAVI at same time points. Second, although the mechanism of increased CAVI was presumed to be due to the activation of the sympathetic nervous system, the precise mechanism has not been verified. Third, the number of cases in which CAVI decreased during HD was small, and we were not able to show mechanisms for the decrease in CAVI. Therefore, further studies with a large number of patients and on precise mechanism are needed to clarify.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the ethics committee of the Toho University Sakura Medical Center (S18019) and Mihama Hospital (17-010). The patients/participants provided their written informed consent to participate in this study.
Author contributions
SS, KaS, and KoS. conceived, desined the study. KaS and KoS supervised the study. SS, JU, TS, YS and NI collected data. SS, KaS, and MT analyzed data and interpreted results. SS and KaS. drafted the manuscript. MM and ON clinically supervied the study. All authers edited and revised the manuscript. All authors approved the published version of the manuscript.
Acknowledgments
We would like to thank the cardiology staff and clinical physiology staff of Toho University Sakura Medical Center and Mihama Hospital.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Blacher J., Guerin A. P., Pannier B., Marchais S. J., Safar M. E., London G. M. (1999). Impact of aortic stiffness on survival in end-stage renal disease. Circulation 99, 2434–2439. doi:10.1161/01.cir.99.18.2434
Cecelja M., Chowienczyk P. (2009). Dissociation of aortic pulse wave velocity with risk factors for cardiovascular disease other than hypertension: A systematic review. Hypertension 54, 1328–1336. doi:10.1161/HYPERTENSIONAHA.109.137653
Chiba T., Yamanaka M., Takagi S., Shimizu K., Takahashi M., Shirai K., et al. (2015). Cardio-ankle vascular index (CAVI) differentiates pharmacological properties of vasodilators nicardipine and nitroglycerin in anesthetized rabbits. J. Pharmacol. Sci. 128, 185–192. doi:10.1016/j.jphs.2015.07.002
Gul A., Miskulin D., Harford A., Zager P. (2016). Intradialytic hypotension. Curr. Opin. Nephrol. Hypertens. 25, 545–550. doi:10.1097/MNH.0000000000000271
Kanda Y. (2013). Investigation of the freely available easy-to-use software ‘EZR’ for medical statistics. Bone Marrow Transpl. 48, 452–458. doi:10.1038/bmt.2012.244
Lacolley P., Regnault V., Segers P., Laurent S. (2017). Vascular smooth muscle cells and arterial stiffening: Relevance in development, aging, and disease. Physiol. Rev. 97, 1555–1617. doi:10.1152/physrev.00003.2017
Lafta A., Ukrainetz J., Davison S., Thompson S., Bello A., Braam B. (2022). Inter- and intradialytic fluid volume changes and vascular stiffness parameters in patients on hemodialysis. PLOS ONE 17, e0262519. doi:10.1371/journal.pone.0262519
Nagasawa Y., Shimoda A., Shiratori H., Morishita T., Sakuma K., Chiba T., et al. (2022). Analysis of effects of acute hypovolemia on arterial stiffness in rabbits monitored with cardio-ankle vascular index. J. Pharmacol. Sci. 148, 331–336. doi:10.1016/j.jphs.2022.01.008
Öğünç H., Akdam H., Alp A., Gencer F., Akar H., Yeniçerioğlu Y. (2015). The effects of single hemodialysis session on arterial stiffness in hemodialysis patients. Hemodial. Int. 19, 463–471. doi:10.1111/hdi.12277
Sars B., van der Sande F. M., Kooman J. P. (2020). Intradialytic hypotension: Mechanisms and outcome. Blood Purif. 49, 158–167. doi:10.1159/000503776
Sato S., Takahashi M., Mikamo H., Kawazoe M., Iizuka T., Shimizu K., et al. (2020). Effect of nicorandil administration on cardiac burden and cardio-ankle vascular index after coronary intervention. Heart Vessels 35, 1664–1671. doi:10.1007/s00380-020-01650-9
Schillaci G., Battista F., Settimi L., Anastasio F., Pucci G. (2015). Cardio-ankle vascular index and subclinical heart disease. Hypertens. Res. 38, 68–73. doi:10.1038/hr.2014.138
Shimizu K., Yamamoto T., Takahashi M., Sato S., Noike H., Shirai K. (2016). Effect of nitroglycerin administration on cardio-ankle vascular index. Vasc. Health Risk Manag. 12, 313–319. doi:10.2147/VHRM.S106542
Shirai K., Hiruta N., Song M., Kurosu T., Suzuki J., Tomaru T., et al. (2011a). Cardio-ankle vascular index (CAVI) as a novel indicator of arterial stiffness: Theory, evidence and perspectives. J. Atheroscler. Thromb. 18, 924–938. doi:10.5551/jat.7716
Shirai K., Song M., Suzuki J., Kurosu T., Oyama T., Nagayama D., et al. (2011b). Contradictory effects of β1-and α1-aderenergic receptor blockers on cardio-ankle vascular stiffness index (CAVI)--CAVI independent of blood pressure. J. Atheroscler. Thromb. 18, 49–55. doi:10.5551/jat.3582
Shirai K., Utino J., Otsuka K., Takata M. (2006). A novel blood pressure-independent arterial wall stiffness parameter; cardio-ankle vascular index (CAVI). J. Atheroscler. Thromb. 13, 101–107. doi:10.5551/jat.13.101
Takahashi M., Shiba T., Hirano K., Hitsumoto T., Shirai K. (2012). Acute decrease of cardio-ankle vascular index with the administration of Beraprost sodium. J. Atheroscler. Thromb. 19, 479–484. doi:10.5551/jat.9266
Keywords: CAVI, arterial stiffness, hemodialysis, intra-dialytic hypotension, vascular function
Citation: Sato S, Shimizu K, Takahashi M, Masai M, Nagakawa O, Uchino J, Suzuki T, Sato Y, Iwai N and Shirai K (2023) Changes in blood pressure and arterial stiffness monitored using the cardio–ankle vascular index during hemodialysis. Front. Physiol. 14:1133037. doi: 10.3389/fphys.2023.1133037
Received: 28 December 2022; Accepted: 07 February 2023;
Published: 20 February 2023.
Edited by:
Antonio Colantuoni, University of Naples Federico II, ItalyReviewed by:
Romeo Martini, University Hospital of Padua, ItalyDominga Lapi, University of Naples Federico II, Italy
Copyright © 2023 Sato, Shimizu, Takahashi, Masai, Nagakawa, Uchino, Suzuki, Sato, Iwai and Shirai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shuji Sato, c2h1dWppLnNhdG91QG1lZC50b2hvLXUuYWMuanA=