Journal of Clinical Medicine Research, ISSN 1918-3003 print, 1918-3011 online, Open Access
Article copyright, the authors; Journal compilation copyright, J Clin Med Res and Elmer Press Inc
Journal website https://jocmr.elmerjournals.com

Original Article

Volume 18, Number 5, May 2026, pages 326-335


Obesity Severity Differentially Shapes Diabetes-Related Impairment in Cardiorespiratory Fitness: A Cross-Sectional Propensity Score–Weighted Analysis of Middle-Aged Adults

Figure

↓  Figure 1. Study enrollment flowchart.
Figure 1.

Tables

↓  Table 1. Baseline Characteristics of Participants According to Obesity–Diabetes Phenotype
 
VariableNone DMDMP value
Moderate obesity (24 ≤ BMI < 35) (n = 213)Severe obesity (BMI ≥ 35) (n = 190)Moderate obesity (24 ≤ BMI < 35) (n = 215)Severe obesity (BMI ≥ 35) (n = 205)
Values are presented as mean ± standard deviation or number (percentage). P values were calculated using one-way ANOVA for continuous variables and the χ2 test for categorical variables. Treatment includes antihypertensive and cardiometabolic medications (β-blockers, ACE inhibitors, CCBs, and diuretics), and does not refer exclusively to antidiabetic therapy. aParticipants may receive more than one type of medication; therefore, counts for individual medication categories are not mutually exclusive. ACE: angiotensin-converting enzyme; ANOVA: analysis of variance; BMI: body mass index; CCB: calcium channel blocker; DM: diabetes mellitus; HbA1c: glycated hemoglobin; HDL: high-density lipoprotein; LDL: low-density lipoprotein.
Age (years)29.8 ± 8.028.5 ± 8.334.9 ± 10.731.7 ± 9.0< 0.001
Female sex, n (%)72 (28)101 (45)204 (74)120 (48)< 0.001
Height (cm)168.6 ± 9.6169.4 ± 8.8169.8 ± 8.7169.1 ± 8.00.49
Body weight (kg)88.4 ± 11.0118.3 ± 19.282.3 ± 12.5116.1 ± 18.1< 0.001
BMI (kg/m2)31.1 ± 2.841.1 ± 4.928.5 ± 3.040.5 ± 4.6< 0.001
Lean body mass (kg)53.5 ± 10.364.1 ± 12.756.7 ± 9.463.7 ± 11.8< 0.001
Duration of diabetes (years)--3.3 ± 4.51.0 ± 2.0-
Family history of diabetes, n (%)24 (9)0 (0)96 (35)45 (20)< 0.001
Fasting plasma glucose (mmol/L)4.6 ± 0.54.8 ± 0.48.6 ± 2.98.4 ± 3.0< 0.001
2-h postprandial glucose (mmol/L)5.9 ± 1.67.3 ± 1.615.5 ± 4.214.0 ± 4.1< 0.001
HbA1c (%)5.3 ± 0.35.5 ± 0.38.9 ± 2.58.6 ± 3.9< 0.001
Triglycerides (mmol/L)1.5 ± 0.71.9 ± 0.82.5 ± 2.02.6 ± 2.0< 0.001
Total cholesterol (mmol/L)4.1 ± 0.85.7 ± 0.94.6 ± 1.15.2 ± 1.00.040
HDL cholesterol (mmol/L)1.2 ± 0.31.1 ± 0.31.1 ± 0.21.0 ± 0.2< 0.001
LDL cholesterol (mmol/L)2.4 ± 0.62.7 ± 0.82.5 ± 0.92.9 ± 0.90.70
Hypertension, n (%)48 (18)78 (35)124 (45)130 (58)< 0.001
Pulmonary diseases, n (%)148 (57)155 (70)196 (71)105 (47)< 0.001
Current smoking, n (%)0 (0)8 (4)56 (20)20 (9)< 0.001
Treatment (cardiometabolic medications)a, n (%)76 (29)64 (29)76 (28)60 (24)0.555
β-blocker, n (%)26 (10)18 (8)24 (9)0 (0)< 0.001
ACE inhibitor, n (%)39 (15)31 (14)28 (10)15 (7)0.018
CCB, n (%)18 (7)15 (7)32 (12)50 (22)< 0.001
Diuretic, n (%)31 (12)28 (13)28 (10)10 (4)0.015

 

↓  Table 2. Unadjusted Cardiopulmonary Exercise Variables Across Obesity–Diabetes Phenotypes
 
VariableNone DMDMP value
Moderate obesity (24 ≤ BMI < 35) (n = 213)Severe obesity (BMI ≥ 35) (n = 190)Moderate obesity (24 ≤ BMI < 35) (n = 215)Severe obesity (BMI ≥ 35) (n = 205)
Values are presented as mean ± standard deviation. Variables with missing data exceeding 10% (peak respiratory rate, breathing reserve, and peak blood pressure) were not included as covariates in the propensity score weighting models. Missing values were handled using multiple imputation, and sensitivity analyses comparing original and imputed datasets yielded consistent results. P values were calculated using one-way analysis of variance or the Kruskal–Wallis test, as appropriate. P values were calculated using the original dataset. For most variables, results were consistent after multiple imputation. Breathing reserve was the only variable showing discrepant significance (original dataset: P = 0.191; after imputation: P < 0.001) and should therefore be interpreted with caution. Implausible values were identified based on predefined physiological ranges and set to missing prior to analysis. Peak VO2: maximal oxygen consumption; RER: respiratory exchange ratio; VE: minute ventilation; VCO2: carbon dioxide output; FEV1: forced expiratory volume in 1 second; FVC: forced vital capacity; MVV: maximal voluntary ventilation.
MVV (L/min)93.5 ± 34.1103.8 ± 28.297.4 ± 30.896.0 ± 29.70.004
Forced vital capacity (L)3.6 ± 1.13.8 ± 0.93.6 ± 0.83.5 ± 0.7< 0.001
FEV1 (L)3.0 ± 0.93.1 ± 0.72.9 ± 0.72.9 ± 0.70.004
FEV1/FVC (%)0.8 ± 0.10.8 ± 0.10.8 ± 0.10.8 ± 0.10.355
Resting heart rate (beats/min)93.5 ± 14.195.2 ± 12.786.0 ± 12.896.1 ± 11.5< 0.001
Resting systolic blood pressure (mm Hg)122.7 ± 17.1128.2 ± 21.2131.2 ± 26.5138.3 ± 14.8< 0.001
Resting diastolic blood pressure (mm Hg)83.8 ± 10.685.8 ± 16.182.2 ± 12.292.6 ± 10.1< 0.001
Peak RER1.1 ± 0.11.1 ± 0.11.1 ± 0.11.1 ± 0.1< 0.001
Peak VO2 (L/min)1.7 ± 0.42.0 ± 0.51.4 ± 0.42.0 ± 0.5< 0.001
Peak VO2 (mL/kg/min)18.5 ± 3.816.9 ± 3.317.9 ± 3.516.8 ± 3.0< 0.001
Peak metabolic equivalents (METs)5.7 ± 1.15.1 ± 0.95.2 ± 1.05.1 ± 0.9< 0.001
Peak work rate (W)136.6 ± 37.3158.2 ± 41.9125.7 ± 33.8151.2 ± 40.6< 0.001
Peak work rate (W/kg)1.5 ± 0.31.3 ± 0.31.5 ± 0.31.3 ± 0.3< 0.001
Peak ventilation (L/min)53.3 ± 12.361.0 ± 16.248.8 ± 12.858.5 ± 15.3< 0.001
Peak VCO2 (L/min)2.0 ± 0.52.3 ± 0.61.8 ± 0.42.3 ± 0.6< 0.001
Peak respiratory rate (breaths/min)32.6 ± 5.634.6 ± 8.330.1 ± 5.734.4 ± 8.3< 0.001
Breathing reserve (%)69.3 ± 20.967.3 ± 19.570.9 ± 20.166.7 ± 22.20.191
VE/VCO2 slope26.7 ± 2.326.3 ± 2.327.3 ± 2.926.0 ± 2.9< 0.001
VO2 at anaerobic threshold (L/min)1.1 ± 0.21.3 ± 0.31.0 ± 0.21.3 ± 0.3< 0.001
Peak O2 pulse (mL/beat)11.1 ± 2.714.6 ± 5.110.5 ± 2.213.4 ± 3.0< 0.001
ΔVO2/ΔWork rate (mL/min/W)12.3 ± 1.512.8 ± 1.511.7 ± 1.813.2 ± 1.8< 0.001
Peak heart rate (beats/min)152.7 ± 15.2149.6 ± 15.7144.7 ± 21.4147.0 ± 12.8< 0.001
Percent predicted maximal heart rate (%)80 ± 1181 ± 1480 ± 580 ± 100.001
Peak systolic blood pressure (mm Hg)176.3 ± 25.8195.8 ± 26.0179.4 ± 34.3193.3 ± 27.7< 0.001
Peak diastolic blood pressure (mm Hg)85.1 ± 18.789.8 ± 17.987.3 ± 17.493.2 ± 14.4< 0.001

 

↓  Table 3. Propensity Score–Weighted Outcome Analyses: PSW-Adjusted Main Effects on Peak VO2 (L/min)
 
GroupPeak VO2 (mean ± SE)
Data are presented as mean ± standard error (SE), predicted values with 95% confidence intervals, or regression coefficients, as appropriate. All analyses were conducted using propensity score weighting to balance baseline covariates across obesity–diabetes phenotypes. Predicted values were derived from a propensity score–weighted linear regression model including a BMI × diabetes interaction term (P for interaction < 0.001). Predicted values and regression estimates were derived from the same propensity score–weighted linear regression model including BMI category, T2DM status, their interaction term, and all covariates. Medication classes were encoded as separate binary indicators because treatment categories were not mutually exclusive. BMI: body mass index; DM: diabetes mellitus; PSW: propensity score weighting; VO2: oxygen uptake. The footnote applies to Tables 3–7.
Moderate obesity, no DM1.69 ± 0.03
Severe obesity, no DM1.95 ± 0.04
Moderate obesity, DM1.44 ± 0.02
Severe obesity, DM1.95 ± 0.03

 

↓  Table 4. Propensity Score–Weighted Outcome Analyses: BMI × DM Interaction: Predicted Peak VO2 (L/min)
 
BMI categoryDiabetes statusPredicted peak VO2 (95% CI)
Moderate obesityNo DM1.54 (1.50–1.58)
Severe obesityNo DM1.74 (1.68–1.80)
Moderate obesityDM1.12 (1.05–1.18)
Severe obesityDM1.72 (1.67–1.77)

 

↓  Table 5. Propensity Score–Weighted Outcome Analyses: Full PSW-Adjusted Interaction Model for Peak VO2 (L/min)
 
TermEstimate (β)Std. error95% CIP value
Intercept0.30730.1308(0.0509, 0.5637)0.019
Severe obesity0.01870.0339(−0.0478, 0.0851)0.582
T2DM−0.21860.0412(−0.2993, −0.1379)< 0.001
Age−0.00190.0014(−0.0048, 0.0009)0.186
Sex (male)−0.04820.0523(−0.1507, 0.0543)0.357
Smoking−0.07990.0530(−0.1838, 0.0240)0.132
Hypertension0.00090.0273(−0.0526, 0.0544)0.974
β-blocker use (yes)−0.01110.044(−0.0974, 0.0752)0.801
ACEI use (yes)−0.01210.0337(−0.0781, 0.0539)0.719
CCB use (yes)0.04850.0293(−0.0089, 0.1059)0.098
Diuretic use (yes)−0.03320.033(−0.1578, 0.0286)0.062
Diabetes duration−0.01450.0061(−0.0265, −0.0026)0.017
Lean body mass0.02710.0025(0.0222, 0.0320)< 0.001
Severe obesity × T2DM0.27850.0540(0.1727, 0.3843)< 0.001

 

↓  Table 6. Propensity Score–Weighted Outcome Analyses: Weighted Correlation With Peak VO2 (L/min)
 
VariableWeighted rP value
2-h postprandial glucose (mmol/L)−0.27< 0.001
Lean body mass (kg)0.7< 0.001

 

↓  Table 7. Propensity Score–Weighted Outcome Analyses: Weighted Linear Regression for Peak VO2 (L/min)
 
VariableEstimateStd. errort valueP value
Intercept−0.130.09−1.470.143
2-h postprandial glucose (mmol/L)−0.010−6.63< 0.001
Lean body mass (kg)0.04027.21< 0.001
Age (years)00−1.550.121
Sex (reference: female)----