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

Association between insulin resistance and lung function trajectory over 4 years in South Korea: community-based prospective cohort

Abstract

Background

Hyperglycemic conditions are associated with respiratory dysfunction. Although several studies have reported that insulin resistance (IR) is related to decreased lung function, the association between IR and change in lung function has been rarely studied. This study aimed to investigate the potential association of IR on annual change in lung function using a community-based prospective cohort in Korea.

Methods

We selected 4827 Korean participants whose serial lung functions were assessed over 4 years using 1:3 propensity score matching. Exposure was baseline IR estimated with homeostatic model assessment (HOMA-IR), and outcomes were annual changes in lung function determined by calculating the regression coefficient using least-square linear regression analysis.

Results

In the multivariate linear regression, per one unit increased log transformed HOMA-IR was associated with decline in FEV1%-predicted (β: − 0.23, 95% CI: − 0.36 to − 0.11) and FVC %-predicted (β: − 0.20, 95% CI: − 0.33 to − 0.08), respectively. In the generalized additive model plot, HOMA-IR showed a negative linear association with annual changes in FEV1%-predicted and FVC %-predicted. The suggested threshold of HOMA-IR for decline in lung function was 1.0 unit for annual change in FEV1%-predicted and 2.2 unit for annual change in FVC %-predicted. Age showed statistically significant effect modification on the relationship between HOMA-IR and annual change in FEV1%-predicted. Increased HOMA-IR was associated with the decreased annual change in FEV1%-predicted, particularly in older people.

Conclusions

In South Korea, increased HOMA-IR was associated with decline in lung function. Since IR was related to decline in FEV1%-predicted, particularly in older people, tailored approaches are needed in these populations. The potential pulmonary hazard of IR needs to be confirmed in future studies.

Peer Review reports

Background

Diabetes mellitus is a prevalent chronic disease in modern society, with many experiencing diabetes and its complications [1, 2]. Lung function decline has been considered a major complication of diabetes [3,4,5]. Several studies have found that diabetes and hyperglycemia are associated with the development of various pulmonary diseases, including chronic obstructive pulmonary disease, asthma, and interstitial lung disease [6,7,8]. Therefore, good control of diabetes may prevent future lung diseases.

Increased insulin resistance (IR) has been associated with the future development of obesity, diabetes [9], cardiovascular diseases [10], neurologic impairment [11], and kidney dysfunction [12]. Although several studies have reported the potential association between IR and decreased lung function, none of these studies have evaluated the longitudinal effect of IR on future lung function changes. The aging process can cause diverse physiological, immunological and structural changes in the respiratory system [13, 14]. Thus, analyzing consecutive measurements of lung function may help to interpret the effects of aging on IR and lung function changes. In addition, it can clarify a causal association between IR and lung function.

Since diabetes is associated with future decline in lung function and IR is a fundamental pre-condition of diabetes, we hypothesized that increased IR may be associated with future decline in lung function. Therefore, this study aimed to investigate the association between IR and lung function change using data from a community-based prospective Ansan-Ansung cohort in Korea.

Methods

Participants

The Ansan-Ansung cohort comprise participants aged 40–69 years who lived in Ansan (urban) and Ansung (rural) areas, and this cohort was analyzed to determine the factors affecting the incidence of chronic diseases. Baseline measurements were performed between May 2001 and February 2003, and two subsequent lung function tests were performed biennially thereafter. More detailed information about the Ansan-Ansung cohort can be obtained in previous reports [15, 16]. Of the 10,030 participants, we excluded 983 missing baseline measurements. Among the 9047 subjects, 4827 subjects with additional spirometry were selected using 1:3 propensity score matching with 1609 subjects without additional spirometry. The differences in clinical characteristics of participants with or without additional spirometry were described in Additional file 3: Table S1. Therefore, 4827 matched participants were finally included in this analysis (Fig. 1).

Fig. 1
figure 1

Flow chart of the study subject selection

Exposure: IR

IR was assessed using the homeostatic model assessment for IR (HOMA-IR). HOMA-IR was calculated using the following equation: fasting insulin (μIU/mL) × fasting glucose (mg/dL)/405 [17].

Outcomes: change in lung function

The main outcomes were annual changes in lung function. Lung function tests were performed using a spirometer (VMAX2130; Sensormedics Corporation, Yorba, CA, USA). Pre-bronchodilator values were measured three times by trained technicians, and the best scores were recorded. Predicted values of functional vital capacity (FVC %-predicted) and functional expiratory volume in 1 s (FEV1%-predicted) were calculated using the Korean formula [18]. The FEV1/FVC ratio was expressed as the percentage of raw value of FEV1 divided by the raw value of FVC, according to previous literature [19, 20]. Annual changes in FEV1%-predicted, FVC %-predicted, and FEV1/FVC ratio was determined by calculating the regression coefficient using least-square linear regression analysis, with FEV1%-predicted, FVC %-predicted, and FEV1/FVC ratio as functions of time in years. It was applied to all values of FEV1%-predicted, FVC %-predicted, and FEV1/FVC ratio that was obtained during the follow-up period (Additional file 2: Figure S1).

Measurements and other definitions

Trained investigators interviewed participants regarding their socioeconomic status and lifestyle habits using Additional file 1: KoGES Baseline Core Questionnaire. High income was defined as the highest quintile of monthly household income (≥ 3 million won a month). Active physical activity was defined as 60 min/day of moderate-intensity activity or 30 min/day of vigorous-intensity activity [21]. Blood pressure (BP) was defined as the average BP on both arms using a standard mercury sphygmomanometer (Baumanometer-Standby; W. A. Baum Co., Inc., Copiague, NY, USA). High BP was defined as systolic BP > 140 mmHg or diastolic BP > 90 mmHg. Body mass index (BMI) was defined by dividing the weight by the square of the height (kg/m2). Waist circumference was measured at the narrowest point between the lower rib and the iliac crest (measured to the nearest 0.1 cm). All blood samples were examined after fasting for at least 8 h. Hemoglobin and white blood cell (WBC) count were analyzed using enzymatic methods with ADVIA 120 (Bayer Diagnostics, Tarrytown, NY, USA). Fasting glucose, hemoglobin A1c (HbA1c), triglyceride, high density lipoprotein (HDL) cholesterol, C-reactive protein (CRP), and serum creatinine levels were measured using ADVIA 1650 (Siemens, Tarrytown, NY, USA). Estimated glomerular filtration rate (eGFR) was defined using the Chronic Kidney Disease Epidemiology Collaboration equation [22]. Annual averaged values were obtained for BMI, waist circumference, systolic and diastolic BP, triglyceride, HDL cholesterol, eGFR, WBC count, fasting glucose, HbA1c, hemoglobin, and CRP level.

Statistical analyses

We compared individuals without additional spirometry or with missing values to the rest of participants. Variables that showed significant differences were selected for matching. As a result, age, sex, high income, college graduate, smoking status, active physical activity, baseline values of BMI, systolic BP, diastolic BP, WBC, CRP, and FEV1%-predicted were used in the propensity score matching. The normality of the distribution of continuous variables was assessed using histogram and Q–Q plots. Mean ± standard deviation (SD) was used to express normally distributed continuous variables, whereas median with interquartile range (IQR) was used for non-normally distributed continuous variables. The P-trend was analyzed using linear regression for the normally distributed continuous variables, the Jonckheere–Terpstra test for the non-normally distributed continuous variables, and the Cochran–Armitage test for categorical variables. We further analyzed the difference between individuals with and without additional spirometry using student’s t-test for a normally distributed continuous variable, Wilcoxon rank-sum test for non-normally distributed continuous variable, and Pearson’s Chi-squared test for categorical variables. The correlation between age and HOMA-IR was performed using Spearman’s rank correlation analysis.

The simple associations between HOMA-IR and spirometric parameters were expressed using the locally weighted scatterplot smoothing (LOWESS) function. We performed the linear regression analysis between log-transformed HOMA-IR and annual changes in FEV1%-predicted and FVC %-predicted. Odds ratio (OR) and 95% confidence interval (CI) for the lowest quartile of annual changes in FEV1%-predicted and FVC %-predicted were analyzed using logistic regression analysis.

In multivariate analyses, covariates were chosen based on clinical and statistical relevance. The non-linear association between HOMA-IR and annual changes in FEV1%-predicted and FVC %-predicted was analyzed using the multivariate generalized additive model (GAM) for Gaussian distributions with the “mgcv” package, and calculation of Akaike information criterion (AIC) was used for the model fitting. In AIC, lower scores within the data set indicate a better model fit [23]. The threshold point of HOMA-IR was chosen based on the best fit determined by AIC among the models.

In the subgroup analysis, all control variables were used, and participants were divided by the median for continuous variable. The modulatory effect was confirmed by adding the interaction term of variables defining the subgroup to the multivariate linear regression analysis. To visualize potential interaction of age on the association between IR and annual change in lung function, the LOWESS regression with interaction was implemented using the “predict3d” package. A P-value < 0.05 was considered statistically significant. All statistical analyses were performed using R version 3.6.2 (R Core Team 2019; R Foundation for Statistical Computing, Vienna, Austria).

Results

The mean age of the 4827 participants was 52.6 years, 45.3% were men, 48.0% had active physical activity, and 41.7% were current or past smokers. The mean ± SD of annual averaged BMI, systolic and diastolic BP, waist circumference, HDL cholesterol, eGFR, hemoglobin, and WBC count were 24.4 ± 3.0 kg/m2, 119.6 ± 15.6 mmHg, 79.0 ± 9.4 mmHg, 83.2 ± 8.6 cm, 44.9 ± 8.8 mL/dL, 82.9 ± 10.9 mL/min/1.73m2, 13.6 ± 1.4 g/dL, and 6.5 ± 1.5 × 103/μL, respectively. Moreover, the median (IQR) of annual averaged fasting glucose, HbA1c, triglyceride, and CRP level were 88 (83–94) mg/dL, 5.6 (5.3–5.8) %, 130 (97–178) mg/dL, and 0.12 (0.07–0.22) mg/dL, respectively. At baseline, the mean ± SD values of FEV1%-predicted, FVC %-predicted, and FEV1/FVC ratio were 94.9 ± 14.1% predicted, 95.4 ± 13.2% predicted, and 79.5 ± 7.9%, respectively. The mean of FEV1%-predicted, FVC %-predicted, and FEV1/FVC ratio decreased in subsequent lung function measurement (Fig. 2). The mean ± SD values of annual changes in FEV1%-predicted, FVC %-predicted, and FEV1/FVC ratio were − 0.42 ± 2.62%-predicted/year, − 0.49 ± 2.72% predicted/year, and − 0.35 ± 1.56%/year, respectively.

Fig. 2
figure 2

Mean lung function trajectories over 4 years. Abbreviations: FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity

The clinical characteristics of the study population according to HOMA-IR quartile are shown in Table 1. Although age was not associated with HOMA-IR quartile, the increased HOMA-IR quartile was associated with a decreased proportion of men. The correlation analysis between age and HOMA-IR also was not statistically significant (Additional file 2: Fig. S2). Although income and educational status were not associated with HOMA-IR quartile, the increased HOMA-IR quartile was associated with a decreased rate of active physical activity and current or past smoking. With the increase of HOMA-IR quartile, BMI, systolic and diastolic BP, waist circumference, WBC count, fasting glucose, HbA1c, triglyceride, hemoglobin, and CRP level increased, whereas HDL cholesterol decreased. With the increase of HOMA-IR quartile, baseline FVC %-predicted decreased, but FEV1/FVC increased. Baseline FEV1%-predicted was not associated with HOMA-IR.

Table 1 Clinical characteristics of the study population according to HOMA-IR quartile

We explored the association between HOMA-IR and annual changes in lung function (Fig. 3). As HOMA-IR increased, annual changes in FEV1%-predicted and FVC %-predicted change decreased, whereas annual change in FEV1/FVC ratio change increased. In multivariate linear regression analysis, per one unit increased log transformed HOMA-IR was associated with decline in FEV1%-predicted (β: − 0.23, 95% CI: − 0.36 to − 0.11, p < 0.001) and FVC %-predicted (β: − 0.20, 95% CI: − 0.33 to − 0.08, p = 0.001) (Table 2). In the categorical analysis, the odds for the lowest quartile of annual changes in FEV1%-predicted (Q1: <  − 1.73%-predicted/year) and FVC %-predicted (Q1: − 1.90%-predicted/year) showed increasing trend (P-trend < 0.001) as quartiles of HOMA-IR were increased. The OR of HOMA-IR for lowest quartile of annual changes in FEV1%-predicted and FVC %-predicted were significantly higher in the upper quartiles compared to the lower quartile. In the multivariate GAM plot analysis, we identified that the association between HOMA-IR and annual changes in FEV1%-predicted and FVC %-predicted were negative linear, and the lowest AICs were found in HOMA-IR 1.0 unit for annual change in FEV1%-predicted and 2.2 unit for annual change in FVC %-predicted (Fig. 4).

Fig. 3
figure 3

The association between HOMA-IR and annual changes in lung function. The red line represented the LOWESS regression curve. Abbreviations: HOMA-IR, homeostatic model assessment-insulin resistance; LOWESS, locally weighted scatter-plot smoothing

Table 2 The effects of HOMA-IR on the lung function change
Fig. 4
figure 4

GAM plot between HOMA-IR and annual changes in lung function. The dashed lines indicate 95% CIs for values of smoothed annual changes in FEV1%-predicted and FVC %-predicted, using the multivariate GAM analysis after adjusting for age, sex, high income, college graduate, active physical activity, smoking status, and annual averaged values of BMI, waist circumference, systolic and diastolic BP, triglyceride, HDL cholesterol, eGFR, WBC count, hemoglobin, and CRP level. Abbreviations: GAM, general additive model; HOMA-IR, homeostatic model assessment-insulin resistance; CI, confidence interval; FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; BMI, body mass index; BP, blood pressure; HDL, high density lipoprotein; eGFR, estimated glomerular filtration rate; WBC, white blood cell; CRP, C-reactive protein

There were no subgroups showing the effect modification between HOMA-IR and annual changes in FVC % predicted (Additional file 3: Table S3). However, participant age significantly modified the effect of HOMA-IR on annual change in FEV1%-predicted (p for interaction = 0.003). Although HOMA-IR was not associated with annual change in FEV1%-predicted in people aged < 50 years, increased HOMA-IR was significantly associated with the decreased annual change in FEV1%-predicted in individuals aged ≥ 50 years (Table 3). We also identified similar effect modification of age on the association between HOMA-IR and annual change in FEV1%-predicted in the LOWESS regression plot with interaction (Fig. 5).

Table 3 Subgroup analysis for the effects of HOMA-IR on annual change in lung function
Fig. 5
figure 5

LOWESS regression plot with interaction for the effect modification of age on the association between HOMA-IR and annual change in FEV1%-predicted. Solid lines represented the univariate LOWESS regression curves modulated by age. The darkest blue, and lightest blue lines represented lines of mean minus one SD (44 years) and plus one SD (62 years). Abbreviations: LOWESS, locally weighted scatter-plot smoothing; HOMA-IR, homeostatic model assessment-insulin resistance; FEV1, forced expiratory volume in 1 s; SD, standard deviation

Discussion

Our study revealed that increased IR associated with decline in lung function using community-based prospective cohort in Korea. Recently, decreased lung function has been suggested as a complication of IR. In their large-scale cross-sectional study, Lawlor et al. found that HOMA-IR decreased by 5% and 8% for a 1 SD increase in log FEV1 and FVC, respectively [24]. Sagun et al. also showed that decreased FEV1%-predicted was correlated with increased insulin resistance (HOMA-IR ≥ 2.5) in an age-adjusted analysis (OR 0.97, 95% CI 0.94–1.00, P = 0.028) [25]. However, most previous studies were cross-sectional, and a causal association between insulin resistance and decreased lung function needed to be evaluated in a prospective cohort. Therefore, we conducted the current study using a prospective community-based cohort and found that increased insulin resistance was associated with a decline in lung function, particularly in the elderly population.

In the current study, increased HOMA-IR was independently associated with decline in FEV1%-predicted and FVC %-predicted. Cigarette smoking is a well-established risk factor for lung function deterioration. In the large scale randomized clinical trial, annual change in FEV1%-predicted was approximately − 1.4%-predicted/year in the continuing smoker group [26]. Compared to this result, IR can be considered a risk factor for rapid decline in lung function. In LOWESS regression analysis, increased HOMA-IR presented negative association with annual changes in FEV1%-predicted and FVC %-predicted, but not with FEV1/FVC ratio. In the natural history of chronic obstructive pulmonary disease, the rate of FEV1 decline is slow at high FEV1 and faster as FEV1 decreases [27]. Due to the nature of the community-based cohort, the majority of participants had a high baseline FEV1. These factors could explain why HOMA-IR was not associated with annual change in FEV1/FVC ratio.

We also found negative linear associations between HOMA-IR and annual changes in FEV1%-predicted and FVC %-predicted using the multivariate GAM plot, and the suggested threshold of HOMA-IR for decline in lung function was 1.0 unit for annual change in FEV1%-predicted and 2.2 unit for annual change in FVC %-predicted. As a surrogate marker of IR, several studies have shown that the cutoff value of HOMA-IR for the risk of developing type 2 diabetes mellitus and metabolic syndrome in the Korean population is approximately 2.5 [28,29,30]. Our findings showed lower threshold than those of previous studies, indicating that decline in lung function may be more sensitive to the systemic effects of IR. Moreover, these results suggest that early treatment of IR can be an effective intervention to prevent decline in lung function.

There were two possible explanations for the potential pulmonary hazard of IR. First, chronic low-grade inflammation in adipose tissue contributes to the development of IR [31, 32]. Levels of cytokines, such as tumor necrosis factor-alpha (TNF-α) and interleukin-6, are increased with IR and cause lung fibrosis and inflammation [33, 34]. Second, respiratory muscle weakness is related to IR. In high IR status, TNF-α and interferon-gamma cause muscle wasting through the degeneration of myotubes and modulation of myogenesis [35, 36].

In subgroup analysis, age showed a statistically significant effect modification on the association between HOMA-IR and annual change in FEV1%-predicted. Although HOMA-IR was not associated with annual change in FEV1%-predicted in the younger group, increased HOMA-IR was significantly related to decline in FEV1%-predicted in the older group.

IR is associated with chronic low-grade inflammation [31]. The risk of chronic low-grade inflammation increased with the aging process [37, 38]. Lung function can decline with increased inflammation [39, 40]. Therefore, we assumed that increased chronic low-grade inflammation during the aging process may potentiate the metabolic hazard of IR in lung function deterioration. In this regard, early intervention to reduce IR in the elderly population may have a beneficial effect on future lung function preservation.

Several limitations should be considered when interpreting our findings. First, two or fewer additional lung function measurements were performed in this study. Therefore, our study period was not sufficiently long. Further studies performing more spirometric tests will be needed to clarify the long-term effect of IR on decline in lung function. Second, our study was conducted in single ethnicity and selected areas. Therefore, the results should be elucidated cautiously in accordance with racial and regional differences.

Despite these limitations, our study has several strengths. First, to the best of our knowledge, this is the first study to investigate the association between IR and lung function using a large-scale prospective cohort based on the general population. Second, we calculated the threshold of HOMA-IR for decline in lung function using a nonlinear analytic method. Finally, the interaction between age and IR to annual change in FEV1%-predicted was investigated and presented using visualization.

Conclusion

In South Korea, increased IR was independently associated with decline in FEV1%-predicted and FVC %-predicted. We also found that age modified the relationship between IR and FEV1%-predicted, particularly valid in older people. The potential pulmonary hazard of IR needs to be confirmed in future interventional studies.

Availability of data and materials

The data of our study is fully available when manuscript is accepted for publication.

Abbreviations

HOMA-IR:

Homeostatic model assessment-insulin resistance

FEV1:

Forced expiratory volume in 1 s

FVC:

Functional vital capacity

OR:

Odds ratio

CI:

Confidence interval

BP:

Blood pressure

BMI:

Body mass index

WBC:

White blood cell

HbA1c:

Hemoglobin A1c

HDL:

High density lipoprotein

CRP:

C-reactive protein

eGFR:

Estimated glomerular filtration rate

SD:

Standard deviation

IQR:

Interquartile range

LOWESS:

Locally weighted scatter-plot smoothing

GAM:

Generalized additive model

AIC:

Akaike’s information criterion

TNF-α:

Tumor necrosis factor-alpha

References

  1. Zheng Y, Ley SH, Hu FB. Global aetiology and epidemiology of type 2 diabetes mellitus and its complications. Nat Rev Endocrinol. 2018;14(2):88.

    Article  Google Scholar 

  2. van Susan D, Beulens JW, Yvonne T. van der S, Grobbee DE, Nealb B: The global burden of diabetes and its complications: an emerging pandemic. Eur J Cardiovasc Prev Rehabil. 2010;17(1_suppl):s3–8.

    Article  Google Scholar 

  3. Klein O, Krishnan J, Glick S, Smith LJ. Systematic review of the association between lung function and Type 2 diabetes mellitus. Diabet Med. 2010;27(9):977–87.

    Article  CAS  Google Scholar 

  4. Yeh H-C, Punjabi NM, Wang N-Y, Pankow JS, Duncan BB, Cox CE, Selvin E, Brancati FL. Cross-sectional and prospective study of lung function in adults with type 2 diabetes: the Atherosclerosis Risk in Communities (ARIC) study. Diabetes Care. 2008;31(4):741–6.

    Article  Google Scholar 

  5. Khateeb J, Fuchs E, Khamaisi M. Diabetes and lung disease: a neglected relationship. Rev Diabetic Stud RDS. 2019;15:1.

    Article  Google Scholar 

  6. Mamillapalli C, Tentu R, Jain NK, Bhandari R. COPD and type 2 diabetes. Curr Respir Med Rev. 2019;15(2):112–9.

    Article  CAS  Google Scholar 

  7. Enomoto T, Usuki J, Azuma A, Nakagawa T, Kudoh S. Diabetes mellitus may increase risk for idiopathic pulmonary fibrosis. Chest. 2003;123(6):2007–11.

    Article  Google Scholar 

  8. Perez MK, Piedimonte G. Metabolic asthma: is there a link between obesity, diabetes, and asthma? Immunol Allergy Clin. 2014;34(4):777–84.

    Article  Google Scholar 

  9. Gallagher EJ, LeRoith D, Karnieli E. The metabolic syndrome—from insulin resistance to obesity and diabetes. Endocrinol Metab Clin North Am. 2008;37(3):559–79.

    Article  CAS  Google Scholar 

  10. Ginsberg HN. Insulin resistance and cardiovascular disease. J Clin Investig. 2000;106(4):453–8.

    Article  CAS  Google Scholar 

  11. Kim B, Feldman EL. Insulin resistance in the nervous system. Trends Endocrinol Metab. 2012;23(3):133–41.

    Article  CAS  Google Scholar 

  12. De Cosmo S, Menzaghi C, Prudente S, Trischitta V. Role of insulin resistance in kidney dysfunction: insights into the mechanism and epidemiological evidence. Nephrol Dial Transplant. 2013;28(1):29–36.

    Article  Google Scholar 

  13. Sharma G, Goodwin J. Effect of aging on respiratory system physiology and immunology. Clin Interv Aging. 2006;1(3):253.

    Article  CAS  Google Scholar 

  14. Skloot GS. The effects of aging on lung structure and function. Clin Geriatr Med. 2017;33(4):447–57.

    Article  Google Scholar 

  15. Kim Y, Han B-G, Group K. Cohort profile: the Korean genome and epidemiology study (KoGES) consortium. Int J Epidemiol. 2017;46(2):e20–e20.

    Article  Google Scholar 

  16. Woo A, Lee SW, Koh HY, Kim MA, Han MY, Yon DK. Incidence of cancer after asthma development: 2 independent population-based cohort studies. J Allergy Clin Immunol. 2021;147(1):135–43.

    Article  CAS  Google Scholar 

  17. Wallace TM, Levy JC, Matthews DR. Use and abuse of HOMA modeling. Diabetes Care. 2004;27(6):1487–95.

    Article  Google Scholar 

  18. Choi JK, Paek D, Lee JO. Normal predictive values of spirometry in Korean population. Tuberc Respir Dis. 2005;58(3):230–42.

    Article  Google Scholar 

  19. Global Strategy for Asthma Prevention and Treatment. 2020 Update, https://ginasthma.org/wp-content/uploads/2020/06/GINA-2020-report_20_06_04-1-wms.pdf.

  20. Singh D, Agusti A, Anzueto A, Barnes PJ, Bourbeau J, Celli BR, Criner GJ, Frith P, Halpin DM, Han M: Global strategy for the diagnosis, management, and prevention of chronic obstructive lung disease: the GOLD science committee report 2019. European Respiratory Journal 2019, 53(5).

  21. Yang YJ. An overview of current physical activity recommendations in primary care. Korean J Fam Med. 2019;40(3):135–42.

    Article  Google Scholar 

  22. Levey AS, Stevens LA, Schmid CH, Zhang Y, Castro AF III, Feldman HI, Kusek JW, Eggers P, Van Lente F, Greene T. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.

    Article  Google Scholar 

  23. Anderson D, Burnham K. Model selection and multi-model inference. Second NY: Springer-Verlag. 2020;2004(63):10.

    Google Scholar 

  24. Lawlor D, Ebrahim S, Smith GD. Associations of measures of lung function with insulin resistance and type 2 diabetes: findings from the British Women’s Heart and Health Study. Diabetologia. 2004;47(2):195–203.

    Article  CAS  Google Scholar 

  25. Sagun G, Gedik C, Ekiz E, Karagoz E, Takir M, Oguz A. The relation between insulin resistance and lung function: a cross sectional study. BMC Pulm Med. 2015;15(1):139.

    Article  Google Scholar 

  26. Anthonisen NR, Connett JE, Murray RP. Smoking and lung function of Lung Health Study participants after 11 years. Am J Respir Crit Care Med. 2002;166(5):675–9.

    Article  Google Scholar 

  27. Vestbo J, Lange P. Natural history of COPD: Focusing on change in FEV1. Respirology (Carlton, Vic). 2016;21(1):34–43.

    Article  Google Scholar 

  28. Kim B, Choi HY, Kim W, Ahn C, Lee J, Kim JG, Kim J, Shin H, Yu JM, Moon S. The cut-off values of surrogate measures for insulin resistance in the Korean population according to the Korean Genome and Epidemiology Study (KOGES). PLoS ONE. 2018;13(11):e0206994.

    Article  Google Scholar 

  29. Yun K-J, Han K, Kim MK, Park Y-M, Baek K-H, Song K-H, Kwon H-S. Insulin resistance distribution and cut-off value in Koreans from the 2008–2010 Korean National Health and Nutrition Examination Survey. PLoS ONE. 2016;11(4):e0154593.

    Article  Google Scholar 

  30. Baek JH, Kim H, Kim KY, Jung J. Insulin resistance and the risk of diabetes and dysglycemia in Korean general adult population. Diabetes Metab J. 2018;42(4):296–307.

    Article  Google Scholar 

  31. De Luca C, Olefsky JM. Inflammation and insulin resistance. FEBS Lett. 2008;582(1):97–105.

    Article  Google Scholar 

  32. Zatterale F, Longo M, Naderi J, Raciti GA, Desiderio A, Miele C, Beguinot F. Chronic adipose tissue inflammation linking obesity to insulin resistance and type 2 diabetes. Front Physiol. 2020;10:1607–1607.

    Article  Google Scholar 

  33. Lundblad LK, Thompson-Figueroa J, Leclair T, Sullivan MJ, Poynter ME, Irvin CG, Bates JH. Tumor necrosis factor–α overexpression in lung disease: a single cause behind a complex phenotype. Am J Respir Crit Care Med. 2005;171(12):1363–70.

    Article  Google Scholar 

  34. Rincon M, Irvin CG. Role of IL-6 in asthma and other inflammatory pulmonary diseases. Int J Biol Sci. 2012;8(9):1281.

    Article  CAS  Google Scholar 

  35. Laghi F, Tobin MJ. Disorders of the respiratory muscles. Am J Respir Crit Care Med. 2003;168(1):10–48.

    Article  Google Scholar 

  36. Londhe P, Davie JK. Gamma interferon modulates myogenesis through the major histocompatibility complex class II transactivator, CIITA. Mol Cell Biol. 2011;31(14):2854–66.

    Article  CAS  Google Scholar 

  37. Fagiolo U, Cossarizza A, Scala E, Fanales-Belasio E, Ortolani C, Cozzi E, Monti D, Franceschi C, Paganelli R. Increased cytokine production in mononuclear cells of healthy elderly people. Eur J Immunol. 1993;23(9):2375–8.

    Article  CAS  Google Scholar 

  38. Krabbe KS, Pedersen M, Bruunsgaard H. Inflammatory mediators in the elderly. Exp Gerontol. 2004;39(5):687–99.

    Article  CAS  Google Scholar 

  39. Jiang R, Burke GL, Enright PL, Newman AB, Margolis HG, Cushman M, Tracy RP, Wang Y, Kronmal RA, Barr RG. Inflammatory markers and longitudinal lung function decline in the elderly. Am J Epidemiol. 2008;168(6):602–10.

    Article  Google Scholar 

  40. Hancox RJ, Poulton R, Greene JM, Filsell S, McLachlan CR, Rasmussen F, Taylor DR, Williams MJ, Williamson A, Sears MR. Systemic inflammation and lung function in young adults. Thorax. 2007;62(12):1064–8.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to express our gratitude to the NRF for funding to conduct Ansan-Ansung cohort (Grant No: NRF-2019R1G1A1099753). We also thank all participants and staffs of the Ansan-Ansung cohort.

Funding

This study was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MIST) (Grant No: NRF-2019R1G1A1099753). The funder had no role in the design and conduct of the study, data curation and analysis, or drafting of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

S.W.L. conceptualization. S.H.K. data curation/formal analysis. H.S.K. investigation/ methodology. S.H.K. writing – original draft. H.K.M. and H.S.K. writing – review and editing. S.W.L. supervision/validation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Sung Woo Lee.

Ethics declarations

Ethics approval and consent to participate

Using the Ansan-Ansung cohort, this study was conducted according to the Declaration of Helsinki. The institutional review boards of the Nowon Eulji Medical Center, Eulji University approved the protocol and ethics of the study (Approval Number: 2019-06-014). Participants provided written informed consent, and all data were completely anonymized prior to access.

Consent for publication

Not applicable.

Competing interests

All authors have no conflict of interest to declare.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1

. KoGES baseline Core Questionnaire.

Additional file 2

: Figure S1. Definitions of outcome and exposure variables. Abbreviations: HOMA-IR, homeostatic model assessment-insulin resistance. Figure S2. Correlation between age and HOMA-IR. Abbreviations: HOMA-IR, homeostatic model assessment-insulin resistance.

Additional file 3

: Tables S1. Comparison of participants' baseline characteristics with or without additional spirometry and Table S2. Additional subgroup analysis for the effects of HOMA-IR on annual change in lung function.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kim, S.H., Kim, H.S., Min, H.K. et al. Association between insulin resistance and lung function trajectory over 4 years in South Korea: community-based prospective cohort. BMC Pulm Med 21, 110 (2021). https://doi.org/10.1186/s12890-021-01478-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12890-021-01478-7

Keywords