Front. Immunol. Frontiers in Immunology Front. Immunol. 1664-3224 Frontiers Media S.A. 10.3389/fimmu.2021.758386 Immunology Original Research TNF-α Carried by Plasma Extracellular Vesicles Predicts Knee Osteoarthritis Progression Zhang Xin 1 2 * Hsueh Ming-Feng 1 2 Huebner Janet L. 1 Kraus Virginia B. 1 2 3 1 Duke Molecular Physiology Institute, Duke University School of Medicine, Duke University, Durham, NC, United States 2 Department of Orthopaedic Surgery, Duke University School of Medicine, Duke University, Durham, NC, United States 3 Department of Medicine, Duke University School of Medicine, Duke University, Durham, NC, United States

Edited by: Matthew William Grol, Western University, Canada

Reviewed by: Ahuva Nissim, Queen Mary University of London, United Kingdom; Mark Hamrick, Augusta University, United States

*Correspondence: Xin Zhang, xin.zhang193@duke.edu

This article was submitted to Inflammation, a section of the journal Frontiers in Immunology

 06 10 2021 2021 12 758386 13 08 2021 21 09 2021 Copyright © 2021 Zhang, Hsueh, Huebner and Kraus 2021 Zhang, Hsueh, Huebner and Kraus

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Objectives

To identify plasma extracellular vesicles (EVs) associated with radiographic knee osteoarthritis (OA) progression.

 Methods

EVs of small (SEV), medium (MEV) and large (LEV) sizes from plasma of OA participants (n=30) and healthy controls (HCs, n=22) were profiled for surface markers and cytokine cargo by high-resolution flow cytometry. The concentrations of cytokines within (endo-) and outside (exo-) EVs were quantified by multiplex ELISA. EV associations with knee radiographic OA (rOA) progression were assessed by multivariable linear regression (adjusted for baseline clinical variables of age, gender, BMI and OA severity) and receiver operating characteristic (ROC) curve analysis.

 Results

Based on integrated mean fluorescence intensity (iMFI), baseline plasma MEVs carrying CD56 (corresponding to natural killer cells) predicted rOA progression with highest area under the ROC curve (AUC) 0.714 among surface markers. Baseline iMFI of TNF-α in LEVs, MEVs and SEVs, and the total endo-EV TNF-α concentration, predicted rOA progression with AUCs 0.688, 0.821, 0.821, 0.665, respectively. In contrast, baseline plasma exo-EV TNF-α (the concentration in the same unit of plasma after EV depletion) did not predict rOA progression (AUC 0.478). Baseline endo-EV IFN-γ and exo-EV IL-6 concentrations were also associated with rOA progression, but had low discriminant capacity (AUCs 0.558 and 0.518, respectively).

 Conclusions

Plasma EVs carry pro-inflammatory cargo that predict risk of knee rOA progression. These findings suggest that EV-associated TNF-α may be pathogenic in OA. The sequestration of pathogenic TNF-α within EVs may provide an explanation for the lack of success of systemic TNF-α inhibitors in OA trials to date.

 extracellular vesicles knee osteoarthritis progression immune cells cytokines TNF-α R56AG060895, R01AG070146, P30AG028716 National Institute on Aging10.13039/100000049 Introduction

Extracellular vesicles (EVs) are released by almost all mammalian cells. Due to their cargo (cytoplasmic proteins, DNA, mRNA, miRNA, small non-coding RNAs, mitochondria, and cytokines), EVs are believed to be able to act as mediators of cell-to-cell communication and as paracrine effectors (1–4). Studies in OA have focused on the beneficial effects of mesenchymal stem cell (MSC)-derived EVs (5, 6), the detrimental effects of subchondral bone osteoblast-derived small EVs (SEVs) (7), and the surface markers and cytokine cargo of SEVs (also known as exosomes) in OA synovial fluid (SF) (8, 9). SEVs from knee OA SF carry surface markers CD9, CD81 and CD63 and cytokines (IL-1β, IL-2, IL-4, IL-5, IL-6, IL-13, IL-17, TNF-α and IFN-γ (8–10), and are associated with OA disease severity (8). SEVs from OA SF induce the release of pro-inflammatory cytokines (IL-1β, IL-6 and TNF-α), chemokines and metalloproteases in vitro by M1 macrophages (9), enhance chemotaxis of peripheral blood mononuclear cells, promote inflammatory responses, and inhibit chondrocyte proliferation (10). SEVs from IL-1β stimulated human synovial fibroblasts significantly up-regulate articular chondrocyte expression of MMP-13 and ADAMTS-5 (11). Taken together, these studies suggest a pathogenic role of SF-derived EVs in OA.

Despite the evidence for a role of EVs in the pathogenesis of OA, to our knowledge, there have been no previous studies evaluating the role of EV subpopulations in OA progression. To fill this important knowledge gap, we profiled OA and healthy control (HC) plasma for EV surface markers, ‘endo-EV’ (within EV) cytokine cargo, and exo-EV cytokine concentrations (in EV-depleted supernatant after ExoQuick precipitation of EVs from the unit of plasma), to evaluate and compare their associations with knee radiographic (r)OA severity and progression.

 Methods Study Participants

Sixty plasma specimens were analyzed including: (1) healthy controls (HCs, n=16), non-progressive OA (OA-NP, n=16) and progressive OA (OA-P, n=14) from the completed Genetics of Generalized Osteoarthritis (GOGO) study (12); and (2) additional HC samples (n=6) from the commercial vendor (Zenbio). The samples were matched for gender, race and decade of age ( Table 1 ). Samples were stored at -80°C until analysis. All samples and data were acquired with informed consent under IRB approval of Duke University or the commercial vendor (Zenbio).

Demographic information of the study participants.

HC* HC OA-NP OA-P
Sample number n=6 n=16 n=16 n=14
Mean (SD) age at enrollment (years) 55.5 (11) 68.4 (8) 68.7 (8) 69.3 (8)
Gender (Female) % 50% 50% 50% 57%
Mean (SD) BMI at enrollment (kg/m2) 28 (8) 27 (3) 32 (8) 29 (5)
Median (range) Summed Baseline K/L grade N/A 0 (0) 2 (1-6) 2 (1-6)
Median (range) Summed Change in K/L grade N/A 0 (0) 0 (0) 3 (1-6)

HC, healthy control; *obtained from Zenbio; OA-P, radiographic knee osteoarthritis progressor; OA-NP, radiographic knee osteoarthritis non-progressor; SD: standard deviation; BMI, body mass index; K/L grade, Kellgren and Lawrence grade; N/A, not applicable.

 Radiographic Procedures and Grading

Knee radiographic imaging was performed as reported previously for the GOGO study (12) and scored for Kellgren and Lawrence (K/L) grade (0–4) (12–14). K/L scores from both knees of a participant were summed yielding scores with range 0-8. HCs from GOGO were defined as participants having knee K/L grade 0 bilaterally. Participants having knee OA was defined as having summed K/L grade ≥1 at baseline. The change of K/L scores from both knees of a participant were summed yielding scores with range 0-6. Radiographic knee OA progression was defined as K/L grade increase ≥1 unit in at least one knee during follow-up (mean 3.8 years, range 1.1-8.6 years) ( Table 1 ). HCs from Zenbio were defined as no self-reported diseases or medical conditions.

 EV Separation From Plasma Samples

EV isolation required 50 μl plasma for each marker panel as previously reported (4, 15). Blood samples were centrifuged at 3000 rpm for 15 min at 4 °C to separate plasma from cells and debris; plasma samples were aliquoted and frozen at − 80 °C until analysis. Frozen plasma samples were thawed followed by centrifugation at 2000 g for 10 minutes at 4°C to remove remaining debris. EVs in plasma were separated by ExoQuick (System Biosciences) following the manufacturer’s instructions (4, 15, 16). As described below, EVs were profiled for surface markers and cytokines; endo-EV and exo-EV cytokine concentrations were also measured.

 Profiling EV-Carried Surface Markers and Cytokines by High Resolution Multicolor Flow Cytometry

As previously described (4, 15), EVs were profiled for the following surface markers to identify EV subpopulations by cell of origin ( Supplementary Table 1 ) ( 4, 15, 17–25): CD81, CD9, CD29, CD63, CD8, CD68, CD14, CD56, CD15, CD235a, CD41a, CD34, CD31, major histocompatibility complex (MHC)-class I antigens HLA-A, HLA-B and HLA-C (HLA-ABC), MHC-class II antigens HLA-DR, -DP and -DQ (HLA-DRDPDQ) (BD Biosciences), CD4, CD19 and MHC-class I antigen HLA-G (ThermoFisher Scientific). We also profiled endo-EV cytokines IL-1β, TNF-α, IFN-γ (BD Biosciences), and IL-6 (ThermoFisher Scientific). The percentages (%) and geometric mean fluorescence intensity (MFI) of EVs carrying each tested marker were determined using a high-resolution multicolor BD LSR Fortessa X-20 Flow Cytometer with the BD FACSDiVa software (BD Biosciences). The integrated MFI (iMFI) of surface markers and cytokines was calculated by multiplying percentage of positive population with the MFI of that population (15, 26, 27).

 Multiplex ELISA for Cytokine Quantification

To confirm the efficiency of the EV precipitation, the concentration and size distribution of particles (in the ExoQuick precipitate of EVs and the remaining EV-depleted supernatants) derived from 50 µl plasma were measured by nanoparticle tracking analysis and dynamic light scattering as previously reported (4, 15) ( Supplementary Figure 1 ). As previously described (15), EV pellets were lysed in NP40 lysis buffer (Thermo Fisher Scientific) in the same volume as the EV-depleted supernatants. The concentrations of endo-EV and exo-EV (remaining EV-depleted supernatant after ExoQuick precipitation) cytokines were measured by multiplex ELISA using the Custom Pro-inflammatory Panel (IL-1β, IL-6, TNF-α, IFN-γ, Meso Scale Diagnostics) following the manufacturer’s instructions (15, 28, 29).

 Receiver Operating Characteristic (ROC) Curve Analysis

Multivariable logistic regression and ROC curve analyses (30, 31) were performed to evaluate the discriminant ability of baseline EV-related variables for knee OA progression, defined as change in knee K/L grade (summed across knees). Model stabilities were validated using a 2,500 non-parametric bootstrap resampling approach; 95% bias-corrected confidence intervals (CIs) for area under the ROC curve (AUC) are reported. Likelihood ratio (32) test was used to assess model fit; the R-square value (RSq) (33), corrected Akaike’s Information Criterion (AICs) (34) and the Bayesian Information Criterion (BIC) (35) are reported. Rsq closer to 1 indicates a better fit to the data; while for AIC and BIC, the model having the smaller value is considered better. Specificity was determined at sensitivity 80%. The analyses were performed using JMP Pro 15 (SAS). AUC is interpreted as follows: AUC ≤0.5 indicates no better than a random classifier; AUC >0.5 is considered validated; AUC >0.65 is considered a moderate discriminant capability.

 Statistical Analyses

Data in this study were not normally distributed based on D’Agostino-Pearson omnibus normality test; therefore, nonparametric analyses were performed. Comparisons between HC, OA-NP and OA-P were performed using Kruskal-Wallis test. Comparisons of endo-EV and exo-EV cytokines were performed using Wilcoxon matched-pairs signed rank test. Comparisons between the tested cytokines in each participant group were performed using Friedman test. False Discovery Rate (FDR) was generated using the Benjamini and Yekutieli method with significant results defined by FDR (q value) <0.05. Multivariable linear regression modeling was performed with adjustment for baseline clinical variables (age, gender, body mass index [BMI] and summed knee OA K/L score) to identify associations of endo-EV and exo-EV biomarkers with knee rOA progression, defined by change in summed K/L score from baseline to follow-up. A p value <0.05 was considered statistically significant. GraphPad Prism 8.0 software (GraphPad) and JMP Pro 15 software were used for statistical analyses.

 Results Multiple Immune Cell-Related LEVs and MEVs Were Associated With Knee rOA Progression

There are currently no specific markers for different EV subtypes (36). Following the recommendations from the International Society for Extracellular Vesicles (36), we used operational terms (range of sizes along with descriptions of cell of origin defined by their surface markers) to describe EV subsets as in our previous work (4, 15). With our newly developed high-resolution multicolor flow cytometry-based methodology, we identified three major subsets of plasma EVs in human HCs: large EVs (LEVs), 1000-6000 nm; medium sized EVs (MEVs), 100-1000 nm; and SEVs, <100 nm; these major subsets based on size were confirmed using dynamic light scattering (4, 15). In this study, these three major subsets of plasma EVs were also identified in plasma of participants with knee OA. Plasma EVs from HC and OA participants all carried surface markers of human stem cells and progenitor cells, immune cells, activated pro-inflammatory fibroblasts, epithelial and endothelial cells indicative of their cell origins ( Supplementary Figure 2 and Supplementary Table 1 ) ( 4, 15, 17–25).

Baseline iMFI of multiple plasma EV subpopulations were associated with baseline clinical variables including age, gender, BMI and summed knee OA K/L score ( Supplementary Figure 3 ). Adjusting for these baseline clinical variables, baseline iMFI of multiple plasma LEVs and MEVs were statistically significantly associated with knee rOA progression. LEV subpopulations associated with rOA progression included the following: CD29+, CD63+, CD8+, CD15+, CD14+, CD19+, and ratio of CD15/CD8 (reflecting neutrophil-EV to T cell-EV ratio) ( Figure 1A ). MEV subpopulations associated with rOA progression included the following: CD81+, CD9+, CD31+, CD29+, CD63+, CD56+, CD68+, and HLA-DRDPDQ+ ( Figure 1B ). No SEV subpopulations based on EV surface markers were associated with knee rOA progression. These LEV and MEV surface markers yielded AUCs >0.5 ( Figures 1C, D ); CD56+ MEVs (AUC 0.714), corresponding to NK cell-EVs, yielded the highest AUC ( Figure 1D ).

Multiple immune cell-associated LEVs and MEVs at baseline were associated with knee radiographic (r)OA progression. Plasma EVs from OA participants (n = 30, 53% OA-NP and 47% OA-P) at baseline were profiled with the indicated surface markers by high resolution multicolor flow cytometry. (A, B) Multivariable linear regression modeling was performed with adjustment for the baseline clinical variables (age, gender, BMI and summed knee OA K/L score) to identify the associations of the knee OA radiographic severity changes from baseline to follow-up with the baseline iMFI of each surface marker in gated LEVs, MEVs or SEVs. The Leverage Plots represent plots of leverage of summed knee rOA progression (x axis, scaled in units of the response) vs the baseline iMFI (y axis) of the indicated surface markers in LEVs (A) and MEVs (B) with an asterisks indicate the p < 0.05. No SEV subpopulations based on EV surface markers were associated with knee rOA progression. (C, D) Multivariable logistic regression and ROC curve analysis were used to assess AUC of the knee rOA progression-associated surface markers in LEVs (C) and MEVs (D) to predict progression of radiographic knee OA. Since all tested surface markers in SEVs were not significantly associated with knee rOA progression, they were excluded from ROC analysis. AUC > 0.5 was considered validated and AUC > 0.65 (red font) was considered a moderate discriminant capability. 95% bias-corrected confidence intervals (CIs) for AUCs are reported. LEVs, large EVs; MEVs, medium EVs; iMFI, integrated mean fluorescence intensity. *p < 0.05.

 Endo-EV TNF-EV TNF-α in Plasma Associated With Knee OA Progression

TNF-α is a classical pro-inflammatory cytokine playing critical roles in OA pathogenesis (29, 37–41). Recently, we found that endo-EV concentrations of TNF-α (in lysates of EV pellets) were significantly higher than the exo-EV concentrations (in EV-depleted supernatants) in the same unit volume of plasma from OA participants, suggesting that EVs are a major source of systemic TNF-α(15). Here, we report, the concentration of endo-EV TNF-α was significantly higher in OA-P than HC participants, while exo-EV TNF-α concentrations did not differ among the study groups ( Figure 2A ). Adjusting for baseline clinical variables, the baseline concentration of endo-EV TNF-α was statistically significantly associated with knee rOA progression ( Figure 2B ). In addition, the baseline concentration of endo-EV TNF-α predicted knee rOA progression with moderate discriminant ability (AUC 0.665, Figure 2C ). In contrast, exo-EV TNF-α was neither associated with nor predictive of knee rOA progression ( Figures 2B, C ). In addition to EVs being a major source of systemic TNF-α, this suggests that plasma endo-EV TNF-α represents a promising systemic biomarker for predicting risk of knee rOA progression.

The concentration of endo-EV TNF-α at baseline was associated with knee rOA progression. The concentrations of exo-EV and endo-EV TNF-α in plasma from HC (n = 16), OA-NP (n = 16) and OA-P (n = 14) participants at baseline were measured by multiplex immunoassay. (A) The graphs represent the differential concentrations of endo-EV and exo-EV TNF-α between HC, OA-NP and OA-P. Comparisons between HC, OA-NP and OA-P were performed using Kruskal-Wallis test with significant results defined by FDR q< 0.05; asterisks indicate the p value as * <0.05. (B) Multivariable linear regression modeling was performed with adjustment for the baseline clinical variables to identify the associations of summed knee rOA progression with baseline concentrations of TNF-α in OA participants (n = 30, 53% OA-NP and 47% OA-P). The Leverage Plots represent plots of leverage of summed knee rOA progression (x axis) vs the baseline concentrations (y axis) of endo-EV and exo-EV TNF-α. (C) Multivariable logistic regression and ROC curve analysis were used to assess AUC and specificity (at sensitivity 80%) of endo-EV and exo-EV TNF-α in OA participants (n = 30, 53% OA-NP and 47% OA-P) to predict progression of radiographic knee OA. Rsq, AICs, and BICs were based on likelihood ratio test. 95% bias-corrected confidence intervals (CIs) for AUCs are reported. Endo-EV, within EV; Exo-EV, outside EV; HC, healthy control; OA-P, radiographic knee osteoarthritis progressor; OA-NP, radiographic knee osteoarthritis non-progressor.

Plasma EVs from HCs, OA-NP and OA-P participants were assessed by flow cytometry for an array of cytokines including TNF-α. TNF-α was present in all sizes of plasma EVs and was the most abundant endo-EV cytokine of the four analyzed and the most abundant cytokine in all three participant groups based on mean percentage of TNF-α+ LEVs, MEVs and SEVs ( Supplementary Figure 4 ). Baseline iMFI of TNF-α in LEVs was significantly higher in OA-P than HCs ( Figure 3A ), which is consistent with the differential concentration of endo-EV TNF-α observed between these two groups by ELISA based analyses ( Figure 2A ). Moreover, baseline iMFI of TNF-α in MEVs and SEVs was significantly higher in OA-P than OA-NP ( Figure 3A ). With adjustment for baseline clinical variables, baseline iMFI of TNF-α in EVs of all sizes was associated with knee rOA progression (P<0.1); this was statistically significant for TNF-α in MEVs ( Figure 3B ). Although the baseline iMFI of TNF-α in EVs of all sizes predicted knee rOA progression, MEVs and SEVs had higher AUCs and specificity (both AUC 0.821 and specificity 81% at sensitivity 80%) than LEVs (AUC 0.688 and specificity 63% at sensitivity 80%) ( Figure 3C ).

The iMFI of endo-EV TNF-α at baseline was associated with knee rOA progression. Plasma EVs from HC (n = 6), OA-NP (n = 16) and OA-P (n = 14) participants at baseline were profiled for intra-vesicle TNF-α by high-resolution multicolor flow cytometry. (A) The graphs represent the iMFI of TNF-α in gated LEVs, MEVs and SEVs between HC, OA-NP and OA-P. Comparisons between HC, OA-NP and OA-P were performed using Kruskal-Wallis test with significant results defined by FDR q< 0.05; asterisks indicate the p value as * <0.05 and ** <0.01. (B) Multivariable linear regression modeling was performed with adjustment for the baseline clinical variables to identify the associations of the knee OA radiographic severity changes from baseline to follow-up with the baseline iMFI of TNF-α in OA participants (n = 30, 53% OA-NP and 47% OA-P). The Leverage Plots represent plots of leverage of summed knee rOA progression (x axis) vs the baseline iMFI (y axis) of TNF-α in LEVs, MEVs and SEVs. (C) Multivariable logistic regression and ROC curve analysis were used to assess AUC and specificity (at sensitivity 80%) of TNF-α in LEVs, MEVs and SEVs in OA participants (n = 30, 53% OA-NP and 47% OA-P). Rsq, AICs, and BICs were based on likelihood ratio test. 95% bias-corrected CIs for AUCs are reported. LEVs, large EVs; MEVs, medium EVs; SEVs, small EVs; iMFI, integrated mean fluorescence intensity. HC, healthy control obtained from Zenbio; OA-P, radiographic knee osteoarthritis progressor; OA-NP, radiographic knee osteoarthritis non-progressor.

 Other Pro-Inflammatory Cytokines Associated With Knee rOA Progression

Endo-EV IFN-γ concentrations were significantly higher in OA-P than both OA-NP and HC participants; exo-EV IFN-γ concentrations did not differ among the study groups ( Figure 4A ). In contrast, exo-EV IL-6 concentrations were significantly higher in OA-P than both OA-NP and HC participants; endo-EV IL-6 concentrations did not differ among the study groups ( Figure 4B ). There were no statistically significant differences between the study groups based on concentrations of endo-EV or exo-EV IL-1β, or iMFIs of IL-1β, IL-6 or IFN-γ in EVs of all sizes (data not shown). Adjusting for baseline clinical variables, baseline concentrations of endo-EV IFN-γ ( Figure 4C ) and exo-EV IL-6 ( Figure 4D ) were statistically significantly associated with knee rOA progression and yielded validated AUCs of 0.558, and 0.518, respectively.

The concentrations of endo-EV IFN-γ and exo-EV IL-6 at baseline were associated with knee rOA progression. The concentrations of exo-EV and endo-EV IFN-γ and IL-6 in plasma from HC (n = 16), OA-NP (n = 16) and OA-P (n = 14) participants at baseline were measured by multiplex immunoassay. (A, B) The graphs represent the differential concentrations of endo-EV and exo-EV IFN-γ (A) and IL-6 (B) between HC, OA-NP and OA-P. Comparisons between HC, OA-NP and OA-P were performed using Kruskal-Wallis test with significant results defined by FDR q < 0.05; asterisks indicate the p value as * < 0.05. (C, D) Multivariable linear regression modeling was performed with adjustment for the baseline clinical variables to identify the associations of summed knee rOA progression from baseline to follow-up with the baseline concentrations of IFN-γ and IL-6 in OA participants (n = 30, 53% OA-NP and 47% OA-P). The Leverage Plots represent plots of leverage of summed knee rOA progression (x axis) vs the baseline concentrations (y axis) of endo-EV and exo-EV IFN-γ (C) and IL-6 (D). Endo-EV, within EV; Exo-EV, outside EV; HC, healthy control; OA-P, radiographic knee osteoarthritis progressor; OA-NP, radiographic knee osteoarthritis non-progressor.

 Discussion

All major types of immune cells (neutrophils, macrophages, NK cells, T cells, B cells, dendritic cells, and granulocytes) actively infiltrate OA synovial tissues, releasing cytokines and EVs into the joint that may contribute to structural progression of OA (11, 42–45). In this study, we identified plasma EVs from OA participants carrying surface markers indicative of multiple cell origins including human stem cells and progenitor cells, all major types of immune cells, activated pro-inflammatory fibroblasts, epithelial and endothelial cells (4, 15, 17–25). We identified plasma EVs in OA that carry the major pro-inflammatory cytokines, TNF-α, IFN-γ, IL-6 and IL-1β, demonstrating their pro-inflammatory phenotype. We found that the concentration of endo-EV TNF-α and the iMFI of TNF-α in EV subsets in OA plasma were associated with and strong predictors of knee rOA progression, while plasma exo-EV TNF-α was neither associated with nor predictive of knee rOA progression. Based on this prior work, we know that these TNF-α+ EVs are present and abundant in OA synovial fluid (15). This suggests that EV associated TNF-α may play a role in OA progression.

TNF-α is a classical pro-inflammatory cytokine produced by a broad variety of cell types, including, but not limited to, macrophages, monocytes, T cells, B cells, NK cells, mast cells, keratinocytes, astrocytes, microglial cells, muscle cells, intestinal paneth cells, tumor cells, synoviocytes and articular chondrocytes (46–48). Based on human in vitro data, TNF-α may play a critical role in the pathogenesis of OA (38, 40, 49). Previous human studies reported that serum TNF-α concentrations predicted knee rOA progression (50), and were associated with joint space narrowing (51). SF TNF-α concentrations were associated with knee pain (29). In addition, TNF-α polymorphisms were associated with susceptibility to OA in a Korean population (52).

Although in an animal study, TNF-α inhibition significantly decreased pro-inflammatory cytokines (IL-1β, IL-17a and IL-8), MMP-3 and MMP-9, inflammatory cell infiltration and bone destruction in joints and cartilage of rats with OA (53), the reported human trials of TNF-α inhibitors for OA are limited to three and all have been negative. These three trials (n=60 hand OA, n=84 hand OA, and n=20 knee OA) all used human TNF antibody, adalimumab (54). Only a single case report indicated that neutralizing TNF-α using adalimumab (40 mg subcutaneously every other week) decreased pain and improved joint function, visibly decreased synovial effusion and synovitis and bone marrow edema, and dramatically decreased nocturnal pain and improved walking distance (48). One possible explanation for these conflicting results is that TNF-α binds to two types of receptors with opposing functions — TNFRI mediating pro-inflammatory signals, and TNFRII mediating anti-inflammatory signals. Therefore, neutralization of TNF-α could result in a mixed clinical response, similar to prior observations resulting from treatments targeting both receptors: inhibiting signaling of both TNFRI and TNFRII reduced collagen-induced arthritis, but increased pro-inflammatory cytokine levels and reduced Treg cell activity. In contrast, selective blockade of TNFRI signaling reduced collagen-induced arthritis without the major side effects observed with both TNFRI and TNFRII blockade (55, 56). Another possible explanation is that TNF-α carried by EVs is sequestered from systemic anti-TNF-α. For flow cytometric quantification of intra-vesicular cytokines including TNF-α, it is necessary to fix and permeabilize the EVs; thus, permeability to antibodies is not expected to be a property of circulating EVs in vivo. We observed that the concentration of endo-EV TNF-α was higher than the matched exo-EV TNF-α in the same unit volume of OA plasma. Therefore, neutralizing soluble TNF-α in blood may not be sufficient when a large amount of pathogenic TNF-α is sequestered in EVs and potentially inaccessible to the neutralizing antibody. Nevertheless, these EVs would be expected to be able to fuse their membrane to the plasma membrane of specific target cells, followed by discharge of their luminal cargo [e.g. TNF-α) to the cytoplasm of their target cell (57)]. Given the dual biological functions of TNF-α in arthritis mediated by TNFRI and TNFRII, and potential technical difficulties depleting EV-carried TNF-α in vivo, neutralizing TNF-α with biologics may not be an ideal therapeutic approach for treating OA. Instead, it may be important to lower the production of TNF-α+ EVs. Monitoring of TNF-α+ EVs in plasma could also be a promising companion diagnostic for OA.

Similar to TNF-α, the concentration of endo-EV, but not exo-EV, IFN-γ was significantly associated with knee OA progression, although it was not as strong a predictor of knee OA progression. In contrast, the concentration of exo-EV, but not endo-EV, IL-6 was significantly associated with knee OA progression. These findings suggest that endo-EV and exo-EV pro-inflammatory cytokines may play different roles in the pathogenesis and worsening of OA and represent different biological processes in OA progression.

Among the test surface markers, CD56+ MEVs, related to NK cells, were the strongest predictor of knee rOA progression with the highest AUC. NK cells are one of the principal leukocyte subsets that infiltrate OA synovia (58); NK cells in both peripherial blood and SF of patients with OA produce pro-inflammatory cytokines,TNF-α, IL-6, and IFN-γ. The frequency of NK cells producing these pro-inflammatory cytokines is significantly higher in OA SF than plasma (59). In vivo antibody-mediated depletion of NK cells ameliorates disease in experimental OA, demonstrating their pathogenicity, which is likely mediated by their ability to promote inflammation and bone destruction (60).

The neutrophil to lymphocyte ratio (NLR) in peripheral blood is already a non-invasive and cost-effective biomarker of various systemic diseases, including cancers, cardiovascular diseases and rheumatologic diseases, and has been shown to be associated with OA severity (61–63). Recently, we found that the ratios of neutrophil EVs to lymphocyte EVs (from T cells, B cells and NK cells), and to T cell EVs, were highly correlated between SF and plasma in OA (15). Here we show that the ratio of neutrophil EVs to T cell EVs was associated with knee OA progression (with a validated AUC). These data taken together suggest that these NLR ratios, based on EV quantification, could serve as systemic biomarkers of OA joint pathology and possibly other comorbidities, such as cardiovascular disease. The breakthrough advantage of determining NLR ratio based on EVs is the ability to use frozen archival samples in contrast to the current method requiring fresh cells.

The first limitation of this exploratory study was the limited number of available samples related to difficulties collecting human specimens of matched gender, race, decade of age, and diversity of OA disease severity and progression. The second limitation was the lack of fresh cells, so we could not directly link EV profiles to the inflammatory phenotype of their parent cells. In addition, OA is often associated with comorbidities that might increase the prevalence of pro-inflammatory EVs. Given the sample sizes, assessing effects of comorbidities was outside the scope of these analyses. However, the GOGO study recruited cognitively intact, ambulatory older adults. In this regard, there is potential selection bias for individuals without major comorbidities other than OA. Although these limitations exist, this study identified several potential new knee OA progression-associated EV biomarkers and predictors worthy of further study as new systemic predictors of knee OA progression. This study also provides insights to optimize therapies targeting TNF-α pathways in OA. Our findings would encourage more studies to explore the roles of these EV biomarkers in OA pathogenesis and disease development.

In summary, we identified several immune cell- and pro-inflammatory cytokine-associated EV subpopulations that were significant independent predictors of knee rOA progression. Among these EV biomarkers, baseline iMFI of TNF-α in MEVs and SEVs was significantly higher in plasma of OA-Progressors than OA-Non-progressors, and the best predictor of knee rOA progression with highest AUC and specificity. Interestingly, baseline iMFI of TNF-α in EVs of all sizes, and total concentration of endo-EV TNF-α all predicted rOA progression. In contrast, baseline exo-EV TNF-α concentration was not associated with nor predictive of rOA progression. These data suggest that EV-associated TNF-α may be pathogenic in OA. Sequestration of large amounts of TNF-α in plasma EVs of OA patients, as shown here, may explain the disappointing results to date of TNF-α inhibitors in OA disease modifying trials.

 Data Availability Statement

The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.

 Ethics Statement

The plasma samples used for this study were acquired under approval of Duke University (n=46) or acquired through a commercial entity, ZenBio (n=6) who provided an assurance, as indicated on their website (https://www.zen-bio.com/products/cells/), that all samples are acquired under donor consent and IRB approval. The patients/participants provided their written informed consent to participate in this study.

 Author Contributions

XZ and VK designed the study. XZ and JH performed the experiments and analyzed the data. M-FH performed ROC analyses. VK performed data analyses for radiographic procedures and grading. XZ and VK drafted the manuscript. All authors contributed to the article and approved the submitted version.

 Funding

This research was funded by the National Institute on Aging at National Institutes of Health, grant numbers R56AG060895 (VK and XZ), R01AG070146 (VK, XZ, and JH) and P30AG028716 (VK and JH).

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Acknowledgments

The authors wish to acknowledge all participants who donated specimens for this study; the Duke Cancer Institute Flow Cytometry Shared Resource for providing access to the LSR Fortessa X-20 Flow Cytometer; and the Duke Human Vaccine Institute Research Flow Cytometry Shared Resource Facility for providing the FCS Express 5 software. Nanoparticle tracking analyses were performed by Marina Sokolsky-Papkov at the Nanomedicines Characterization Core Facility, Center of Nanotechnology in Drug Delivery, UNC School of Pharmacy.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.758386/full#supplementary-material

References Samanta S Rajasingh S Drosos N Zhou Z Dawn B Rajasingh J . Exosomes: New Molecular Targets of Diseases. Acta Pharmacol Sin (2018) 39(4):501–13. doi: 10.1038/aps.2017.162 Thery C Ostrowski M Segura E . Membrane Vesicles as Conveyors of Immune Responses. Nat Rev Immunol (2009) 9(8):581–93. doi: 10.1038/nri2567 Phinney DG Di Giuseppe M Njah J Sala E Shiva S St Croix CM . Mesenchymal Stem Cells Use Extracellular Vesicles to Outsource Mitophagy and Shuttle microRNAs. Nat Commun (2015) 6:8472. doi: 10.1038/ncomms9472 Zhang X Hubal MJ Kraus VB . Immune Cell Extracellular Vesicles and Their Mitochondrial Content Decline With Ageing. Immun Ageing (2020) 17(1):1. doi: 10.1186/s12979-019-0172-9 D’Arrigo D Roffi A Cucchiarini M Moretti M Candrian C Filardo G . Secretome and Extracellular Vesicles as New Biological Therapies for Knee Osteoarthritis: A Systematic Review. J Clin Med (2019) 8(11):1867. doi: 10.3390/jcm8111867 Wu X Wang Y Xiao Y Crawford R Mao X Prasadam I . Extracellular Vesicles: Potential Role in Osteoarthritis Regenerative Medicine. J Orthop Translat (2020) 21:73–80. doi: 10.1016/j.jot.2019.10.012 Wu X Crawford R Xiao Y Mao X Prasadam I . Osteoarthritic Subchondral Bone Release Exosomes That Promote Cartilage Degeneration. Cells (2021) 10(2):251. doi: 10.3390/cells10020251 Gao K Zhu W Li H Ma D Liu W Yu W . Association Between Cytokines and Exosomes in Synovial Fluid of Individuals With Knee Osteoarthritis. Mod Rheumatol (2020) 30(4):758–64. doi: 10.1080/14397595.2019.1651445 Domenis R Zanutel R Caponnetto F Toffoletto B Cifu A Pistis C . Characterization of the Proinflammatory Profile of Synovial Fluid-Derived Exosomes of Patients With Osteoarthritis. Mediators Inflammation (2017) 2017:4814987. doi: 10.1155/2017/4814987 Li Z Li M Xu P Ma J Zhang R . Compositional Variation and Functional Mechanism of Exosomes in the Articular Microenvironment in Knee Osteoarthritis. Cell Transplant (2020) 29:963689720968495. doi: 10.1177/0963689720968495 Kato T Miyaki S Ishitobi H Nakamura Y Nakasa T Lotz MK . Exosomes From IL-1beta Stimulated Synovial Fibroblasts Induce Osteoarthritic Changes in Articular Chondrocytes. Arthritis Res Ther (2014) 16(4):R163. doi: 10.1186/ar4679 Kraus VB Jordan JM Doherty M Wilson AG Moskowitz R Hochberg M . The Genetics of Generalized Osteoarthritis (GOGO) Study: Study Design and Evaluation of Osteoarthritis Phenotypes. Osteoarthritis Cartilage (2007) 15(2):120–7. doi: 10.1016/j.joca.2006.10.002 Zuo X Luciano A Pieper CF Bain JR Kraus VB Kraus WE . Combined Inflammation and Metabolism Biomarker Indices of Robust and Impaired Physical Function in Older Adults. J Am Geriatr Soc (2018) 66(7):1353–9. doi: 10.1111/jgs.15393 Daghestani HN Pieper CF Kraus VB . Soluble Macrophage Biomarkers Indicate Inflammatory Phenotypes in Patients With Knee Osteoarthritis. Arthritis Rheumatol (2015) 67(4):956–65. doi: 10.1002/art.39006 Zhang X Huebner JL Kraus VB . Extracellular Vesicles as Biological Indicators and Potential Sources of Autologous Therapeutics in Osteoarthritis. Int J Mol Sci (2021) 22(15):8351. doi: 10.3390/ijms22158351 Eitan E Green J Bodogai M Mode NA Baek R Jorgensen MM . Age-Related Changes in Plasma Extracellular Vesicle Characteristics and Internalization by Leukocytes. Sci Rep (2017) 7(1):1342. doi: 10.1038/s41598-017-01386-z Chou CH Jain V Gibson J Attarian DE Haraden CA Yohn CB . Synovial Cell Cross-Talk With Cartilage Plays a Major Role in the Pathogenesis of Osteoarthritis. Sci Rep (2020) 10(1):10868. doi: 10.1038/s41598-020-67730-y Wright MD Moseley GW van Spriel AB . Tetraspanin Microdomains in Immune Cell Signalling and Malignant Disease. Tissue Antigens (2004) 64(5):533–42. doi: 10.1111/j.1399-0039.2004.00321.x Park KR Inoue T Ueda M Hirano T Higuchi T Maeda M . CD9 Is Expressed on Human Endometrial Epithelial Cells in Association With Integrins Alpha(6), Alpha(3) and Beta(1). Mol Hum Reprod (2000) 6(3):252–7. doi: 10.1093/molehr/6.3.252 Lin KK Rossi L Boles NC Hall BE George TC Goodell MA . CD81 Is Essential for the Re-Entry of Hematopoietic Stem Cells to Quiescence Following Stress-Induced Proliferation via Deactivation of the Akt Pathway. PloS Biol (2011) 9(9):e1001148. doi: 10.1371/journal.pbio.1001148 Iwasaki T Takeda Y Maruyama K Yokosaki Y Tsujino K Tetsumoto S . Deletion of Tetraspanin CD9 Diminishes Lymphangiogenesis In Vivo and In Vitro . J Biol Chem (2013) 288(4):2118–31. doi: 10.1074/jbc.M112.424291 Jin Y Takeda Y Kondo Y Tripathi LP Kang S Takeshita H . Double Deletion of Tetraspanins CD9 and CD81 in Mice Leads to a Syndrome Resembling Accelerated Aging. Sci Rep (2018) 8(1):5145. doi: 10.1038/s41598-018-23338-x Privratsky JR Newman PJ . PECAM-1: Regulator of Endothelial Junctional Integrity. Cell Tissue Res (2014) 355(3):607–19. doi: 10.1007/s00441-013-1779-3 Sarangi PP Hyun YM Lerman YV Pietropaoli AP Kim M . Role of Beta1 Integrin in Tissue Homing of Neutrophils During Sepsis. Shock (2012) 38(3):281–7. doi: 10.1097/SHK.0b013e31826136f8 Togarrati PP Dinglasan N Desai S Ryan WR Muench MO . CD29 Is Highly Expressed on Epithelial, Myoepithelial, and Mesenchymal Stromal Cells of Human Salivary Glands. Oral Dis (2018) 24(4):561–72. doi: 10.1111/odi.12812 Shooshtari P Fortuno ES 3rd Blimkie D Yu M Gupta A Kollmann TR . Correlation Analysis of Intracellular and Secreted Cytokines via the Generalized Integrated Mean Fluorescence Intensity. Cytomet A (2010) 77(9):873–80. doi: 10.1002/cyto.a.20943 Brummelman J Mazza EMC Alvisi G Colombo FS Grilli A Mikulak J . High-Dimensional Single Cell Analysis Identifies Stem-Like Cytotoxic CD8(+) T Cells Infiltrating Human Tumors. J Exp Med (2018) 215(10):2520–35. doi: 10.1084/jem.20180684 Hsueh MF Zhang X Wellman SS Bolognesi MP Kraus VB . Synergistic Roles of Macrophages and Neutrophils in Osteoarthritis Progression. Arthritis Rheumatol (2020) 73(1):89–99. doi: 10.1002/art.41486 Leung YY Huebner JL Haaland B Wong SBS Kraus VB . Synovial Fluid Pro-Inflammatory Profile Differs According to the Characteristics of Knee Pain. Osteoarthritis Cartilage (2017) 25(9):1420–7. doi: 10.1016/j.joca.2017.04.001 Zweig MH Campbell G . Receiver-Operating Characteristic (ROC) Plots: A Fundamental Evaluation Tool in Clinical Medicine. Clin Chem (1993) 39(4):561–77. doi: 10.1093/clinchem/39.4.561 Kraus VB Feng S Wang S White S Ainslie M Brett A . Trabecular Morphometry by Fractal Signature Analysis Is a Novel Marker of Osteoarthritis Progression. Arthritis Rheumatol (2009) 60(12):3711–22. doi: 10.1002/art.25012 Wilks SS . The Large-Sample Distribution of the Likelihood Ratio for Testing Composite Hypotheses. Ann Math Statist (1938) 9(1):60–2. doi: 10.1214/aoms/1177732360 McFadden D . Conditional Logit Analysis of Qualitative Choice Behavior. Berkeley, Calif: Univ. of California (1973). Akaike H . A New Look at the Statistical Model Identification. IEEE Trans Automatic Control (1974) 19(6):716–23. doi: 10.1109/TAC.1974.1100705 Burnham KP Anderson DR . Multimodel Inference : Understanding AIC and BIC in Model Selection. Sociologic Methods Res (2004) 33(2):261–304. doi: 10.1177/0049124104268644 Thery C Witwer KW Aikawa E Alcaraz MJ Anderson JD Andriantsitohaina R . Minimal Information for Studies of Extracellular Vesicles 2018 (MISEV2018): A Position Statement of the International Society for Extracellular Vesicles and Update of the MISEV2014 Guidelines. J Extracell Vesicles (2018) 7(1):1535750. doi: 10.1080/20013078.2018.1535750 Westacott CI Sharif M . Cytokines in Osteoarthritis: Mediators or Markers of Joint Destruction? Semin Arthritis Rheum (1996) 25(4):254–72. doi: 10.1016/s0049-0172(96)80036-9 Lopez-Armada MJ Carames B Lires-Dean M Cillero-Pastor B Ruiz-Romero C Galdo F . Cytokines, Tumor Necrosis Factor-Alpha and Interleukin-1beta, Differentially Regulate Apoptosis in Osteoarthritis Cultured Human Chondrocytes. Osteoarthritis Cartilage (2006) 14(7):660–9. doi: 10.1016/j.joca.2006.01.005 Kennedy MI Whitney K Evans T LaPrade RF . Platelet-Rich Plasma and Cartilage Repair. Curr Rev Musculoskelet Med (2018) 11(4):573–82. doi: 10.1007/s12178-018-9516-x Bondeson J Wainwright SD Lauder S Amos N Hughes CE . The Role of Synovial Macrophages and Macrophage-Produced Cytokines in Driving Aggrecanases, Matrix Metalloproteinases, and Other Destructive and Inflammatory Responses in Osteoarthritis. Arthritis Res Ther (2006) 8(6):R187. doi: 10.1186/ar2099 Wojdasiewicz P Poniatowski LA Szukiewicz D . The Role of Inflammatory and Anti-Inflammatory Cytokines in the Pathogenesis of Osteoarthritis. Mediators Inflammation (2014) 2014:561459. doi: 10.1155/2014/561459 Caruso S Poon IKH . Apoptotic Cell-Derived Extracellular Vesicles: More Than Just Debris. Front Immunol (2018) 9:1486. doi: 10.3389/fimmu.2018.01486 de Lange-Brokaar BJ Ioan-Facsinay A van Osch GJ Zuurmond AM Schoones J Toes RE . Synovial Inflammation, Immune Cells and Their Cytokines in Osteoarthritis: A Review. Osteoarthritis Cartilage (2012) 20(12):1484–99. doi: 10.1016/j.joca.2012.08.027 Deligne C Casulli S Pigenet A Bougault C Campillo-Gimenez L Nourissat G . Differential Expression of Interleukin-17 and Interleukin-22 in Inflamed and non-Inflamed Synovium From Osteoarthritis Patients. Osteoarthritis Cartilage (2015) 23(11):1843–52. doi: 10.1016/j.joca.2014.12.007 Li YS Luo W Zhu SA Lei GH . T Cells in Osteoarthritis: Alterations and Beyond. Front Immunol (2017) 8:356. doi: 10.3389/fimmu.2017.00356 Vassalli P . The Pathophysiology of Tumor Necrosis Factors. Annu Rev Immunol (1992) 10:411–52. doi: 10.1146/annurev.iy.10.040192.002211 Fernandes JC Martel-Pelletier J Pelletier JP . The Role of Cytokines in Osteoarthritis Pathophysiology. Biorheology (2002) 39(1-2):237–46. Grunke M Schulze-Koops H . Successful Treatment of Inflammatory Knee Osteoarthritis With Tumour Necrosis Factor Blockade. Ann Rheum Dis (2006) 65(4):555–6. doi: 10.1136/ard.2006.053272 Xue J Wang J Liu Q Luo A . Tumor Necrosis Factor-Alpha Induces ADAMTS-4 Expression in Human Osteoarthritis Chondrocytes. Mol Med Rep (2013) 8(6):1755–60. doi: 10.3892/mmr.2013.1729 Bastick AN Belo JN Runhaar J Bierma-Zeinstra SM . What Are the Prognostic Factors for Radiographic Progression of Knee Osteoarthritis? A Meta-Analysis. Clin Orthop Relat Res (2015) 473(9):2969–89. doi: 10.1007/s11999-015-4349-z Stannus O Jones G Cicuttini F Parameswaran V Quinn S Burgess J . Circulating Levels of IL-6 and TNF-Alpha Are Associated With Knee Radiographic Osteoarthritis and Knee Cartilage Loss in Older Adults. Osteoarthritis Cartilage (2010) 18(11):1441–7. doi: 10.1016/j.joca.2010.08.016 Han L Song JH Yoon JH Park YG Lee SW Choi YJ . TNF-Alpha and TNF-Beta Polymorphisms Are Associated With Susceptibility to Osteoarthritis in a Korean Population. Korean J Pathol (2012) 46(1):30–7. doi: 10.4132/KoreanJPathol.2012.46.1.30 Li H Xie S Qi Y Li H Zhang R Lian Y . TNF-Alpha Increases the Expression of Inflammatory Factors in Synovial Fibroblasts by Inhibiting the PI3K/AKT Pathway in a Rat Model of Monosodium Iodoacetate-Induced Osteoarthritis. Exp Ther Med (2018) 16(6):4737–44. doi: 10.3892/etm.2018.6770 Chevalier X Eymard F Richette P . Biologic Agents in Osteoarthritis: Hopes and Disappointments. Nat Rev Rheumatol (2013) 9(7):400–10. doi: 10.1038/nrrheum.2013.44 Bluml S Binder NB Niederreiter B Polzer K Hayer S Tauber S . Antiinflammatory Effects of Tumor Necrosis Factor on Hematopoietic Cells in a Murine Model of Erosive Arthritis. Arthritis Rheumatol (2010) 62(6):1608–19. doi: 10.1002/art.27399 McCann FE Perocheau DP Ruspi G Blazek K Davies ML Feldmann M . Selective Tumor Necrosis Factor Receptor I Blockade Is Antiinflammatory and Reveals Immunoregulatory Role of Tumor Necrosis Factor Receptor II in Collagen-Induced Arthritis. Arthritis Rheumatol (2014) 66(10):2728–38. doi: 10.1002/art.38755 Prada I Meldolesi J . Binding and Fusion of Extracellular Vesicles to the Plasma Membrane of Their Cell Targets. Int J Mol Sci (2016) 17(8):1296. doi: 10.3390/ijms17081296 Huss RS Huddleston JI Goodman SB Butcher EC Zabel BA . Synovial Tissue-Infiltrating Natural Killer Cells in Osteoarthritis and Periprosthetic Inflammation. Arthritis Rheumatol (2010) 62(12):3799–805. doi: 10.1002/art.27751 Jaime P Garcia-Guerrero N Estella R Pardo J Garcia-Alvarez F Martinez-Lostao L . CD56(+)/CD16(-) Natural Killer Cells Expressing the Inflammatory Protease Granzyme A Are Enriched in Synovial Fluid From Patients With Osteoarthritis. Osteoarthritis Cartilage (2017) 25(10):1708–18. doi: 10.1016/j.joca.2017.06.007 Benigni G Dimitrova P Antonangeli F Sanseviero E Milanova V Blom A . CXCR3/CXCL10 Axis Regulates Neutrophil-NK Cell Cross-Talk Determining the Severity of Experimental Osteoarthritis. J Immunol (2017) 198(5):2115–24. doi: 10.4049/jimmunol.1601359 Tasoglu O Boluk H Sahin Onat S Tasoglu I Ozgirgin N . Is Blood Neutrophil-Lymphocyte Ratio an Independent Predictor of Knee Osteoarthritis Severity? Clin Rheumatol (2016) 35(6):1579–83. doi: 10.1007/s10067-016-3170-8 Galizia G Lieto E Zamboli A De Vita F Castellano P Romano C . Neutrophil to Lymphocyte Ratio Is a Strong Predictor of Tumor Recurrence in Early Colon Cancers: A Propensity Score-Matched Analysis. Surgery (2015) 158(1):112–20. doi: 10.1016/j.surg.2015.02.006 Adamstein NH MacFadyen JG Rose LM Glynn RJ Dey AK Libby P . The Neutrophil-Lymphocyte Ratio and Incident Atherosclerotic Events: Analyses From Five Contemporary Randomized Trials. Eur Heart J (2021) 42(9):896–903. doi: 10.1093/eurheartj/ehaa1034