Association of Body Mass Index and Fracture Risk Varied by Affected Bones in Patients with Diabetes: A Nationwide Cohort Study
Article information
Abstract
Background
Body mass index (BMI) is a risk factor for the type 2 diabetes (T2DM), and T2DM accompanies various complications, such as fractures. We investigated the effects of BMI and T2DM on fracture risk and analyzed whether the association varied with fracture locations.
Methods
This study is a nationwide population-based cohort study that included all people with T2DM (n=2,746,078) who received the National Screening Program during 2009–2012. According to the anatomical location of the fracture, the incidence rate and hazard ratio (HR) were analyzed by dividing it into four categories: vertebra, hip, limbs, and total fracture.
Results
The total fracture had higher HR in the underweight group (HR, 1.268; 95% CI, 1.228 to 1.309) and lower HR in the obese group (HR, 0.891; 95% CI, 0.882 to 0.901) and the morbidly obese group (HR, 0.873; 95% CI, 0.857 to 0.89), compared to reference (normal BMI group). Similar trends were observed for HR of vertebra fracture. The risk of hip fracture was most prominent, the risk of hip fracture increased in the underweight group (HR, 1.896; 95% CI, 1.178 to 2.021) and decreased in the obesity (HR, 0.643; 95% CI, 0.624 to 0.663) and morbidly obesity group (HR, 0.627; 95% CI, 0.591 to 0.665). Lastly, fracture risk was least affected by BMI for limbs.
Conclusion
In T2DM patients, underweight tends to increase fracture risk, and overweight tends to lower fracture risk, but association between BMI and fracture risk varied depending on the affected bone lesions.
INTRODUCTION
Prevalence of diabetes has increased steadily during the last decades, leading to worldwide health concerns [1]. According to a World Health Organization (WHO) report [2], an estimated 422 million adults with diabetes were living in 2014, compared to 108 million in 1980. The age-standardized global prevalence of diabetes in 2014 has approximately doubled since 1980, rising from 4.7% to 8.5% in adults aged 18 years and older. The International Diabetes Federation (IDF) predicts the global number of diabetic patients and the prevalence to be 642 million and 10.4%, respectively in 2040 [3]. Diabetes being a substantial cause of mortality and morbidity worldwide has become a major public health problem. In 2010, there were approximately 1.3 million diabetes-related deaths worldwide, doubling compared to 1990 [4]. Diabetes is a metabolic disorder that greatly alters the body system affecting nearly all organs and resulting in various health complications.
Previous studies have shown that individuals with type 1 diabetes mellitus (T1DM) have lower bone mineral density and higher risk of hip fractures [5]. Although type 2 diabetes mellitus (T2DM) also has been reported to show increased risk of hip fractures, the results are controversial [6-9]. Recently, even prediabetic conditions have been reported to be associated with the occurrence of hip fractures [10]. Some studies suggested that various factors, including abnormal bone turnover, poor glycemic control and bone loss, may explain the increased risk of fractures in T2DM [11-14]. It is recommended by the WHO to use body mass index (BMI) as an important indicator for obesity [15]. A higher BMI level is associated with a higher prevalence of T2DM, which is known as an important risk for hip fractures [16,17].
Most studies have focused on hip fractures, and the fracture risk according to location of the skeletal lesion is not clear. Few studies have investigated the impact of BMI and T2DM in limited number of bones, and others have investigated whether the fracture risk differed depending on skeletal sites. In addition, the association between BMI and fracture risk in Asian populations have been reported, in whom obesity is not as extensive among diabetic patients as in Western countries [18, 19]. To fill in this gap, we first aimed to investigate the association between BMI and risk of fractures using a large nationwide study from health insurance claims data. We further analyzed whether the association varied by fracture location.
METHODS
Data source and study population
In this retrospective longitudinal nationwide population-based cohort study, we used the Korean National Health Insurance System (NHIS) database, a mandatory social medical insurance system run by the Korean government. Those enrolled in the health insurance service are recommended to receive health check-ups at least biennially. Among the patients who underwent health screenings from January 2009 to December 2012 we focused on patients with T2DM. At that time, enrolled diabetes mellitus (DM) patients were linked claim data and confirmed by medical details. In other words, definition of T2DM included the following criteria: (1) at least one claim per year under International Classification of Diseases, 10th Edition (ICD-10) codes E11–14 and at least one claim per year for the prescription of antidiabetic medication (sulfonylureas, metformin, meglitinides, thiazolidinediones, dipeptidyl peptidase-4 inhibitors, α-glucosidase inhibitors, or insulin) or (2) fasting glucose level ≥126 mg/dL. Patients who were admitted to the hospital more than once or visited the outpatient clinic more than twice were also included as T2DM subjects [20].
All subjects who underwent health screenings from January 2009 to December 2012 (n=2,746,078) and conformed to our criteria were defined as T2DM. The date of the NHIS medical checkup during that window was considered as baseline. Subjects younger than 40 years old and subjects with missing data were excluded. Because health information data existed from 2002, subjects with a history of fracture during a washout period from January 2002 to December 2008 were excluded. In addition, those who developed fracture or died within 1 year from the health screening day were excluded. The reason for excluding the 1-year lag is follows. First, if the period from the index date to the occurrence of hip fracture is too short, it is difficult to determine whether it is caused by weight change. Second, problem of reverse causation may also be raised. After all, a total of 2,086,187 participants were included in our study. The cohorts were followed from the day a patient received health screening, to the occurrence of any fracture, or the last follow-up day (December 31, 2018), whichever came first. The mean duration of follow-up was 6.56±1.92 years. Flow chart depicting the design of study cohorts are shown in Fig. 1. This study protocol was reviewed and approved by the Institutional Review Board of Soongsil University (SSU-202003-HR-201- 01), and was conducted in accordance with the tenets of the Helsinki Declaration. Because all data provided by the NHIDSS to researchers were anonymized, he need for informed consent was exempted.
Data collection
Information regarding health-related lifestyle (including income level, smoking, drinking, and physical activity), medical history, and clinical tests results are included in the questionnaire obtained during the mandatory health screening. In addition, anthropometric measurements including height, body weight, waist circumference and blood pressure are checked manually by medical staff on duty during the mandatory health screening. Comorbidities were largely defined using a combination of past history (ICD-10 code and self-reported) and use of medication history for the corresponding disease. Hypertension was defined when at least one claim per year under ICD10 codes I10 or I11 and at least one claim per year for the prescription of an antihypertensive agent or systolic/diastolic blood pressure ≥140/90 mm Hg was present. Blood samples for the measurement of serum glucose, total cholesterol, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol levels and triglycerides were collected after fasting overnight. BMI was calculated by dividing the weight (in kg) by the square of the height (in m2). We defined obesity as BMI ≥25 kg/m2. The participants were then classified into the following categories according to The Korean Society for the Study of Obesity recommendations: Underweight (BMI <18.5 kg/m2), normal (≥18.5 to <23 kg/m2), overweight (≥23 to <25 kg/m2), obese (≥25 to 30 kg/m2), and morbidly obese (≥30 kg/m2) [21,22].
Fracture identification depending on the skeletal lesions
According to the anatomical location of the fracture, the incidence rate (IR) and hazard ratio (HR) were analyzed by dividing it into four categories: vertebra fracture, hip fracture, limbs fracture, and total fracture. Participants having ICD-10 codes for vertebral fracture (S22.0, S22.1, S32.0, M48.4, and M48.5) and who have visited the hospital more than twice due to same codes were defined as having vertebral fracture. Likewise, those having ICD-10 codes and who have visited the hospital more than twice with fracture of upper arm (S42.0, S42.2, and S42.3), forearm (S52.5 and S52.6), or lower leg (S82.3, S82.5, and S82.6) were classified as having a limb fracture. Because almost all elderly with hip fracture requires either surgical treatment or supportive care via hospital admission, hip fracture was defined as having ICD-10 codes of hip fracture (S72.0 and S72.1) and history of hospitalization due to hip fracture. Finally, total fractures included fractures of vertebra, limbs, hip, and others not listed above (i.e., S02.X for skull fracture, S12. X for neck fracture, S62.X for hand, S92.X for foot) [23].
Statistical analysis
Baseline characteristics are presented as mean±standard deviation, median value (interquartile range) or number (%). Difference between groups in baseline demographic and clinical characteristics were compared using analysis of variance (ANOVA) for continuous variables and chi-square test for categorical variables. The incidence of fracture occurrence was calculated by dividing the number of incident cases by the total follow-up duration (person-years). We performed Cox proportional hazard regression, with normal BMI (≥18.5 to <23 kg/m2) as reference category, to evaluate the risk of fractures. In addition, Cox proportional- hazard regression for fracture of vertebral, hip, and limbs were conducted to assess whether the association between BMI and fracture risk differed depending on the affected or fracture lesions. The cox proportional-hazard model was adjusted for potential confounding variables known to predict risk of fracture, which included age, sex, smoke, drink, regular exercise, hypertension, dyslipidemia, chronic kidney disease (CKD), insulin use, diagnosis of T2DM for more than 5 years, use of three or more oral hypoglycemic agents, fasting glucose level, and height. The proportional hazard assumption was checked by visual inspection of stratified –log(–log) survival curves using tertile classification of continuous variables and by the test of Harrel and Lee on Schoenfeld residuals in univariate models [24]. The disease-free probability of primary outcomes related to the BMI was calculated using Kaplan– Meier curves, and log-rank test was performed to analyze the differences between the groups. HRs and 95% confidence intervals (CIs) for fractures were calculated using a Cox proportional hazards model for each category. In addition, the mortality rate as a competing risk was analyzed (Supplementary Table 1). For all statistical analysis, we used SAS version 9.3 (SAS Institute, Cary, NC, USA) with P<0.05 considered as significant.
RESULTS
Characteristics of each group related to the prevalence of obesity and T2DM
Table 1 shows the characteristics of the study population classified into the five BMI categories. The mean baseline age was 58.6 years, and 60.2% of the sample were men. The proportion of patients with underweight (BMI <18.5 kg/m2) was 1.4% (n=28,229) among a total of 2,086,187 subjects, while the proportion of subjects with obesity (BMI >25 kg/m2) was 48.8%. The percentages of subjects with obesity among ‘obese’ (≥25 to 30 kg/m2) and ‘morbidly obese’ (≥30 kg/m2) patients were 41.6% and 7.1%, respectively. In our data analysis, the lipid profile measurements (total cholesterol, triglyceride, low-density lipoprotein and high-density lipoprotein) worsened and the waist circumference, the prevalence of hypertension, dyslipidemia, and cardiovascular disease increased as BMI level increased. On the other hand, the proportion of current smoker, subjects with cancer, insulin use, use of three or more oral hypoglycemic agents, diabetic duration more than 5 years, and fasting glucose level increased as BMI level decreased. The proportion of male subjects was low only in the ‘morbidly obese’ category, and similar results were seen in the rest of the categories. The proportion of low income level and subjects with CKD increased among diabetic patients with low and high BMI. The proportion of subjects with regular exercise habits decreased among diabetic patients with low and high BMI. Baseline characteristics of subjects are shown in Table 1.
Associations of BMI with fracture risk
To analyze the impact of BMI on the incidence of fracture, we divided the population into five groups according to their BMI results. The HR was calculated using subjects with normal BMI range (18.5 to 23 kg/m2) as the reference group (Table 2). We also used the Cox proportional hazard model to analyze the effects of BMI and T2DM on fracture risk. Results of the three models are shown in Table 2. The first model is unadjusted, the second is adjusted for age, sex, smoke, drink, and regular exercise, and the third is adjusted for hypertension, dyslipidemia, CKD, insulin use, duration of DM for more than 5 years, use of three or more oral hypoglycemic agents, fasting glucose level, and height in addition to model 2. As for cox regression, total fractures were significantly higher in the underweight group (adjusted hazard ratio [aHR], 1.268; 95% CI, 1.228 to 1.309) and lower in the obese group (aHR, 0.891; 95% CI, 0.882 to 0.901) and the morbidly obese group (aHR, 0.873; 95% CI, 0.857 to 0.89). In terms of specific fractures, similar trends were observed for vertebra fractures. However, the risk of hip fracture was significantly increased in the ‘underweight’ group (aHR, 1.896; 95% CI, 1.78 to 2.021) and decreased in the ‘obese’ (aHR, 0.643; 95% CI, 0.624 to 0.663) and ‘morbidly obese’ group (aHR, 0.627; 95% CI, 0.591 to 0.665), compared to vertebra fractures or total fractures. Lastly, risk of limb fractures showed a linear relationship with BMI increase (Table 2). The Kaplan–Meier estimations for the incidence probability of the risk of fractures according to BMI category are shown in Fig. 2.
Subgroup analyses and the fracture risk according to BMI by age and sex
For subgroup analyses, the subjects were divided into groups by age (40–64 and ≥65 years) and sex. The IRs and HRs of fractures according to affected bone were investigated for each group, and is demonstrated in Table 3. Additional data on P value for interaction with CKD, ischemic heart disease (IHD), duration of DM ≥5 years, use of three or more oral hypoglycemic agents excluding age and sex are shown in the Supplementary Table 2.
Age and sex affected total fractures, vertebral fractures, hip fractures, and limb fractures; its P value for interaction was 0.0006, <0.0001, <0.0001, 0.0012 for age group (40–64, more than 65), and <0.0001, <0.0001, <0.0001, <0.0001 for sex group, respectively. The presence of IHD and use of three or more oral hypoglycemic agents had no effect on total fractures, vertebral fractures, hip fractures, or fracture limbs; its P value for interaction was 0.4555, 0.9369, 0.4574, 0.1224 for IHD group, and 0.109, 0.8252, 0.9454, 0.2813 for ‘use of three or more oral hypoglycemic agents’ group, respectively. The presence of CKD had an effect on vertebra fracture (P for interaction=0.0053) and hip fracture (P for interaction <0.0001) compared with the cases without CKD, but did not affect total fractures (P for interaction=0.1779) and limb fracture (P for interaction=0.5165). Duration of DM ≥5 years had effect on the incidence of total fractures (P for interaction=0.015), but had little effect on vertebra fracture (P for interaction=0.0708), hip fracture (P for interaction=0.1407), and limb fractures (P for interaction=0.0586).
Among the for subgroups, male aged <65, female aged <65, male aged ≥65, and female aged ≥65 (Fig. 3), fracture risk according to the affected bone, showed the highest in the underweight group (BMI <18.5 kg/m2), and the higher the BMI was, the lower the fracture risk was. Exceptionally, the risk of fractures tended to increase as BMI increased for vertebra fractures in the female aged <65 group. In subjects of female aged <65, the fully aHRs of total fractures (Model 3) were 1 in the normal group (≥18.5 to <23 kg/m2), 0.984 (95% CI, 0.824 to 1.174) in the underweight group (BMI <18.5 kg/m2), 1.079 (95% CI, 1.028 to 1.132) in the overweight group (≥23 to <25 kg/m2), 1.133 (95% CI, 1.084 to 1.183) in the obese group (≥25 to 30 kg/m2), and 1.232 (95% CI, 1.16 to 1.31) in the morbidly obese group (≥30 kg/m2). In terms of hip fractures, BMI lower than 18.5 kg/m2 significantly increased the risk of hip fracture, regardless of the affected bone. Among them, the HR (3.12; 95% CI, 2.673 to 3.642) was especially increased in the male aged <65 group.
DISCUSSION
To the best of our knowledge, this study is the most large-scale nationwide population study to show the association between BMI and fracture risk varied by affected site in T2DM cohort. Low BMI is a both well-established and important risk factor for osteoporosis and osteoporotic fractures in both men and women, whereas high BMI appears to act as a protective factor [25]. In this population-based cohort study of a total of 2,086,187 subjects >40 years of age with diagnosis of T2DM, we described that there is a reverse J-shaped association between BMI and fracture risk. In the reverse J shape curve, in terms of amplitude, hip fracture was the largest, followed by vertebra, all, and limb fractures. In this T2DM cohort, the HR for hip fractures was highest in the underweight group (HR, 1.896; 95% CI, 1.178 to 2.021).
De Laet et al. [25] conducted a meta-analysis study comprising nearly 60,000 men and women from 12 cohorts of both Asian and Western participants. The study showed that low BMI was a significant risk factor for all types of fractures in both Asian and Western populations. Also, the authors reported that there were only 83 diabetes patients with hip fractures in their study, of which only six individuals had low BMI (BMI <18.5) [25]. Thus, the authors concluded that the result was not significant. Nonetheless, because a sufficiently large number of subjects were included in our diabetic cohort study, the event number of hip fractures, vertebra fractures, and limbs fractures at low BMI (BMI <18.5) were large enough for statistical analysis.
In addition to considering low BMI as a risk factor for fractures, the effect of high BMI scores on fracture risk has also been discussed recently. Some reports in the previous literature have shown that obesity or morbid obesity has a protective effect on the risk of hip fracture [26,27]. Sogaard et al. [28] stated that the risk of hip fracture decreased with increasing BMI. Johansson et al. [26] also analyzed over 300,000 women from more than 25 countries and found that 87% of hip fractures occurred in those without obesity (defined as BMI ≥30 kg/m2). The same results were seen in this study; higher BMI was protective against fracture risk regardless of affected site. However, in terms of amplitude, hip fracture was the largest, followed by vertebra, all, and limb fracture. Some articles stated that obesity acts as a protective effect on fractures due to consequent higher BMD and reduced impact of falls as a result of increased soft-tissue padding [27-29].
Some interesting results were seen in the subgroup analyses. According to the subgroup analysis by age, in the group younger than 65 years old, there was a higher HR for fracture risk among subjects who were underweight (BMI <18.5), compared with the group older than 65 years old (Supplementary Table 2). This tendency was evident in hip fracture, followed by vertebra fracture, then total fracture. This tendency was least evident in limbs fracture. Limited to hip fractures, it is known that the risk of hip fractures in the younger age group is higher than in the older age group [10,19,30]. Our study showed that although the amplitude was relatively smaller, total fractures and vertebral fractures showed similar trends to hip fractures. Some authors stated that the difference of the risk of hip fractures between young and old age groups comes from age- or menopause-related changes which overshadow the effect of diabetes itself on fracture risk [31]. We presume that the reason may be similar not only for hip fractures, but also for vertebra fractures and total fractures.
A subgroup analysis by sex showed that the risk of fractures in underweight (BMI <18.5) subjects was higher in men than in women (Supplementary Table 2). This trend appeared in the order of hip fractures, vertebra fractures, and total fractures, and the least in limb fractures. However, in the case of vertebral and total fractures as well as hip fractures, the effect on fracture risk due to low BMI was significantly higher in men than in women, which is a contrary result to the previous literature [10]. One possible reason for the disparity in fracture risk by sex is that men have larger bones and less adipose tissue than women [32]. Since the weight of bone remains relatively constant compared to muscle or fat tissue, the effect of low BMI on fracture risk may be more significant in males. Additionally, previous reports have shown that underweight men lose more muscle mass than women, which increases mechanical stress on the bones [33]. For instance, male subjects with low BMI would be more prone to have sarcopenia, which is related to physical disability and consequent incident falls [34]. Conversely, the increase in fracture risk was greater in females than in males when overweight, obese or morbidly obese. Similarly, since bone occupies a low proportion of body weight in females, the proportion of adipose tissue and muscle is relatively high, so higher proportion of adipose tissue in female can be estimated that the increase in fracture risk in overweight, obese or morbidly obese is greater in females than in males. One meta-analysis study demonstrated that abdominal obesity itself was associated with an increased risk of hip fracture [35].
In addition, among the four age-sex categories (male ≥65 years old, female ≥65 years old, male <65 years old, and female <65 years old), the HR for fracture risk was highest in the category of male <65 years old, followed by female <65 years old (Fig. 3). Judging from these results, among age and sex, the effect of BMI on fracture risk has greater interaction with age than with sex.
Our study has several notable strengths. It has a large sample size of >2,000,000 subjects and a long follow-up period of >6 years. Various subgroup analyses were possible using this large group of data, which provided interesting conclusions. The analyses were performed after adjusting for substantial confounding variables, including age, sex, smoke, drink, regular exercise, hypertension, dyslipidemia, CKD, insulin use, duration of DM for more than 5 years, use of three or more oral hypoglycemic agents, fasting glucose level, and height, which may cause BMI changes. Because the Korean population is a single ethnic society, it was possible to involve a homogeneous group in a nationwide study.
Despite these advantages, our study also has limitations. First, because data was extracted solely from the NHIS for analysis, not all health behaviors and other factors were included in the analysis, such as eating habits, body composition/muscle mass, muscle function/sarcopenia, and history of falls. In addition, osteopenia, osteoporosis, and bone mineral density are all important risk factors of fracture, but were not included in the analysis. Second, there may be reverse causation in our results because of the retrospective cohort design. However, washout period was considered to address this issue when study outcomes were assessed. Thus, subjects with a history of fracture during the washout period from January 2002 to December 2008 were excluded. Third, because the analytic sample was limited to Korean individuals, additional studies in other ethnic groups are required to generalize our results. Fourth, by using only the NHIS claim database to search for fractures, without reviewing medical and radiologic records, traumatic fractures may have been included in our analyses. The concept of major osteoporotic fracture was not used and that fractures were divided according to anatomical location. Major osteoporotic fracture refers to hip fracture, vertebra fracture, and forearm and humerus fracture, and represents osteoporotic fracture. Although limbs fractures include fractures of the humerus and forearm, limbs fracture is not often analyzed as a clinical outcome. The fact that limbs fractures in this article show less prominent results compared to other sites can also be presumed to have such an effect. Fifth, there may be a misclassification bias. In concern of exclusion of insulindependent T1DM, we did not use the ICD-10 code E10 representing T1DM as an exclusion criterion. However, there is a possibility that subjects with T1DM still may have been enrolled in our study. Fortunately, East Asian countries, including the Republic of Korea, have the lowest incidence of T1DM in the world [36]. The prevalence of T1DM in the entire population of the Republic of Korea was reported to be from 0.017% to 0.021% [37]. The majority of diabetes is T2DM, and the proportion increases with age. We excluded subjects younger than 40 years in this study. Consequently, the possibility of inclusion of T1DM is minimal and was regarded insignificant.
In conclusion, we further demonstrated that fracture risk differed depending on the skeletal sites. For the hip, the risk of fracture was most pronounced in the underweight group, with a slight increase in the risk in obese patients as well. The fracture risk of vertebra also showed a similar trend to that of the hip, but the amplitude was relatively smaller. For vertebra fractures, the increase in fracture risk in the underweight group was the second highest following hip fractures, and the decrease in fracture risk in obese or morbidly obese subjects was also the second highest following hip. For limbs, fracture risk was least affected by BMI, showing a linear relationship.
SUPPLEMENTARY MATERIALS
Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2022.0001.
Notes
CONFLICTS OF INTEREST
No potential conflict of interest relevant to this article was reported.
AUTHOR CONTRIBUTIONS
Conception or design: S.W.L., K.H., H.S.K.
Acquisition, analysis, or interpretation of data: S.W.L., K.H., H.S.K.
Drafting the work or revising: S.W.L., H.S.K.
Final approval of the manuscript: S.W.L., K.H., H.S.K.
FUNDING
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2021R1A2C2013890).
Acknowledgements
None