sMAdCAM:IL-6 ratio influences disease progression and anti-viral responses in SARS-CoV-2 infection

Study population, setting, and data collection

A total of 60 in-patients and 24 convalescent individuals were recruited (from April-July, 2020), following informed consent, for the study from the isolation ward at Kasturba Hospital for Infectious Diseases (Municipal Corporation of Greater Mumbai). Pregnant women, prisoners, and children (less than 18 years of age) were excluded from the study. The Kasturba Hospital and ICMR-NIRRH institutional ethics committees approved this study. We obtained demographic data, clinical history at presentation, and laboratory results during admission.

Blood samples for the study were handled in accordance with ICMR guidelines for biosafety. Whole blood (1-3 ml) was collected in EDTA vacutainers and plasma was separated by centrifugation at 425 g for 10 minutes. IgG and IgM antibodies against SARS-CoV-2 were detected in fresh plasma samples using Rapid test from either Voxpress (Voxtur Bio LTD, India) or Tell me fast (Biocan Diagnostics Inc., Canada) kits and also by ARCHITECT™ Abbott™ (Abbott Diagnostics, USA) chemiluminescence immunoassay (CLIA) directed against SARS -CoV-2 anti-NC IgG. Remaining plasma samples were aliquoted and stored at −80° C until batch analysis of cytokines, Soluble MAdCAM and antibody binding profile using Surface Plasmon Resonance.

Cytokine Bead Array

Human Th1/Th2/Th17 cytokines kit (BD CBA kit Cat no. 560484, NJ, USA) was used to measure IL-2, IL-4, IL-6, IL-10, TNF, IFN-γ and IL-17 cytokine levels in plasma of study participants. Briefly, 50 µl of plasma sample and standard were added to 50 µl of capture bead mixture followed by addition of 50 µl of Human Th1/Th2/ Th17 PE Detection reagent. Following incubation at room temperature in dark for 3 hours, the assay tubes were washed with 1 ml of wash buffer. The tubes were centrifuged at 200 g for 5 min and bead pellet was re-suspended in 300 µl wash buffer and acquired using BD Accuri™ C6. Five thousand events were acquired gated on selected bead population and analysis was carried out using BD C6 Accuri analysis software. GraphPad Prism was used to extrapolate individual cytokine concentrations in each sample.

MAdCAM ELISA

Human MAdCAM-1 DuoSet ELISA kit (R&D Systems-DY6056-05, Minneapolis, USA) was used for estimation of soluble MAdCAM-1 in plasma in accordance with manufacturer’s protocol. Briefly, wells coated with capture antibody and reagent diluent (1% BSA in PBS) was added for blocking. MAdCAM standard or 1:100 diluted plasma was added to wells and incubated overnight at 4°C. Incubation with detection antibody was followed by addition of streptavidin-HRP which was allowed to react with substrate solution to give chromogenic readout. Reaction was stopped using 2N H₂SO₄ and absorbance measured at 450 nm and 570 nm using Synergy H1 microplate reader (Biotek, Vermont, USA). Readings at 570 nm were subtracted from 450 nm and used to plot standard curve. Plasma concentration of sMAdCAM was calculated by interpolating absorbance values and multiplying by dilution factor.

Surface Plasmon Resonance (SPR) studies: Binding of COVID-19 patients plasma to immobilized SARS-CoV-2 proteins-Spike, Receptor Binding Domain (RBD) and Nucleocapsid

The binding analysis of real time interaction between COVID-19 patients plasma and SARS-CoV-2 proteins was carried out using the SPR spectrometer (Biacore 3000, GE Healthcare Bio-Sciences AB, Uppsala, Sweden). Flow cells of CM5 (carboxmethylated dextran matrix) sensor chip were first activated with 0.2 M 1-ethyl-3-(3-dimethyl aminopropyl) carbodiimide hydrochloride and 0.05 M N-hydroxysuccinimide. SARS-CoV-2 proteins (Spike, RBD and Nucleocapsid) at a concentration of 20 μg/ml each in 10 mM sodium acetate, pH 4.5, was immobilized on the activated CM5 chip to about 1400 response units (RUs), and the unreacted groups were blocked with 1.0 M ethanolamine-HCl, pH 8.5. Flow cell, which was not immobilized with any protein, served as negative control.

Human COVID-19 patients plasma was divided into four groups based on their IgM/IgG status (Rapid Antibody Tests) and was used to determine the binding of antibodies to three SARS-CoV-2 proteins. Diluted plasma (1:10) in 10 mM HBS-EP, pH 7.4 running buffer was injected over Spike, RBD & Nucleocapsid immobilized CM5 chip for 150 seconds (association phase) followed by subsequent dissociation for 150 seconds and recording of their spectra. All the experiments were performed at 25 °C, and the flow rate used was 20 μl/min. The surface at the end of each experiment was regenerated using a 25 μl pulse of 10 mM Glycine-HCL, pH 2.5 to remove any bound analyte. Sensorgrams were corrected by subtraction of the signal from the negative control surface as well as pre COVID-19 plasma.

Data were evaluated using the BIAevaluation software program (4.1 version, GE Healthcare Bio-Sciences AB). The 1:1 Langmuir model was used for fitting the experimental data and to calculate the association response differences and dissociation rates (kd).

Statistical analysis

Statistical analysis was performed in GraphPad Prism 8 using non parametric tests. Statistical significance of differences between groups were assessed using Mann-Whitney U-test. Follow up data from same individual was compared with Wilcoxon matched-paired signed rank test. Spearman’s rank-order correlation was used to analyse the association between participant attributes and statistical significance was accepted at p<0.05.

Demographic and clinical characteristics of in-patients and convalescent individuals

During the period from April through July 2020, 60 in-patients with confirmed Covid-19 infection and under isolation were recruited at Kasturba Hospital. Of these 35 were followed up 7 days post admission. A total of 24 individuals, clinically documented to have been recovered and discharged were also recruited as convalescents. Convalescent individuals had 2 consecutive negative respiratory swabs or were cleared of symptoms prior to discharge and their blood samples were collected atleast 2 weeks post discharge. The demographic, serological and clinical characteristics of the patients are shown and summarized in Supplemental Table 1 and Table 1 respectively. As apparent from these tables our study cohort included individuals that were analysed in a cross-sectional, longitudinal (0 day and 7 day time points) manner and at clinically distinct (in patients and convalescent) stages of disease progression. Also, as expected, males pre-dominated and individuals with co-morbidities constituted a significant proportion (28/66; 42%) of the cohort. Serological data evaluated through both a rapid test (IgG, IgM) and CLIA (anti-NC IgG) showed an overall concordance of 77% for IgG and also allowed us to group these individuals based on IgG and/or IgM positivity. When grouped in this fashion it was interesting to note that while a majority of convalescent individuals ∼88% (21 out of 24) were seropositive for SARS-CoV-2 about 12% (3 out of 24) were seronegative both by CLIA and rapid test (Table 1 and Supplemental Table 1). Also, as a group, convalescent individuals seemed to have lower levels of circulating anti-viral antibodies following infection.

Inflammatory cytokine profile and gut migration

As shown in (Supplemental Figure 1A and 1B), elevated levels of IL-6 and TNF were observed both in-patients as well as, of note, in convalescent individuals compared to healthy individuals. Also, these levels spanned a range of sero-positivity in terms of virus specific IgG and/or IgM. Also, we noticed a significant reduction of IL-6, but not TNF, as seropositivity transitioned from double negative to IgG, IgM double positive and finally to class switched (IgG single positive) responses in convalescent individuals. Increased levels of these cytokines did not correlate with initial Ct values, each other and that of any other cytokine measured (Supplemental Figures 2 and 3). When these levels were examined longitudinally (Supplemental Figure 4) in a group of 11 individuals that were followed at day 0 (day of diagnosis) and day 7 (post diagnosis) we noticed a declining trend for IL-6 but not TNF.

Additionally and for the first time we report on gut inflammatory/migratory marker sMAdCAM (Figure 1 and Supplemental Figure 1C) which appeared to be circulating at levels lower than those seen in healthy individuals and convalescents. Also there was no difference between sMAdCAM levels in 34 individuals that were followed at day 0 and day 7. Interestingly however (Supplemental Figure1c), in contrast and conversely with respect to IL-6, sMAdCAM levels seemed to reduce (compared to healthy controls) following early infection (IgM single positive and IgG, IgM double positive stages) and attain restoration in class switched individuals (IgG single positive and convalescents). Further, LPS levels, evaluated in a subset of in-patients and convalescents, did not vary signifying absence of microbial translocation associated gut inflammation (Supplemental Figure 5).

• Download figure
• Open in new tab

Figure 1:

Dynamic changes in plasma inflammatory markers of SARS-CoV-2 infected individuals and convalescent individuals. Variation in levels of (A) IL-6 at day 0 (n=11), day 7 (n=11), in convalescent (n=22) and in healthy (n=19) group (B) TNF at day 0 (n=11), day 7 (n=11), in convalescent (n=22) and in healthy (n=3) group and (C) Soluble MAdCAM at day 0 (n=35), day 7 (n=35), in convalescent (n=24) and in healthy (n=19) group. Pink coloured symbols indicate IgG-/IgM-individuals in convalescent group. Statistical significance was calculated by Wilcoxon matched-pairs signed rank test and Mann-Whitney U-test; *, p < 0.05; **, p < 0.01; and ***, p<0.001.

Next, analysing these individuals based on disease progression clearly demonstrated a declining circulating IL-6 but not TNF concentration (Figure 1). Stratification of the in-patient data (Supplemental Figure 6) on the basis of gender revealed a marked difference in circulating levels of IL-6 and sMAdCAM where women had significantly lower and higher levels of these respectively.

sMAdCAM:IL-6 ratio (sMIL index): A novel biomarker for disease progression

Results obtained from the heretofore described analyses (Figure 1 and Supplemental Figure 6) suggested that sMAdCAM and IL-6 levels independently seemed to track (Supplemental Figures 7) with disease progression, albeit in opposing directions. Suspecting an underlying association between both of these apparently unrelated signatures, we examined the utility of their being combined as an integrated marker that could better correlate and discriminate between various stages of infection.

Remarkably, when levels of both these analytes were expressed as a ratio we observed (Figure 2) a clear delineation of in-patient and convalescent data into discrete categories based on disease progression (Figure 2A and 2B), gender (Figure 2C), time from infection (Figure 2D) and intriguingly serological status (Figure 2E). These observations seemed to posit a role for this index in the development of antibodies against SARS-CoV-2.

• Download figure
• Open in new tab

Figure 2:

Analysis of soluble MAdCAM and IL-6 ratio (sMIL index) with SARS-CoV2 disease progression A. Variation in sMIL index at day 0 (n=11), day 7 (n=11), in convalescent (n=22) and in healthy (n=19) group. Pink coloured symbols indicate IgG-/IgM-individuals in convalescent group. B. Change in sMIL index between day 0 (n=11) and day 7 (n=11) expressed per individual. C. Differences in sMIL index between male (n=18) and female (n=7). D. Comparison of sMIL index between healthy (n=19),IgG-/IgM-(n=22), IgG+/IgM+ (n=11), IgG+/IgM-(n=14) and convalescent (n=22) groups.. E. Association of sMIL index with days since SARS-CoV-2 confirmation by PCR (n=47). Red coloured symbols represent asymptomatic individuals Statistical analysis was performed using Graphpad Prism 8.0 Wilcoxon matched-pairs signed rank test was used to compare paired samples of Day 0 and Day 7. Mann-Whitney U-test was used to compare unpaired groups. *, p < 0.05; **, p < 0.01, ***, p<0.001; and ****, p<0.001. Correlation analysis was performed using non parametric Spearman Rank Correlation test.

Anti-viral responses and the role of inflammatory markers

Anti-SARS-CoV-2 responses were measured using surface plasmon resonance (SPR) that allowed us to obtain a sensitive and real-time evaluation of antibody binding to viral Spike, RBD and NC proteins. First, we demonstrated that SPR was indeed picking up antibody binding through the observed concordance between CLIA profile and viral protein binding association rates observed in our prospective cohort (Supplemental Figure 8). Further, we observed that individuals that were below the level of detection in the CLIA assay (Y-axis value 0) did have detectable binding using SPR.

We then proceeded to stratify plasma binding data (association rates) based on disease progression (Supplemental Figure 9) as well as IgM/IgG serostatus (Figure 3). In the case of the former we observed significant increases by day 7 for all 3 viral targets with increased responses against NC and RBD compared to S followed by steep declines of these responses in convalescent plasma.

• Download figure
• Open in new tab

Figure 3:

Distribution of Peak RU across various serological groups of SARS-CoV-2 infected individuals. Comparison of peak RU for (A) Spike(B) RBD and (C) Nucleocapsid between healthy (n=3), IgG-/IgM-(n=22), IgG+/IgM+ (n=11), IgG+/IgM-(n=14) and convalescent (n=22) groups. Red coloured symbols represent asymptomatic individuals. Pink coloured symbols indicate IgG-/IgM-individuals in convalescent group. Statistical significance was calculated by Mann-Whitney U-test; *, p < 0.05; **, p < 0.01; and ***, p<0.001

When grouped cross-sectionally by sero-status we observed (Figure 3), as would be expected, an increase in binding from the IgG-IgM-stage to the single positive, post class-switched stage (IgM-IgG+). Also, as observed in the smaller prospective cohort, association rates were higher for NC and RBD. For all three viral proteins, these rates were significantly lower in convalescent individuals with presumably mature memory B cells and viral clearance.

Having already demonstrated importance of circulating IL-6 and sMAdCAM either in combination (sMIL index, Figure 2) or in isolation (Figure 1 and Supplemental Figures 5 and 6) in influencing disease progression and anti-viral antibody production, we investigated the possibility of identifying a specific viral target (Spike, RBD or NC) whose binding response was associated and thus potentially attributable to these markers. Indeed, we were able to identify (Figure 4) a strong correlation between sMAdCAM levels and peak association rate of responses against the viral neutralization target, RBD in convalescent individuals. Also, we observed a relatively weaker but suggestive association with fold changes (day 0 and day 7) in sMIL index with that of viral Spike protein binding in plasma of the follow-up cohort.

• Download figure
• Open in new tab

Figure 4.

Association of sMAdCAM levels/sMIL index and virus specific binding in plasma(A) Association of plasma MAdCAM levels with Peak RU (association rate) of RBD binding among convalescent study participants. (B) Association of fold change between Peak RU for spike and sMIL index among study participants (n=9 out of 11) who were followed up. Fold change was calculated by dividing Day 7 values with those of Day 0. Fold change in Peak RU of 2 individuals could not be calculated as Day 0 values were undetectable (zero). Correlation analysis performed using non parametric Spearman Rank Correlation test. r indicates strength of correlation, p indicates significance.