Table 4.

The association of tumor MIR21 expression with the density of T cells

Univariable OR (95% CI)Multivariable OR (95% CI)a
Model for CD3+ cell density (n = 518, as an outcome variable)
MIR21 expressionQ1 (lowest)1 (reference)1 (reference)
Q2 (second)0.88 (0.57–1.36)0.85 (0.55–1.31)
Q3 (third)0.67 (0.43–1.04)0.59 (0.37–0.92)
Q4 (highest)0.47 (0.31–0.73)0.44 (0.28–0.68)
Ptrendb0.0004<0.0001
Model for CD8+ cell density (n = 510, as an outcome variable)
MIR21 expressionQ1 (lowest)1 (reference)1 (reference)
Q2 (second)1.14 (0.74–1.77)1.25 (0.80–1.96)
Q3 (third)0.72 (0.46–1.12)0.76 (0.48–1.19)
Q4 (highest)0.89 (0.58–1.38)0.99 (0.63–1.54)
Ptrendb0.270.50
Model for CD45RO+ cell density (n = 522, as an outcome variable)
MIR21 expressionQ1 (lowest)1 (reference)1 (reference)
Q2 (second)0.70 (0.46–1.09)0.72 (0.46–1.12)
Q3 (third)0.57 (0.37–0.89)0.54 (0.34–0.84)
Q4 (highest)0.45 (0.29–0.70)0.41 (0.26–0.64)
Ptrendb0.0002<0.0001
Model for FOXP3+ cell density (n = 495, as an outcome variable)
MIR21 expressionQ1 (lowest)1 (reference)1 (reference)
Q2 (second)0.98 (0.63–1.54)0.93 (0.59–1.46)
Q3 (third)0.73 (0.47–1.14)0.61 (0.39–0.96)
Q4 (highest)0.66 (0.42–1.02)0.55 (0.35–0.86)
Ptrendb0.0320.003

Abbreviations: Q1 to Q4, quartile 1 to quartile 4.

  • aThe multivariable ordinal logistic regression analysis model initially included age; sex; year of diagnosis; family history of colorectal cancer in parent or sibling; tumor location; tumor differentiation; MSI; CpG island methylator phenotype; KRAS, BRAF, and PIK3CA mutations; and LINE-1 methylation level. A backward stepwise elimination with a threshold of P = 0.05 was used to select variables in the final models. Variables remaining in the final multivariable ordinal logistic regression models are shown in Supplementary Table S2.

  • bPtrend value was calculated by the linear trend across the ordinal (first to fourth quartile) categories of MIR21 expression as a continuous variable in the ordinal logistic regression model for the density of CD3+, CD8+, CD45RO+, or FOXP3+ T cells (an ordinal quartile outcome variable). Because we assessed four primary outcome variables, we adjusted the two-sided α level to 0.012 (= 0.05/4) by simple Bonferroni correction.