*Corresponding Author:
P. Limpawattana
Division of Geriatric Medicine, Department of Internal Medicine, Faculty of Medicine, Khon Kaen University, Thailand
E-mail: lpanit@kku.ac.th
This article was originally published in a special issue, “Biomedical applications in Pharmaceutical Sciences”
Indian J Pharm Sci 2020:82(2)Spl issue3;93-100

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License, which allows others to remix, tweak, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms


The objectives of this study were to determine the distributions of weight-adjusted bone mass content and osteoporosis in chronic obstructive pulmonary disease patients, to examine the prevalence of osteoporosis and osteopenia and to determine factors associated with weight-adjusted bone mass content. A crosssectional study of chronic obstructive pulmonary disease patients at a tertiary-care hospital, Thailand was conducted from May 2015 to December 2016. Baseline characteristics were collected and bone mineral density and weight-adjusted bone mass content were measured using dual-energy X-ray absorptiometry. A total of 107 chronic obstructive pulmonary disease patients were included. The prevalence of osteoporosis and osteopenia were 35.5 and 38.2 %, respectively. The distribution of weight-adjusted bone mass content of arms, legs and whole body was significantly correlated with chronic obstructive pulmonary disease severity (p<0.05), but not with the presence of osteoporosis. Factors associated with weight-adjusted bone mass content of the extremities was the severity of chronic obstructive pulmonary disease (adjusted odds ratios of moderate and severe chronic obstructive pulmonary disease were 1.03 and 1.06), waist circumference (adjusted odds ratio 0.99) and appendicular skeletal mass index (adjusted odds ratio 0.98). Both osteoporosis and osteopenia were present in more than a third of the chronic obstructive pulmonary disease patients. The weight-adjusted bone mass content of arms, legs, extremities and whole body were correlated with the degree of airflow limitation but not the presence of osteoporosis. Factors associated with weight-adjusted bone mass content of the extremities were identified. Weight-adjusted bone mass content could be a good index to indicate the degree of bone loss in chronic obstructive pulmonary disease patients.


Bone mass content, Chronic lung disease, Osteopenia, Osteoporosis

Chronic obstructive pulmonary disease (COPD) is a common, preventable and treatable disease that is known mainly to affect the structure and function of the lungs. It is characterized by persistent respiratory symptoms and airflow limitation due to airway and/or alveolar abnormalities, usually caused by significant exposure to noxious particles or gases[1]. It is a major cause of morbidity and premature death and is the third leading cause of death worldwide[2]. According to a systematic analysis for the Global Burden of Disease Study conducted in 2015, the prevalence of COPD increased by 44.2 % from 1990-2015 and there were 3.2 million people who died from COPD globally in 2015 with an increase of 11.6 % compared to 1990[3]. Furthermore, COPD is a systemic disease that is complicated by numerous comorbidities and syndromes including lung cancer, atherosclerosis, weight loss and nutritional abnormalities, sarcopenia, frailty syndrome, hormonal and metabolic disturbances, osteoporosis, anxiety and depression[4-6]. Therefore, management of COPD should be directed toward identifying and treating those systemic effects, in addition to focusing on airway limitations, in order to improve health outcomes in COPD patients[6].

Osteoporosis is one of the systemic effects of COPD, though the causal relationship and molecular association between COPD and osteoporosis has yet to be established[4]. It is commonly found in COPD patients and its prevalence is higher in these patients than in healthy subjects and those with various other chronic lung diseases such as asthma, idiopathic pulmonary fibrosis, pulmonary hypertension, sarcoidosis, scleroderma or Kartagener’s syndrome and bronchiectasis[7-11]. Osteoporosis is also more prevalent in women than in men[11]. Moreover, forced expiratory volume in one second (FEV1 % predict), body mass index (BMI) and free fat mass index (FFMI) have been found to be lower in osteoporotic COPD patients[11]. The reported prevalence of osteoporosis and osteopenia in COPD patients have varied from 8.7 to 69 % and 18.6 to 67 %, respectively[4,6,7,10-17]. Osteoporosis can lead an increased risk of fracture, worsen pulmonary function and impair the activities of daily living in COPD patients. Therefore, the concurrence of COPD and osteoporosis could lead to a vicious cycle, causing a substantial burden to the patients[4,13]. COPD-related risk factors for osteoporosis include low body weight, smoking, age, gender, inactivity, systemic inflammation, pulmonary dysfunction, systemic corticosteroids, and vitamin D deficiency or insufficiency[13,14].

Routine screening and management for osteoporosis in COPD patients enables physicians to diagnose this disorder at an early phase to reduce its negative consequences. However, some studies only have found the degree of reduced bone mineral density (BMD) to be comparable to the severity of airflow limitation[7-10], meaning that measurement of BMD might not be a good index or sensitive enough tool to determine the degree of bone loss in COPD patients. Additionally, there have been few studies conducted in Thailand regarding the prevalence of osteoporosis and osteopenia and their association with COPD severity. One study reported that the prevalence’s of osteoporosis and osteopenia were 31.4 and 32.4 %, respectively and that they were not associated with severity of COPD[16]. Measurement of bone mass content that included parts of the body other than the femoral neck and lumbar spine like BMD might be helpful. One report showed that bone mass content index (BMCI), which is derived from bone mass content/height2, was lower in the legs, trunk and overall in COPD patients than in control subjects. The BMCI of the trunk was associated with VO2 max and BMI, but not with airflow limitation. BMI, which includes body weight, appears to be more closely related to bone mass content than height alone. Therefore, the purpose of this study was to determine the distribution of weightadjusted bone mass content (w-BMC) and osteoporosis in male patients with COPD by COPD severity. This study also aimed to examine the prevalence of osteoporosis and osteopenia and factors that affect w-BMC.

Materials and Methods

Subjects and setting

This was a sub-study of several previously published articles[5,18]. A cross-sectional study of COPD patients was conducted at the Srinagarind Hospital COPD clinic at Khon Kaen University’s Faculty of Medicine, Thailand from May 2015 to December 2016. The clinic is a tertiary care referral center at the university hospital to which patients are sent from the internal medicine outpatient clinic and other community hospitals. The inclusion criteria were, male sex, Thai citizenship, COPD, age over 18 y and attendance at the Srinagarind Hospital COPD clinic. Subjects were excluded if they had other active illnesses such as functional class III-IV heart failure, pneumonia or septicemia, those who had experienced COPD exacerbation within the preceding 4 w, those who were unable to complete the study due to physical limitations such as being bed ridden or suffering from hemiparalysis, those who were unwilling to participate in the study, those who were limited in their ability to undergo dual energy X-ray absorptiometry (DXA) or had conditions that would affect DXA results, those who had taken barium or had had an enema within the previous 2 w, those who had metallic instrumentation such as a vertebral fixation device, had undergone a radionuclide scan within the previous 2 w, or had body weight over 200 kg exceeding the capacity of the DXA. Subjects with incomplete body composition data were also excluded.


COPD was diagnosed in subjects according to the criteria laid out by the Global Initiative for Chronic Obstructive Lung Disease (GOLD): persistent airflow limitation and a post-bronchodilator force expiratory volume in one second (FEV1) less than 70 % of forced vital capacity (FVC). Severity of COPD was based on the GOLD staging system, which consists of 4 levels, mild (FEV1> 80 % normal), moderate (FEV1 50 -79 % normal), severe (FEV1 30-49 % normal) and very severe (FEV1<30 % normal)[1].

The MMRC scale was used to evaluate dyspnea in daily life. The patients were asked to rate their dyspnea on a scale of 0-4 for various physical activities that can provoke the condition as follows: 0, not troubled with breathlessness except with strenuous exercise; 1, troubled by shortness of breath when hurrying on the level or walking up a slight hill; 2, walks slower than people of the same age on the level because of breathlessness or has to stop for breath when walking at own pace on the level; 3, stops for breath after walking about 100 yards or after a few min on the level; and 4, too breathless to leave the house, or breathless when dressing or undressing[19].

DXA (General Electric Lunar-Progidy) was used to estimate bone mass and skeletal muscle mass in patients diagnosed with osteoporosis whose T-scores of the femoral neck or lumbar spine were <-2.5 standard deviation[20]. Weight-adjusted bone mass content (w-BMC) was calculated as bone mass content/weight2 (g/kg2). Appendicular skeletal muscle mass index (ASMI) was calculated as appendicular skeletal mass/ height2 (kg/m2)[21].

Handgrip strength (kg) was measured using a handheld dynamometer (GRIP-D Model: T.K.K.5401). The patients were asked to hold the handgrip dynamometer in one hand with the arm at a right angle and the elbow by the side of the body. When ready, the patients squeezed the dynamometer with maximum isometric effort, which was maintained for about 5 s. Three attempts were measured and the maximal measurement was recorded[22]. The patient’s 6 min walk distance (6MWD) was used to measure physical performance according to the American Thoracic Society guidelines[23]. Gait speed (m/s) was then analyzed based on 6MWD.


Potential subjects were asked by the team of researchers to participate in this study at the COPD clinic. They were informed of the research purpose, protocol, risks and benefits of enrolling. After those who met the inclusion criteria and were willing to participate in the study gave informed consent, their baseline information was collected by trained clinical researchers. The demographic data collected included age, sex, educational level, smoking status, severity of COPD using the most recent pulmonary function test results (within the past 6 mo), MMRC score, comorbid diseases, regular medication use, history of falls over the previous 12 mo, self-reported baseline function with regard to basic and instrumental activities of daily living (ADL) using the Barthel ADL index and Chula index[24], history of non-elective admission over the previous 12 mo, weight, height, waist circumference and systolic and diastolic blood pressure at rest. Weight and height were used to calculate BMI. Waist circumference (WC) was measured after exhalation midway between the lower rib and the iliac crest. Systolic and diastolic blood pressure was presented as mean arterial blood pressure (MAP). Handgrip strength was measured using a grip dynamometer[22], and 6MWD (as previously described) was used to measure physical performance[21,24]. DXA was used to measure body composition including muscle mass, bone mineral density and bone mass content. All measurements were performed on all patients during the same period. The main purpose of DXA is to measure BMD. Using spectral imaging but not the body composition, perhaps DXA may be used to measure body composition only to a certain extent in a few cases.

Statistical analysis:

Baseline data were analyzed using descriptive statistics and were presented as percentage, mean and standard deviation (SD). If the distributions of these data were not normal, then median and inter-quartile range were used instead. The Chi-squared test/Fisher’s exact test was used to compare the proportions of patients with and without osteoporosis between the COPD patients. Differences in w-BMC by COPD severity were determined using a one-way analysis of variance. A p-value <0.05 indicated statistical significance. Effects of factors associated with w-BMC of the extremities were evaluated using univariate and stepwise multivariate linear regression analyses with logistic transformation. In terms of univariate analysis, crude odds ratio (ORs) and 95 % confidence intervals (CIs) were used to determine the strength of association of the factors examined with w-BMC of the extremities. Factors with a p-value of <0.20 or that were found to have clinical significance in previous studies were then entered into a multiple linear regression model. A p-value of <0.05 was considered to indicate statistically significant differences and adjusted OR (AOR) and 95 % CIs were used to determine the strength of association. All data analysis was carried out using STATA version10.0 (StataCorp, College Station, Texas). Ethical approval was provided by the Khon Kaen University Faculty of Medicine Ethics Committee as instituted by the Helsinki Declaration (Number HE581060).

Results and Discussion

Patient characteristics are presented in Table 1. Of the 109 subjects evaluated, 2 were excluded because of incomplete body composition data, which left 107 patients to be enrolled in this study. The mean age was 69.7 y. The majority of the patients had moderate COPD (59.8 %), and inhaled corticosteroids were the most commonly used medication (94.4 %). Most of the subjects were independent with regard to activities of daily living. Osteoporosis was found in about a third (35.5 %) of the subjects.

Variables N=107
Age (y), mean (SD) 69.7 (9.1)
Age (y), n (%) < 65 y 29 (27.1)
65-74 y 43 (40.2)
> 75 y 35 (32.7)
Educational level, n (%) 6 y or less 50 (46.7)
7-12 y 26 (24.3)
>12 y 31 (29.0)
Current smoker, n (%) 8 (7.5)
Severity of COPD, n (%) Mild 26 (24.3)
Moderate 64 (59.8)
Severe 16 (15.0)
Very severe 1 (0.9)
MMRC score, mean (SD) 0.7 (1.1)
Comorbid diseases, n (%) Diabetes mellitus 12 (11.2)
Hypertension 39 (36.5)
Dyslipidemia 10 (9.4)
Chronic arthritis 9 (8.4)
Cancer 6(5.6)
Regular medication, n (%) Inhaled steroids 101 (94.4)
Systemic steroids 2 (1.9)
Oral hypoglycemic drugs 6 (5.6)
Statin 11 (10.3)
NSAIDs 1 (0.9)
Twoor more falls over the past 1 y, n (%) 3 (2.8)
Barthel scores, mean (SD) 19.9 1.1)
Chula iADLs scores, mean (SD) 9 (0.2)
Non-elective admission over the past 1 y, n (%) 28 (26.2)
BMI (kg/m2), mean (SD) 23.2 (4.0)
WC (cm), mean (SD) 89.6 11.4)
MAP (mm Hg), mean (SD) 96.1 (10.8)
Appendicular skeletal mass index (kg/m2), mean (SD) 7.2 (0.9)
6MWD (m), mean (SD) 426.3 (76.2)
Gait speed (m/s), mean (SD) 1.2 (0.2)
Handgrip strength (kg), mean (SD) 26.9 5.9)
Osteopenia, n (%) 41 (38.2)
Osteoporosis, n (%) 38 (35.5)
w-BMC (g/kg2), mean (SD) Arms 0.1 (0.03)
Legs 0.25 (0.08)
Extremities 0.34 (0.11)
Trunk 0.18 (0.05)
Whole body 0.65 (0.20)

Table 1: Baseline Characteristics of the Studied Population

The distribution of osteoporosis in patients with COPD by severity is shown in fig. 1. The prevalence of osteoporosis did not differ significantly with COPD severity, while the w-BMC of the arms, legs, extremities, and overall did (p<0.05), as shown in fig. 2.


Figure 1: The distribution of osteoporosis (%) in patients with COPD by COPD severity
P<0.05 were considered to indicate significance. No- no osteoporosis, Yes- presence of osteoporosis, osteoporosis is defined as a femoral neck or lumbar spine T-score of 2.5 standard deviation or more below the young adult mean


Figure 2: The distribution of the weight adjusted bone mass content The distribution of the w-BMC (g/cm2) in the arms, legs, extremities, trunk and whole body in COPD patients by COPD severity *P<0.05 were considered to indicate significance

w-BMC of the extremities was used as an index of bone loss among COPD patients in this study because w-BMC of the arms, legs, extremities and overall differed significantly by level of COPD severity (fig. 2), which implied that w-BMC of the arms and legs were important regions in predicting bone loss. Thus, w-BMC of the extremities was used to analyze the associations of other factors. Severity of COPD, presence of diabetes, WC, MAP, ASMI, and handgrip strength were associated with w-BMC of the extremities according to univariate analysis as shown in Table 2. After multicollinearity was confirmed, factors with a p<0.2 according to univariate analysis and those previously reported to be clinically significant (including age, smoking status, severity of COPD, MMRC, presence of DM, dyslipidemia, chronic arthritis, regular use of inhaled or systemic steroids, NSAIDs, and statin, fall history, MAP, ASMI, waist circumference, and handgrip strength) were entered into the model for multivariate analysis. These factors were then entered in the stepwise multiple regression models with logistic transformation. Greater severity of COPD, decreased WC and decreased ASMI were factors associated with increased w-BMC in this study (Table 2).

Factors Univariate analysis Multivariate analysis
Crude OR (95 % CI) p-value Adjusted OR (95 % CI) p-value
Age 1.00 (0.99,1.001) 0.09
Educational level
· 6 years or lower 1 - -
· 7-12 years 0.98 (0.94,1.04) 0.72
· > 12 years 0.99 (0.93,1.03) 0.34
Current smoker 0.96 (0.89,1.04) 0.31
Severity of COPD
· Mild 1 - - 1 - -
· Moderate 1.01 (0.97,1.06) 0.56 1.03 (1.003,1.05) 0.03*
· Severe 1.1 (1.04,1.17) 0.003* 1.06 (1.02,1.09) 0.003*
· Very severe 1.35 (1.11,1.65) 0.003* 1.15 (1.03,1.28) 0.02*
MMRC 1.02 (0.99,1.04) 0.05 1.01 (0.99,1.02) 0.19
Diabetes mellitus 0.94 (0.88,0.99) 0.04*
Hypertension 0.98 (0.93,1.02) 0.24
Dyslipidemia 0.94 (088,1.01) 0.11
Chronic arthritis 0.94 (0.87,1.01) 0.09
Cancer 1.03 (0.94,1.13) 0.48
Regular medication
· Inhaled corticosteroids 0.96 (0.87,1.04) 0.31
· Systemic steroids 1 (0.86,1.16) 0.99 0.93 (0.86,1.0.1) 0.06
· Oral hypoglycemic drugs 0.97 (0.88,1.06) 0.45
· Statin 0.95 (0.86,1.01) 0.1 0.93 (0.94,1.01) 0.2
· NSAIDs 0.86 (0.70,1.07) 0.17
Two or more falls over the past 12 months 1.02 (0.89,1.15) 0.78
Barthel scores 1 (0.98,1.02) 0.73
Chula iADLs scores 1.04 (0.93,1.15) 0.47
WC (cm) 0.992 (0.991,0.993) 0.00* 0.993 (0.992,0.994) 0.00*
Non-elective admission over the past 12 months 1.04 (0.99,1.08) 0.09
MAP 0.997 (0.995,0.999) 0.01*
Appendicular skeletal mass index (kg/m2) 0.92 (0.91,0.94) 0.00* 0.98 (0.97,0.99) 0.04*
Gait speed (m/s) 0.99 (0.89,1.09) 0.84
Handgrip strength (kg) 0.993 (0.991,0.996) 0.00* 0.99 (0.99,1.00) 0.14

Table 2: Univariate and Stepwise Multivariate Linear Regression Analyses

The prevalence of osteoporosis and osteopenia among male patients with COPD in this study were 35.5 and 38.2 %, respectively, which is in line with the findings of previous studies, in which these prevalence ranged from 8.7-69 % and 18.6-67 %, respectively[4,6,7,10-17]. Variations in prevalence among studies might be attributable to methodological differences in the assessment of BMD and differences in the characteristics of study populations such as age, sex, ethnicity, severity of COPD and past use of bone medications. Another study conducted in Thai COPD patients found that 31.4 % had osteoporosis and 32.4 % had osteopenia, similar to the findings of this study[16]. The high prevalence of these conditions in COPD patients highlights the importance of early recognition of the risk factors and implementation of strategies to manage osteoporosis in these patients in order to prevent adverse health outcomes such as osteoporotic fracture and poor quality of life[7]. However, this study could not demonstrate the relationship between severity of COPD and the occurrence of osteoporosis found by a study conducted in Italy[9] among others[7,15]. This is partly due to the small number of patients in our study with greater COPD severity.

Rampant use of corticosteroids through inhalational, systemic routes appeared to have blocked the absorption of calcium into systemic circulation eventually causing decrease in the BMD, health of the soft muscle tissues. One interesting finding of this study was the distribution of the w-BMC (g/cm2) in the arms, legs, extremities and whole body in COPD patients, which demonstrated that higher w-BMC was significantly related to greater severity of COPD. Thus, this index might be used to predict the severity of bone loss in patients with COPD with greater accuracy than the presence of osteoporosis, which is based on the BMD of the femoral neck or lumbar spine. A possible explanation for this association is the connection between body weight and BMC found in several previous studies[9,25-27]. Low body weight might be associated with lower calcium intake, protein malnutrition, immobilization, and corticosteroid administration and it can affect the biomechanical stress placed on bones[26,27]. Body weight must thus be considered when examining changes in BMC[26]. Furthermore, body weight is likely to decline to a greater degree than the bone mass in COPD patients[9], meaning that higher w-BMC is associated with greater severity of airflow limitation. A prior report also showed that there was no significant association between height-adjusted BMC and severity of COPD[25], which suggests that height is not a suitable as a parameter to examine changes in BMC. In our study, only w-BMC of the trunk was not significantly related with severity of COPD. This can be explained by the mechanism of BMC preservation in the mechanical stress on the bones like trunk. COPD patients tend to have higher metabolic and ventilatory requirements and dyspnea than non-COPD patients, so the bone mass content in the trunk is inclined to remain stable while the limitations to physical activity, such as walking and standing, in patients with severe COPD leads to less mechanical stress on the bones, resulting in significant changes in the w-BMC of the arms, legs, and overall[25].

Independent factors associated with high w-BMC of the extremities in this study were greater severity of COPD, decreased WC, and ASMI. The association between w-BMC and COPD severity may be explained by the bone loss and decreases in body weight that result from a greater degree of airflow limitation. However, body weight is reduced to a greater degree than bone mass in these cases, as mentioned earlier[9,25-27], meaning that the proportion of bone mass content to body weight2 tends to be higher in more severe cases[6,9,25]. Decreased WC is another factor that was related to high w-BMC, which was a marker of bone loss in COPD patients in this study. This confirms that WC is positively associated with bone mass[28] and supports a prior finding that WC is a key determinant of bone mass[29]. However, this relationship has been found to hold only up to a certain body size threshold and not to be significant in subjects with WCs of 100-110 cm or greater[28,29]. The mean WC of the subjects in this study was below this threshold, at 89.6 cm, indicating that there should be a significant correlation between WC and bone mass in this population. Generally, WC is an indicator of abdominal obesity, which represented visceral adipose tissue[29,30], and several reports have found this surrogate measure to be associated with low bone mass and increased risk of hip fracture when body mass index was taken into account[29,31,32]. This implies that both abdominal obesity and lower waist circumference are associated with bone mass loss. Decreases in ASMI have also been shown to be independently associated with w-BMC. The results presented here support the oppositional relationship between lean mass and bone mass loss[28]. In men, increases in bone and muscle mass result from greater levels of testosterone and IGF-1[33]. In patients with COPD, there are additional cofactors that exacerbate the decline of skeletal mass and bone including (1) physical inactivity due to exercise intolerance, (2) systemic inflammation causing excessive apoptosis of skeletal muscle and increased inflammatory cytokines, such as tumor necrosis factor alpha and IL-6, which contribute to the regulation of bone remodeling, and (3) the possible presence of hypoxia and use of systemic steroids[13,18,33,34].

As osteoporosis is also common in COPD patients and results in various unfavourable health outcomes, early prevention and treatment is crucial. However, this study found that a diagnosis of osteoporosis based on femoral or lumbar bone mineral density was not related to severity of COPD and that w-BMC had several advantages over the traditional diagnosis of osteoporosis as an indicator of bone loss. First, it correlated to the severity of airflow limitation. Second, it includes bone mass content of all extremities, where the traditional method includes only the femoral neck and lumbar spine. Third, it is adjusted by body weight, the effect of which on the occurrence of osteoporosis and COPD is supported by an abundance of data.

There were some limitations to this study. First, the number in the group with COPD severity stage “Very severe” is small, further study with higher number of subjects should be conduct. Second, the subjects were all Asian men, so the result might not be representative of women or other racial or ethnic groups, as these factors can influence bone mineral density[2,35]. Third, some known factors that are associated with osteoporosis were not measured such as physical activity and vitamin D levels. Last, this study was conducted in a tertiary care hospital, and the subjects may have had more comorbidities or greater severity of disease than those in a community hospital setting. Therefore, the prevalence of osteoporosis and osteopenia found in this study might not be generalizable to other settings. Further research in various study populations is required.

The prevalence of osteoporosis and osteopenia was high in male patients with COPD. The traditional diagnosis of osteoporosis was not related to the severity of airflow limitation, where w-BMC of the arms, legs, extremities, and whole body was. Factors associated with w-BMC of the extremities were the severity of COPD, WC, and ASMI. We recommend using w-BMC of the extremities to indicate degree of bone loss as a part of routine COPD management. Further studies are required to investigate the optimal w-BMC cut-off point that is related to adverse health outcomes.


This study was funded by the Neuroscience Research and Development Group (Khon Kaen University, Thailand).

Conflict of interest:

All authors have no conflicts of interest with respect to the authorship.


We would like to thank Dr. Dylan Southard (Research Affairs, Faculty of Medicine, Khon Kaen University, Thailand) for editing the manuscript. The authors also would like to thank Sleep Apnea Research Group (Khon Kaen University, Thailand).

Authors’ contributions:

Dr Jarin Chindaprasirt contributed to the study design, concept of this work, interpretation of data, and approval for final version of the article. Dr Watchara Boonsawat, Dr Daris Theerakulpisut, Dr Manchumad Manjavong, Prof. Kittisak Sawanyawisuth, and Prof. Verajit Chotmongkol contributed to the concept of this study and approval for the final version of the article. Dr Pratchaya Inthasuwan and Dr Siraphong Putraveephong contributed to acquisition of data. Dr Panita Limpawattana contributed to the study design, analysis and interpretation of data, writing drafted manuscript and approval the version to be published. All authors have read and approved the manuscript.