Skip to main content
  • Research article
  • Open access
  • Published:

Estimating the extent of subclinical arteriosclerosis of persons with prediabetes and diabetes mellitus among Japanese urban workers and their families: a cross-sectional study

Abstract

Background

Diabetes mellitus (hereafter called diabetes) is considered to accelerate arteriosclerosis leading to coronary heart disease and stroke. Thus, it is important to quantitatively estimate the extent of subclinical arteriosclerosis. A new method called cardio-ankle vascular index (CAVI) is developed to reflect arterial stiffness independently from blood pressure at the time of measurement. Then, we examined if CAVI scores could discriminate the extent of arteriosclerosis between persons with prediabetes (or borderline diabetes) and with diabetes among Japanese urban workers and their families.

Methods

Subjects were 9881 men and 12033 women of company employees and their families who participated in cardiovascular disease screening in Japan. Persons having diabetes and prediabetes were defined based on the criteria set by American Diabetes Association. CAVI scores were measured by VaSera VS-1000. We applied the established age-sex specific cutoff points of CAVI scores above which were determined to be abnormally high or advanced level of arteriosclerosis. To examine the association of prediabetes and diabetes with CAVI scores, CAVI scores of screening participants were converted to a binary variable: 1 for less than cutoff points and 2 for equal or greater than cutoff points or abnormally high CAVI scores. Logistic regression method was used to examine the association of prediabetes and diabetes with CAVI scores after adjusting for major cardiovascular disease (CVD) risk factors.

Results

Prevalence of abnormally high CAVI scores was significantly higher after 40 years of age among persons with diabetes than either among persons with prediabetes or among normal persons in both genders. Significantly elevated odds ratios (ORs) of abnormally high CAVI scores appeared among persons with prediabetes: 1.29 (95 % confidence interval (CI), 1.11-1.48) for men and 1.14 (CI, 1.01-1.28) for women, and among persons with diabetes: 2.41 (CI, 1.97-2.95) for men and 2.52 (CI, 1.94-3.28) for women.

Conclusions

The extent of subclinical arteriosclerosis (including arterial stiffness and atherosclerosis) was moderately enhanced among persons with prediabetes and was further advanced among persons with diabetes. Thus, it is important to introduce earlier interventions for changing lifestyle and diet of persons with prediabetes in order to prevent them from developing diabetes and further advancing arteriosclerosis.

Peer Review reports

Background

One of the most serious complications among person with diabetes is the deterioration of arterial system or the acceleration of arteriosclerosis which could not be easily measured in the past. One method to quantitatively estimate the extent of arteriosclerosis is the use of the pulse wave velocity (PWV). The idea on the association of PWV with arteriosclerosis is traced back to an experiment using artificial blood vessels conducted by Moens in 1878 [1]. Then, Bramwell and colleagues showed that PWV depends on the modulus of arterial volume elasticity by experiments in 1922-23 [15]. Their experimental results have been a basis for the development of the measurement device PWV-200 (Fukuda-Denshi Co., Tokyo) which measures PWV propagating through the aorta (thorax, abdomen, and part of common iliac artery) from the aortic valve to the femoral pulsation point, as described by Hasegawa in 1970 [6]. Because PWV is highly correlated with diastolic blood pressure, Hasegawa developed a nomogram showing the association between diastolic blood pressure and PWV. He proposed an adjustment to any measured PWV values at 80 mmHg. As a result, such an adjustment was built into the PWV-200 machine. This is an important step allowing clinicians and researchers to compare PWV values between individuals and between populations. Namekata et al. conducted cardiovascular disease prevention screening in Seattle and found that PWV was positively and significantly associated with aging (≥60 years of age), hypertension, diabetes, the ratio of total cholesterol to high density lipoprotein cholesterol, ex-smokers and negatively and significantly with alcohol consumption among Japanese Americans [7]. In addition, they had similar findings among Japanese urban workers [8].

To overcome some technical difficulty for measuring PWV, the cardio-ankle vascular index (CAVI) was developed as a new indicator of arteriosclerosis in 2004 [9]. CAVI scores quantitatively reflect arteriosclerosis of the aorta, femoral and tibial arteries based on Bramwell-Hill’s equation [1] and stiffness parameter [10] which is allowed to be converted from PWV propagating from the aortic valve to ankle. Some researchers proposed to use CAVI scores as an indicator of atherosclerosis. Nakamura et al. found a strong association of CAVI with the presence of severity of coronary atherosclerosis based on their ordinal logistic regression analysis [11]. Kadota et al. suggested the use of CAVI as a screening tool for atherosclerosis based on their findings from the general population study of 1,014 adults showing strong significant associations of CAVI scores with carotid intima-media thickness and with homocysteine after adjustment for age and sex [12]. Thus, it is considered that CAVI scores reflect arterial stiffness, atherosclerosis and arteriosclerosis of which conditions are overlapping and inseparable. We use CAVI scores to represent the extent of subclinical arteriosclerosis in this paper but it is inclusive of arterial stiffness and atherosclerosis.

To practically use CAVI as a diagnostic tool for determining the extent of arteriosclerosis, our previous study established the baseline CAVI scores by age and gender among cardiovascular disease (CVD) risk-free persons [13]. In our present study we measured CAVI scores and identified persons with abnormally high CAVI scores in Japanese urban workers and their families based on the criteria developed by our studies [13, 14], and then examined the extent of subclinical arteriosclerosis among persons with prediabetes and among persons with diabetes mellitus after making adjustment for major cardiovascular disease (CVD) risk factors.

Methods

Subjects

Subjects for the study were recruited from January 2006 to May 2009 through the screening program at Japan Health Promotion Foundation which has been conducting cardiovascular disease and cancer screening throughout major cities in Japan. Subjects were company employees and their family members: 9,881 men and 12,033 women between 17 and 87 years of age. Following Japan’s Personal Information Protection Law, the use of anonymized screening data for research purpose was approved by the board of Japan Health Promotion Foundation. The study was approved by the board of Pacific Rim Disease Prevention Center.

Definition of prediabetes and diabetes mellitus

Following the recommendation from American Diabetes Association [15], persons having diabetes were defined as those having medical history of diabetes and/or taking medication of diabetes and/or whose fasting plasma glucose were equal to or higher than 126 mg/dl and/or whose hemoglobinA1c (HbA1c) were equal to or higher than 6.5 % in NGSP value (48 mmol/mol in IFCC). Persons with the state of prediabetes were defined as those whose fasting plasma glucose was from 100 mg/dl to 125 mg/dl and/or whose HbA1c was from 5.7 % to 6.4 % in NGSP value (39 mmol/mol to 47 mmol/mol in IFCC). Persons defined as “normal” were those without diabetes and prediabetes.

Clinical measurements

Blood was drawn from subjects after a 12 h fast. The following measurements were made: total cholesterol (TC) and triglycerides (TG) by enzymatic assay; high density lipoprotein cholesterol (HDL-C) by modified enzymatic method; glucose by hexokinase glucose-6-phosphate dehydrogenate assay; and glyco-hemoglobin A1c (HbA1c) by latex agglutination.

Following the guideline released by American Heart Association in 2007 [16], persons having hypertension were defined as those having medical history of hypertension and/or taking hypertension drugs and/or whose systolic blood pressure (SBP) was equal to or higher than 140 mmHg and/or whose diastolic blood pressure (DBP) was equal to or higher than 90 mmHg.

Cardio-ankle vascular index

CAVI, a stiffness and arteriosclerosis indicator of thorax, abdomen, common iliac, femoral and tibial arteries, was measured by VaSera VS-1000 manufactured by Fukuda-Denshi Company, LTD (Tokyo, Japan).

As illustrated by Shirai et al. [9], the scale conversion from PWV to CAVI is performed by the following formula:

$$ \mathrm{CAVI}=\mathrm{a}\left\{\left(2\uprho /\Delta \mathrm{P}\right)\ \mathrm{x}\ \ln \left(\mathrm{P}\mathrm{s}/\mathrm{P}\mathrm{d}\right){\mathrm{PWV}}^2\right\} + \mathrm{b} $$

where Ps and Pd are systolic and diastolic blood pressure values, respectively, PWV is the pulse wave velocity between heart and ankle, ∆P is Ps-Pd, ρ is blood density, and a and b are constants. This equation was derived from Bramwell-Hill’s equation [1] and stiffness parameter [10]. Scale conversion constants are determined so as to match CAVI with PWV by Hasegawa’s method [6]. These measurements and calculations are automatically made in VaSera VS-1000.

In the previous study we established the age-sex specific cutoff points of CAVI scores above which were determined to be abnormally high or advanced level of arteriosclerosis [14]. The cutoff points are (mean of CAVI + one standard deviation ) among CVD risk-free subjects: 7.39 for 20-29 years of age, 7.80 for 30-39 years of age, 8.29 for 40-49 years of age, 8.83 for 50-59 years of age, 9.54 for 60-69 years of age, and 10.35 for 70 years of age and over among men; and 7.23 for 20-29 years of age, 7.42 for 30-39 years of age, 7.95 for 40-49 years of age, 8.52 for 50-59 years of age, 8.98 for 60-69 years of age, and 9.46 for 70 years of age and over among women. To apply the logistic regression method for examining the association of prediabetes and diabetes with CAVI scores, CAVI scores of screening participants were converted to a binary variable: 1 for less than cutoff points and 2 for equal or greater than cutoff points or abnormally high CAVI scores.

Questionnaire

A short self-administered questionnaire was filled out by each subject during the screening. It contains questions on medical history and lifestyle factors such as smoking habit and alcohol consumption.

Statistical methods

All statistical analyses were performed gender-specifically. To examine characteristics of study participants by diabetes status, Student’s t-tests and chi-square tests were conducted for detecting significant differences in means and in prevalence, respectively, between persons with normal status and persons with prediabetes and diabetes. Cochran-Armitage test for linear trend was applied to evaluate the dose-dependent association between the degree of glycemic status and the prevalence of abnormally high CAVI scores. Crude, age-adjusted and multivariable-adjusted odds ratios (OR) and the 95 % confidence intervals (CI) of abnormally high CAVI scores according to diabetic status were calculated in logistic regression models, with normal persons who were treated as the reference category. In the age-adjusted model, age was entered as a variable of 10-year interval categories (50-59, 60-69, 70+ vs. <50). The multivariable model was further adjusted for major CVD risk factors including hypertension (yes vs. no), HDL-C (≥40 mg/dl vs. <40 mg/dl for males, ≥50 mg/dl vs. <50 mg/dl for females), triglycerides (150-199 mg/dl, ≥200 mg/dl vs. <150 mg/dl), BMI (20-22.9, 23-24.9, 25-27.9, 28-29.9, 30+ vs. <20), drinking habit (≤4 times/week, ≥5 times/week vs. non-drinkers), and smoking habit (ex-smokers, current smokers vs. non-smokers).

Results

Characteristics of study participants by age and gender were described in our previous paper [14]. We briefly summarize those here. It is observed that the averages for both systolic and diastolic blood pressure and CAVI scores linearly increased as ages advanced in both genders with an exception of men’s diastolic blood pressure after 60 years of age of which averages slightly decreased. Almost all averages of clinical indicators increased until 60 years of age except HDL-C of which averages were at the same level in both genders. Striking differences in averages of clinical indicators between genders are observed and are unfavorable for men in terms of cardiovascular disease risk. BMI averages ranged from 22.3 to 24.2 among men and from 20.4 to 22.9 among women. Greater prevalence of drinkers and smokers was observed in men than in women. Prevalence of drinkers and smokers was greater in younger women than in older women.

Table 1 shows characteristics of study participants by diabetes status. Averages of all variables in both genders were lowest in normal, highest in diabetes, and middle in prediabetes, except averages of HDL-C which were in the reverse order. Prevalence of abnormally high CAVI scores was lowest in normal, highest in diabetes, and middle in prediabetes.

Table 1 Characteristics of study participants by diabetes status

Prevalence of abnormally high CAVI scores is shown according to the status of prediabetes and diabetes by age and gender in Table 2. It is observed that such prevalence was higher after 40 years of age among persons with diabetes than either among persons with prediabetes or among normal persons in both genders.

Table 2 Prevalence (%) of abnormally high CAVI scores by status of prediabetes and diabetes mellitus

Table 3 shows odds ratios (ORs) of abnormally high CAVI scores in association with prediabetes and diabetes as compared with the reference of normal CAVI scores: (1) without an adjustment for confounding factors (crude ORs); (2) age breakdowns were added to logistic regression analysis; and (3) other CVD risk factors were further added to logistic regression analysis. As more confounding factors were adjusted, odds ratios decreased in both prediabetes and diabetes mellitus except ORs of prediabetes in women indicating that ORs after adjusting for ages and for CVD risk factors were almost identical and lower than crude OR. After adjusting for major CVD risk factors, significantly elevated ORs appeared among persons with prediabetes: 1.29 (95 % confidence interval (CI), 1.11-1.48) for men and 1.14 (CI, 1.01-1.28) for women, and among persons with diabetes: 2.41 (CI, 1.97-2.95) for men and 2.52 (CI, 1.94-3.28) for women.

Table 3 Estimated risk of having abnormally high CAVI scores in association with prediabetes and diabetes mellitus among Japanese urban workers and their families

Discussion

In the present cross-sectional study of Japanese urban workers and their families, prevalence of abnormally high CAVI scores were dose-dependently elevated along with the advancing degree of diabetic status both in men and in women. Compared with normal persons, significantly elevated odds ratios of abnormally high CAVI scores were observed not only in persons with diabetes but also in those with prediabetes.

With regard to the validity to use CAVI scores as an indicator of arteriosclerosis, Otsuka examined 72 deceased patients’ antemortem PWV (which is a basis for deriving CAVI scores) and pathological changes measured by the diffuse fibrotic thickening, formation of atheroma and calcification in the wall of their aorta. He reported multiple regression coefficient R = 0.810 between PWV and scores of those pathological changes [17]. In addition, other researchers reported that CAVI scores were significantly associated with coronary atherosclerosis [11], with carotid intima-media thickness and with homocysteine [12]. Thus, the use of CAVI scores derived from PWV values is valid to estimate the extent of arteriosclerosis.

VaSera VS-1000, which was used in our study, was designed to measure CAVI scores independent of blood pressure at the time of CAVI score measurement and CAVI scores represent the extent of arteriosclerosis between the aortic valve and the ankle. We have shown biological aging of the major artery by measuring CAVI scores in the CVD risk-free group and disease-related pathological aging of the major artery in the CVD high-risk group [13]. CAVI scores allow us to evaluate the extent of arteriosclerosis in the major arteries between the aortic valve and the ankle, to screen persons with subclinical stage of CVD, and provide an opportunity to modify diet and lifestyle to improve CAVI scores as reported by Satoh et al [18]. Thus, the use of CAVI scores potentially leads to savings on high treatment costs and to prolonging many productive lives.

Namekata et al. showed significant odds ratios of the abnormally high aortic pulse wave velocity (PWV, an indicator of arteriosclerosis reflecting stiffness of artery and atherosclerosis) in an association with diabetes among Japanese Americans and among native Japanese (3.66 and 2.43, respectively) [8]. Namekata et al. observed that odds ratios of abnormally high CAVI scores in an association with diabetes after making adjustment for ages were 10.02 for men and 8.42 for women, compared with those who did not have diabetes [14]. In our recent study significant ORs of prediabetes and diabetes were observed in association with most CVD risk factors including age, hypertension, triglycerides ≥ 200 mg/dl and BMI ≥ 25 in both genders [19]. Our results imply much faster advancement of arteriosclerosis among persons both with prediabetes and with diabetes than in the normal persons.

In our present study it is shown that prevalence of abnormally high CAVI scores has an increasing trend with an age increase in all three groups of non-diabetes normal persons, prediabetes and diabetes (Table 2). Our results show that comparing with persons with normal CAVI scores, estimated risks of abnormally high CAVI scores in association with prediabetes and with diabetes were 1.29 times higher and 2.41 times higher, respectively, in men and were 1.14 times higher and 2.52 times higher, respectively, among women (Table 3). Our study results have confirmed the association between prediabetes and arteriosclerosis and the much stronger association between diabetes and arteriosclerosis estimated by CAVI scores.

A limitation of this study is that it is an observational and cross-sectional study. Thus, we cannot tell how many persons with prediabetes will develop diabetes or will return to normal in their blood glucose levels or/and HbA1c levels. The strengths of this study are having the large sample size and being able to adjust for many CVD risk factors with enough statistical power to examine the association of CAVI scores with prediabetes and diabetes mellitus.

Conclusion

In conclusion, both prediabetes and diabetes mellitus were significantly associated with CAVI scores. It is implied that the extent of arteriosclerosis (including arterial stiffness and atherosclerosis) was moderately enhanced among persons with prediabetes and was further advanced among persons with diabetes. Thus, it is important to introduce earlier interventions for changing lifestyle and diet of persons with prediabetes in order to prevent them from developing diabetes and further advancing arteriosclerosis.

Abbreviations

BMI:

Body mass index

CAVI:

Cardio-ankle vascular index

CHD:

Coronary heart disease

CI:

Confidence interval

CVD:

Cardiovascular disease

DBP:

Diastolic blood pressure

Diabetes:

Diabetes mellitus

HbA1c:

Glyco-hemoglobin A1c

HDL-C:

High density lipoprotein cholesterol

ORs:

Odds ratios

PWV:

Pulse wave velocity

SBP:

Systolic blood pressure

SD:

Standard deviation

TC:

Total cholesterol

TG:

Triglycerides

References

  1. Bramwell JC, Hill AV. The velocity of the pulse wave in man. Proc Roy Soc of London Series B. 1922;93:298–306.

    Article  Google Scholar 

  2. Bramwell JC, Hill AV. Velocity of transmission of the pulse wave and elasticity of arteries. Lancet. 1922;202:891–2.

    Article  Google Scholar 

  3. Bramwell JC, Hill AV, McSwiney BA. Velocity of transmission of the pulse wave in man as related to age as measured by the hot-wire sphygmograph. Heart. 1923;10:233–49.

    Google Scholar 

  4. Bramwell JC, McDowall RJS, McSwiney BA. The variation of the arterial elasticity with blood pressure in man. Proc Roy Soc of London Series B. 1923;94:450–4.

    Article  Google Scholar 

  5. Bramwell JC, Dowing AC, Hill AV. The effect of blood pressure on the extensibility of the human artery. Heart. 1923;10:289–300.

    Google Scholar 

  6. Hasegawa M. Fundamental research on human aortic pulse wave velocity. Jikei Medical Journal. 1970;85:742–60.

    Google Scholar 

  7. Namekata T, Moore D, Suzuki K, Mori M, Hatano S, Hayashi C, Abe N, Hasegawa M. A study of the association between the aortic pulse wave velocity and atherosclerotic risk factors among Japanese Americans in Seattle, U.S.A. Jpn J Pub Health. 1997;44:942–51. http://www.ncbi.nlm.nih.gov/pubmed/9553384.

  8. Namekata T, Suzuki K, Arai C. Seattle Nikkei health study: Cross cultural surveys between Seattle and Japan. In: Shigematsu K, Okada A, Imaizumi T, Hoshino T, editors. New Trends in Psychometrics. Tokyo: Universal Academy Press, Inc; 2008. p. 339–46.

    Google Scholar 

  9. Shirai K, Utino J, Otsuka K, Takata M. A novel blood pressure-independent arterial wall stiffness parameter; cardio-ankle vascular index (CAVI). J Atherosclerosis Thrombosis. 2006;13:101–7. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3636518/.

    Article  Google Scholar 

  10. Hayashi G, Sato M, Niimi H, Handa H, Moritake K, Okumura A. Analysis of vascular wall constitutive law with finite deformation theory. Medical Electronics & Biological Engineering. 1975;13:293–7.

    CAS  Google Scholar 

  11. Nakamura K, Tomaru T, Yamamura S, Miyashita Y, Shirai K, Noike H. Cardio-ankle vascular index is a candidate predictor of coronary atherosclerosis. Circ J. 2008;72:598–604. http://www.ncbi.nlm.nih.gov/pubmed/18362432.

    Article  PubMed  Google Scholar 

  12. Kadota K, Takamura N, Aoyagi K, Yamasaki H, Usa T, Nakazato M, Maeda T, Wada M, Nakashima K, Abe K, Takeshima F, Ozono Y. Availability of cardio-ankle vascular index (CAVI) as a screening tool for atherosclerosis. Circ J. 2008;72:304–8. http://www.ncbi.nlm.nih.gov/pubmed/18219171.

  13. Namekata T, Suzuki K, Ishizuka N, Shirai K. Establishing baseline criteria of cardio-ankle vascular index as a new indicator of arteriosclerosis: a cross-sectional study. BMC Cardiovasc Disord. 2011; 10.1186/1471-2261-11-51 http://www.ncbi.nlm.nih.gov/pubmed/21831311

  14. Namekata T, Suzuki K, Ishizuka N, Nakata M, Shirai K. 2012: Association of Cardio-Ankle Vascular Index with Cardiovascular Disease Risk Factors and Coronary Heart Disease among Japanese Urban Workers and their Families. J Clinic Experiment Cardioldoi:10.4172/2155-9880.S1-003 http://www.omicsonline.org/association-of-cardio-ankle-vascular-index-with-cardiovascular-disease-risk-factors-and-coronary-heart-disease-among-japanese-urban-workers-and-their-families-2155-9880.S1-003.php?aid=3820

  15. American Diabetes Association. Standards of Medical Care in Diabetes- 2012. Diabetes Care. 2012;35:S11–63. http://care.diabetesjournals.org/content/35/Supplement_1/S11.full.pdf+html.

    Article  Google Scholar 

  16. Rosendorff C, Black HR, Cannon CP, Gersh BJ, Gore J, et al. Treatment of hypertension in the prevention and management of ischemic heart disease: a scientific statement from American Heart Association Council for High Blood Pressure Research and the Councils on Clinical Cardiology and Epidemiology and Prevention. Circulation. 2007;115:2761–88. http://www.ncbi.nlm.nih.gov/pubmed/17502569.

    Article  PubMed  Google Scholar 

  17. Otsuka F. A study of the relationship between pulse wave velocity of human aorta and its postmortem histopathology. Tokyo Jikei Medical College Journal. 1973;88:11–6.

    Google Scholar 

  18. Satoh N, Shimatsu A, Kato Y, Araki R, Koyama K, Okajima T, Tanabe M, Ooishi M, Kotani K, Ogawa Y. Evaluation of the cardio-ankle vascular index, a new indicator of arterial stiffness independent of blood pressure, in obesity and metabolic syndrome. Hypertens Res. 2008;31(10):1921–30. http://www.ncbi.nlm.nih.gov/pubmed/19015600.

  19. Namekata T, Shirai K, Nakata M, Suzuki K, Arai C, Ishizuka N. Association of prediabetes and diabetes mellitus with cardiovascular disease risk factors among Japanese urban workers and their families: A cross- sectional study. Epidemiology Open Access. 2014;4:157. http://www.omicsonline.org/open-access/association-of-prediabetes-and-diabetes-mellitus-with-cardiovascular-disease-risk-factors-among-japanese-urban-workers-and-their-families-a-cross-sectional-study-2161-1165.1000157.pdf.

    Article  Google Scholar 

Download references

Acknowledgment

We express our sincere appreciation to Japan Health Promotion Foundation for their providing us with their data set for this study. The study was in part supported by the grants from Japan Health Promotion Foundation and SR Foundation.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Tsukasa Namekata.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TN, KoS, KeS, NT, and KM participated in the design of the study. KeS, CA, and NI managed the CVD prevention screening and organized the data set. TN and MN performed statistical analyses. TN drafted the manuscript and revised based on other authors’ comments. All other authors reviewed critically and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Namekata, T., Shirai, K., Tanabe, N. et al. Estimating the extent of subclinical arteriosclerosis of persons with prediabetes and diabetes mellitus among Japanese urban workers and their families: a cross-sectional study. BMC Cardiovasc Disord 16, 52 (2016). https://doi.org/10.1186/s12872-016-0230-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12872-016-0230-6

Keywords