Autoantibodies against chemokines post-SARS-CoV-2 infection correlate with disease course

Infection with severe acute respiratory syndrome coronavirus 2 associates with diverse symptoms, which can persist for months. While antiviral antibodies are protective, those targeting interferons and other immune factors are associated with adverse coronavirus disease 2019 (COVID-19) outcomes. Here we discovered that antibodies against specific chemokines were omnipresent post-COVID-19, were associated with favorable disease outcome and negatively correlated with the development of long COVID at 1 yr post-infection. Chemokine antibodies were also present in HIV-1 infection and autoimmune disorders, but they targeted different chemokines compared with COVID-19. Monoclonal antibodies derived from COVID-19 convalescents that bound to the chemokine N-loop impaired cell migration. Given the role of chemokines in orchestrating immune cell trafficking, naturally arising chemokine antibodies may modulate the inflammatory response and thus bear therapeutic potential.

Infection with severe acute respiratory syndrome coronavirus 2 associates with diverse symptoms, which can persist for months. While antiviral antibodies are protective, those targeting interferons and other immune factors are associated with adverse coronavirus disease 2019 (COVID- 19) outcomes. Here we discovered that antibodies against specific chemokines were omnipresent post-COVID-19, were associated with favorable disease outcome and negatively correlated with the development of long COVID at 1 yr post-infection. Chemokine antibodies were also present in HIV-1 infection and autoimmune disorders, but they targeted different chemokines compared with COVID-19. Monoclonal antibodies derived from COVID-19 convalescents that bound to the chemokine N-loop impaired cell migration. Given the role of chemokines in orchestrating immune cell trafficking, naturally arising chemokine antibodies may modulate the inflammatory response and thus bear therapeutic potential.
The spectrum of disease manifestations upon infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is broad and COVID-19 convalescent individuals often lament protracted symptoms over months, a condition referred to as long COVID or PASC (post-acute sequelae of COVID) 1-5 . Some evidence points to a role for immune dysregulation and autoimmunity as contributors to long COVID, although virus persistence has also been proposed 6-8 . Overall, however, there is little understanding of the biology underlying long COVID and the reasons behind the differences in COVID-19 manifestation.
Chemokines are chemotactic cytokines that mediate leukocyte trafficking and activity 9 . In addition to elevated levels of pro-inflammatory cytokines, acute COVID-19 is characterized by higher Letter https://doi.org/10.1038/s41590-023-01445-w The Milan cohort was sampled during the acute phase (day 8, on average) and at month 7 from disease onset (n = 44; 90.5% and 89.5% accuracy, respectively; Extended Data Fig. 2e and Supplementary Tables 2  and 4), and the Zurich cohort was evaluated at month 13 from disease onset (n = 104; 92.9% accuracy; Extended Data Fig. 2f and Supplementary Tables 2 and 4). Thus, COVID-19 was associated with a specific pattern of autoantibodies against chemokines.
To examine the relationship between chemokine antibodies and other serologic features of the Lugano cohort, we used ELISA and a pseudovirus-based neutralization assay to measure the binding and neutralizing capacity of antibodies against SARS-CoV-2 (ref. 27 ). In agreement with previous studies 27 , IgG binding to the Spike receptor binding domain (RBD) of SARS-CoV-2 and plasma half-maximal SARS-CoV-2 neutralizing titers (NT 50 ) were variable (Fig. 1a), but positively correlated with each other and with age (Extended Data Fig. 3a-c) 27 . In contrast, there was no correlation of NT 50 or RBD IgGs with the levels of COVID-19 signature chemokine antibodies (against CCL19, CCL22 and CXCL17), or with the sum of all chemokine IgG reactivities (Extended Data Fig. 3d). A weak negative correlation between age and the sum of all chemokine IgG reactivities was observed (Extended Data Fig. 3d), but there were no significant differences in the levels of the COVID-19 signature chemokine antibodies between males and females (Extended Data Fig. 3e). These observations indicated that, after COVID-19, autoantibodies against specific chemokines were not correlated with the antibodies against SARS-CoV-2.
To document the temporal evolution of chemokine anti bodies following COVID-19, we compared the reactivities of plasma collected from the Lugano cohort at around month 6 and month 12 post-symptom onset side-by-side (Extended Data Fig. 4a and Supplementary Table 2). In agreement with earlier findings 28 , RBD antibodies significantly decreased in unvaccinated COVID-19 convalescents, while they increased in those receiving at least one dose of messenger RNA-based COVID-19 vaccine (Extended Data Fig. 4b and Supplementary Table 2). Conversely, and regardless of vaccination status, antibodies to CCL19, CCL8, CCL13, CCL16, CXCL7 and CX3CL1 significantly increased, those to CXCL17 remained generally stable and those to CCL22 followed variable kinetics at month 12 compared with month 6 in the COVID-19 convalescents ( Fig. 1d and Extended Data Fig. 4c). To further investigate the kinetics of COVID-19 signature chemokine anti bodies, we analyzed COVID-19 convalescents from the Lugano cohort for which acute samples were also available (n = 12; Extended Data Fig. 4d and Supplementary Table 2). During acute COVID-19, IgG antibodies to CCL19, but not to CCL22 or CXCL17, were already higher than in healthy controls, and continued to increase until month 12 (Extended Data Fig. 4e). Similarly, in the Milan cohort, autoantibodies against CCL19, CCL22 and CXCL17 were higher during the acute phase than in healthy controls (Extended Data Fig. 2e). In contrast to natural infection, no significant change in antibody reactivity to any of the chemokines expression of certain chemokines [10][11][12][13][14][15][16] . Accordingly, chemokines recruit neutrophils and monocytes to sites of infection, where they play a key role in the pathophysiology of COVID-19 by sustaining inflammation and causing collateral tissue damage and fibrosis 11,15,17,18 . Autoantibodies to type I interferon (IFN) and other molecules are observed in COVID- 19 and are generally associated with adverse outcome [19][20][21][22][23][24][25] . However, no study comprehensively investigated chemokine antibodies, nor their persistence over time.
Here we devised a peptide-based strategy to measure antibodies that bind to a functional region of each of the 43 human chemokines. By examining three independent COVID-19 cohorts, we found that the presence of autoantibodies against specific chemokines helped identify convalescent individuals with favorable acute and long COVID disease course. Monoclonal antibodies targeting chemokines that were derived from these individuals blocked leukocyte migration and thus may be beneficial through modulation of the inflammatory response.
To evaluate chemokine antibodies after COVID-19, we obtained plasma samples from a cohort of COVID-19 convalescents (n = 71) at month 6 (on average) post disease onset (hereafter, the Lugano cohort; Supplementary Tables 1 and 2). Plasma samples from COVID-19 uninfected individuals, confirmed by negative serologic test, were used as healthy controls (n = 23). Because the N-terminal loop (N-loop) of chemokines is required for receptor binding, we investigated whether biologically active chemokine antibodies targeted this region, whose sequence is specific for each but two of the 43 human chemokines (Extended Data Fig. 1a) 26 . We designed peptides corresponding to the N-loop of each chemokine for use in ELISA assays (Supplementary  Table 3), measured the amount of peptide-specific IgG antibody in serial plasma dilutions and plotted the signal as a heatmap ( Fig. 1a and Extended Data Fig. 1b). Analysis of all parameters by nonlinear dimensionality reduction with t-distributed stochastic neighbor embedding (t-SNE) indicated a clear separation between healthy controls and COVID-19 convalescents (Extended Data Fig. 2a). Some COVID-19 convalescent plasma indicated high levels of IgGs to certain chemokines (for example, CCL8, CXCL13 and CXCL16) compared with healthy controls. For these chemokines, antibody levels to the N-loop significantly correlated with those against the C-terminal region of the same chemokine, suggesting that, when present, antibodies formed against multiple chemokine epitopes ( Fig. 1a and Extended Data Fig. 2b). When considering antibodies against each chemokine individually, a significant difference in reactivity, shown as ratio over healthy controls, was observed for peptides corresponding to 23 of the 43 chemokines (Extended Data Fig. 2c). Antibodies to the three chemokines with P < 10 −4 (CCL19, CCL22 and CXCL17) clustered together and by themselves were sufficient to correctly assign healthy controls and COVID-19 convalescents with high accuracy (96.8%; Fig. 1a-c and Extended Data Fig. 2c,d), so they were defined as a 'COVID-19 signature'. This signature was validated in two independent COVID-19 cohorts.  (COVID-19) of the Lugano cohort at month 6. Samples are ranked according to the level of SARS-CoV-2 RBD reactivity. Chemokine IgGs are ordered by unsupervised clustering analysis of ELISA signal. SARS-CoV-2 pseudovirus neutralizing activity (NT 50 ) and IgG binding to peptides corresponding to negative control, IFN-α2 and SARS-CoV-2 nucleocapsid protein (N) are shown. b, AUC of ELISA showing IgG antibodies to CCL19, CCL22 and CXCL17 (COVID-19 signature) in healthy controls and COVID-19 convalescents at month 6. Two-tailed Mann-Whitney U-test. c, Logistic regression analysis showing the assignment of COVID-19 convalescents and healthy controls based on CCL19, CCL22 and CXCL17 antibodies at month 6. d, AUC of ELISA showing CCL19, CCL22 and CXCL17 antibodies at months 6 and 12 in COVID-19 convalescents (n = 63). Wilcoxon two-tailed signed-rank test. e, Chemokine antibodies in previously hospitalized and outpatient COVID-19 convalescents at month 6, shown as ratio over healthy controls. Circle size indicates significance; colors show the log 2 fold-change increase (red) or decrease (blue), shown as ratio over healthy controls. Kruskal-Wallis test followed by Dunn's multiple comparison test. f, Cumulative AUC of ELISA signal of the IgGs against the 42 chemokine N-loops in healthy controls and previously hospitalized and outpatient COVID-19 convalescents at month 6. Kruskal-Wallis test followed by Dunn's multiple comparison test. g, t-SNE distribution of previously hospitalized and outpatient COVID-19 convalescents at month 6, as determined with the 42 datasets combined. h, Logistic regression analysis showing the assignment of previously hospitalized and outpatient COVID-19 convalescents based on CXCL5, CXCL8 and CCL25 antibodies (COVID-19 hospitalization signature) at month 6. In b and f, horizontal bars indicate median values. In a-h, AUC values are the average from two independent experiments. Healthy controls (n = 23) in a, b, c, e and f; COVID-19 convalescents (n = 71) in a, b and c, of which previously hospitalized (n = 50) and outpatient (n = 21) in e-h. Ab, antibody; AUC, area under the curve; FC, fold-change; ID, identity; m, months; NS, not significant.  Table 2). Thus, unlike the antibodies to SARS-CoV-2 RBD, which decreased over time, the levels of some chemokine antibodies that were present during acute COVID-19 increased over 1 yr of observation.  CCL7  CXCL13  CCL13  CXCL1  CCL27  CXCL2  XCL1/XCL2  CCL26  CCL2  CCL25  CCL21  CXCL14  CXCL6  CCL5  CCL28  CCL11  CXCL11  CCL20  CCL17  CXCL8  CXCL5  CXCL10  CXCL9  CCL3  CCL15  CCL4  CCL24  CCL23  CXCL4  CXCL3  CCL18  CCL14  CCL1  CX3CL1  CXCL7  CXCL12  CCL8  CXCL16   COVID-   The presence of autoantibodies in a portion of hospitalized patients with COVID-19 is linked to severe illness 20,21,24,25 . No significant difference in age distribution was observed between the previously hospitalized (n = 50) and outpatient (n = 21) COVID-19 convalescents in the Lugano cohort (60 ± 14 and 57 ± 15 yr, respectively), while a higher proportion of males was observed in the formerly hospitalized, but not outpatient, group (60% and 38.1%, respectively; Supplementary Table 2). When the most significant differences in autoantibody levels at month 6 were considered (P < 10 −4 ), only the antibodies against CCL19 were higher in previously hospitalized COVID-19 convalescents compared with healthy controls, while antibodies against 8 chemokines (CXCL8, CCL22, CXCL16, CCL27, CXCL7, CCL20, CX3CL1 and CCL19) were increased in outpatient COVID-19 convalescents compared with healthy controls (Fig. 1e and Extended Data Fig. 5a). Only the outpatient COVID-19 convalescents displayed significantly higher cumulative IgG reactivity against the chemokines compared with healthy controls (Fig. 1f).
Thus, a broader pattern and higher overall amounts of autoantibodies against chemokines were observed at month 6 in COVID-19 convalescents who were outpatients during the acute phase of the disease.
Direct comparison by t-SNE analysis of all chemokine antibodies at month 6 separated COVID-19 convalescents in the Lugano cohort that were previously hospitalized from outpatients (Fig. 1g). Antibodies against CXCL5, CXCL8 and CCL25 were all lower in previously hospitalized compared with outpatient COVID-19 convalescents (Extended Data Fig. 5a,b), and this was not linked to the therapy received during hospitalization (Supplementary Table 2). The combination of antibody values against CXCL5, CXCL8 and CCL25 at month 6 alone could correctly assign formerly hospitalized and outpatient COVID-19 convalescents with an accuracy of 77.5% ('COVID-19 hospitalization signature';   In the COVID-19 hospitalization signature, antibodies to CXCL5 and CXCL8 were negatively correlated with RBD IgG and age, but not with sex (Extended Data Fig. 6a,b). Consistent with previous work 27 , both RBD IgG and NT 50 values were significantly higher in previously hospitalized compared with outpatient COVID-19 convalescents and in males compared with females of both hospitalized and outpatient groups (Extended Data Fig. 6c) 27 . Thus, the chemokine antibody signature that distinguished healthy controls from COVID-19 convalescents (autoantibodies to CCL19, CCL22 and CXCL17) was different from the signature associated with COVID-19 severity (autoantibodies to CXCL5, CXCL8 and CCL25). A fraction of individuals who recover from COVID-19 experience long-term sequelae 1,3,4 . To determine whether a specific pattern of chemokine antibodies at month 6 was predictive of long COVID, we used a questionnaire to collect information on self-reported symptoms (for example, pulmonary, systemic, neurological and psychiatric) at month 12 from the Lugano cohort (Supplementary Table 5). 65.1% of all participants reported the persistence of at least one symptom related to COVID-19. Among these, the average number of long-term symptoms was 3.3, and they were more frequent in formerly hospitalized (72.7%) compared with outpatient (47.4%) COVID-19 convalescents (Extended Data Fig. 7a and Supplementary Tables 2 and 5). No differences in age or sex distribution or time from disease onset to month 12 visit were observed between COVID-19 convalescents with protracted symptoms at month 12 (long COVID) or without (no long COVID) (Extended Data Fig. 7b).
COVID-19 convalescents in the Lugano cohort with long COVID, particularly outpatients and females, showed significantly lower cumulative levels of chemokine antibodies compared with those with no long COVID ( Fig. 2a and Extended Data Fig. 7c,d). In contrast, RBD IgG and NT 50 values were comparable between long COVID and no long COVID groups (Fig. 2a). The cumulative amount of chemokine antibodies did not correlate with the number of symptoms (Extended Data Fig. 7e). IgG antibodies against CCL21, CXCL13 and CXCL16 at month 6 distinguished long COVID from no long COVID groups with high significance and were defined as a 'long COVID signature' (Fig. 2b and Extended Data Fig. 7f). Logistic regression analysis using the antibody values for CCL21, CXCL13 and CXCL16 alone predicted the absence of persistent symptoms with 77.8% accuracy (Fig. 2c). Similarly, analysis of the Zurich cohort at month 13 showed 72.1% accuracy of association with lack of long COVID, even though in that cohort only CCL21 antibodies were significantly different between long COVID and no long COVID groups (Extended Data Fig. 7g). These results indicated that specific patterns of chemokine antibodies at month 6 were associated with the longer-term persistence of symptoms after COVID-19.
Because chemokine antibodies to CXCL13 and CXCL16 were associated with decreased likelihood of long COVID, we next derived corresponding memory B cell antibodies from available peripheral blood mononuclear cell (PBMC) samples (Supplementary Table 6). Three N-loop binding monoclonal antibodies were obtained for CXCL16, which blocked migration of the 300.19 preB cell line expressing the CXCL16 cognate receptor CXCR6 (Fig. 2d,e, Extended Data Fig. 8a,b and Supplementary Tables 6 and 7). Similarly, three CXCL13 N-loop antibodies bound in ELISA and inhibited chemotaxis of primary CD19 + human B cells (Fig. 2f,g, Extended Data Fig. 8c and Supplementary Tables 6 and 7). By the same approach, we discovered chemotaxis-blocking antibodies specific for CCL8 and CCL20 (Extended Data Fig. 8d-j and Supplementary Tables 6 and 7). Consistent with these results, polyclonal plasma IgG from COVID-19 convalescents effectively blocked chemotaxis at concentrations 50 times lower than those found in human serum 29 (Extended Data Fig. 8k). Thus, antibodies from COVID-19 convalescents that bound to the N-loop of chemokines were biologically active.
We observed no correlation between the amounts of chemokines and the amounts of corresponding autoantibodies in the acute phase or at month 7 post-infection in the Milan cohort (Fig. 3c). We also did not detect an increase in CCL3 and CCL4 antibodies in COVID-19 convalescents compared with healthy controls in the Lugano cohort at month 6, although the amount of corresponding chemokines was elevated in their plasma ( Fig. 3a and Extended Data Fig. 5a). As such, even though various chemokines rapidly increased and persisted in plasma for at least 6 months post-infection, their level did not correlate with the amount of corresponding autoantibodies in the circulation.
Several chemokines are detected in high amounts in bronchoalveolar and other fluids during COVID-19, fueling a pro-inflammatory environment in the lungs, which likely contributes to critical illness and A c u t e  Letter https://doi.org/10.1038/s41590-023-01445-w hospitalization 10-14 . We found autoantibodies against CXCL5, CXCL8 and CCL25 in COVID-19, but there was no correlation with the amount of the corresponding chemokines in plasma. Because these chemokines attract neutrophils and other cell types that promote inflammation and tissue remodeling, the presence of the corresponding autoantibodies suggests protection through dampening of the damaging inflammatory response associated with severe COVID-19. Autoantibodies to CCL21, CXCL13 and CXCL16 were increased in recovered individuals compared with those with long COVID 1 yr post-infection. These chemokines are important for tissue trafficking and activation of T and B lymphocytes. It is possible that their respective autoantibodies positively impact the long-term outcome of COVID-19 by antagonizing or otherwise modulating the activation, recruitment and retention of these cell types. Persistent immune responses have been proposed as a mechanism for long COVID, and chemokines have been implicated in its pathogenesis 1,33 .
Infection can trigger antibody polyreactivity and autoimmunity which are generally deleterious 34-36 . Because chemokine antibodies are present in plasma after COVID-19 at concentrations able to impair cellular migration, the variety and amount of chemokine antibodies that are present or induced upon infection in each individual may positively modulate the quality and strength of the inflammatory response, which in turn would impact disease manifestation, severity and long COVID. Further studies are needed to determine whether agents that target the chemokine system could impact positively on the early inflammatory phase of COVID-19 and reduce the development of long COVID.

Online content
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41590-023-01445-w.

Material availability
Material used in the present study is available upon request from the lead contact and may require a Material Transfer Agreement (MTA).
A key resources table is provided as Supplementary Table 8.

Study participants and ethical approvals
The Lugano COVID-19 cohort included 71 participants, diagnosed with COVID-19 at the Clinica Luganese Moncucco (CLM, Switzerland) between 8 March 2020 and 22 November 2020, who were enrolled in the study and divided into two groups, according to the severity of the acute disease. The hospitalized group included 50 participants; the outpatient group included 21 close contacts of the hospitalized group, who received only at-home care. Inclusion criteria for the hospitalized group were a SARS-CoV-2-positive nasopharyngeal swab test by quantitative PCR with reverse transcription (RT-qPCR) and age ≥ 18 yr. Inclusion criteria for the outpatient group were being a symptomatic close contact (living in the same household) of an individual enrolled in the hospitalized group and age ≥ 18 yr. Serologic tests confirmed COVID-19 positivity for all the participants ( Fig. 1a and Extended Data Fig. 3a). At the 12-month visit, participants were asked to indicate the presence or absence of persisting symptoms related to COVID-19 according to a questionnaire (Supplementary Table 5 The Zurich COVID-19 cohort 7 included 104 participants, diagnosed with COVID-19 at the University Hospital Zurich, the City Hospital Triemli Zurich, the Limmattal Hospital or the Uster Hospital between April 2020 and April 2021, who were included in the study and divided into two groups, according to the severity of the acute disease. The hospitalized group included 38 participants, whereas the outpatient group included 66 individuals, who received only at-home care. Inclusion criteria for the participants were a SARS-CoV-2-positive nasopharyngeal swab test by RT-qPCR and age ≥ 18 yr. At the 13-month visit, blood was collected and participants were asked by trained study physicians to indicate the presence or absence of persisting symptoms related to COVID-19. Patients who reported at least one symptom at month 13 were included in the long COVID group. The study was performed in compliance with all relevant ethical regulations and the study protocols were approved by the Cantonal Ethics Committee of Zurich (Business Administration System for Ethics Committees (BASEC) no. 2016-01440).
The healthy control cohort included 15 adult participants (≥18 yr) with self-reported absence of previous SARS-CoV-2 infection or vaccination, enrolled between November 2020 and June 2021. An additional eight pre-pandemic samples were obtained from blood bank donors (ECCT: CE-3428). Serologic tests confirmed COVID-19 negativity for all healthy controls (Fig. 1a and Extended Data Fig. 4a). These samples were used as controls for the COVID-19 convalescents (Lugano, Milano and Zurich) and HIV-1, AS, RA and SjS cohorts.
The vaccination cohort included 16 adult participants (≥18 yr) with self-reported absence of previous SARS-CoV-2 infection (confirmed by negative serologic test; Extended Data Fig. 4f) and who received two doses of mRNA-based COVID-19 vaccine 37,38 , enrolled on the day of first vaccine dose or earlier, between November 2020 and June 2021 (ECCT: CE-3428).
The HIV-1 and autoimmune diseases cohorts included pre-pandemic plasma samples obtained from the following participants: 24 HIV-1 positive (ECCT: CE-813) 39  Written, informed consent was obtained from all participants, and all samples were coded. No compensation was provided to the study participants. Demographic and clinical features of the cohorts are reported in Supplementary Table 2.

Blood collection, processing and storage
Blood was collected by venipuncture at approximately 6-month intervals and the PBMCs were isolated using Histopaque density centrifugation (Lugano and healthy control cohorts). Total PBMCs were aliquoted and frozen in liquid nitrogen in the presence of FCS and dimethylsulfoxide. Plasma was aliquoted and stored at −20 °C or less. Before use, plasma aliquots were heat-inactivated (56 °C for 1 h) and then stored at 4 °C. For chemotaxis assays, CD14 + monocytes and CD19 + B cells were enriched from fresh PBMCs derived from blood donors (Swiss Red Cross Laboratory; ECCT: CE-3428) through positive immunoselection (130-050-201 and 130-050-301, respectively, Miltenyi Biotec) according to the manufacturer's instructions. After isolation, CD19 + B cells were rested overnight in RPMI-1640 medium supplemented with 10% (v/v) FBS, 1% (v/v) nonessential amino acids, 1 mM sodium pyruvate, 2 mM GlutaMAX, 50 μM β-mercaptoethanol and 50 U ml −1 penicillin/ streptomycin (all from Gibco) before being used in chemotactic assays. For the other cohorts, see refs. 7,40 .

Reagents
Peptides. Synthetic peptides containing the N-loop or the C-terminal sequence of human chemokines were designed and obtained (>75% purity) from GenScript. All peptides are biotinylated (biotin-Ahx) at the N terminus and amidated at the C terminus. In addition, the first 2-4 amino acids of each peptide (GS, GGS, GGGS or GGK depending on the length of the N-loop/C terminus of the chemokine) consist of a linker between the biotin and the chemokine sequence. Peptides are generally 25 amino acids long, or 22-25 amino acids when synthesis was problematic. The sequence of the IFN-α2 peptide (7-28) was based on a previously described immunoreactive epitope in patients with myasthenia gravis 41 , and that of the SARS-CoV-2 nucleocapsid protein (N) peptide (157-178) was described in ref. 42 . An irrelevant peptide was used as negative control. The amino acid sequences of all peptides in this study are listed in Supplementary Table 3.

Chemotaxis
The migration of primary human monocytes and B cells isolated from buffy coats or of murine preB 300.19 cells stably expressing the human chemokine receptors CCR2 (ref. 46 ), CCR6, CXCR1 (ref. 47 ) and CXCR6 (ref. 48 ) was assayed using 48-well Boyden chambers (Neuro Probe) with polyvinylpyrrolidone-free polycarbonate membranes with pore size of 3 μm for primary human B cells and 5 μm for the other cell types, as previously described 49 . Briefly, 10 5 primary human B cells or 5 × 10 4 primary human monocytes and murine preB 300.19 cells were diluted in RPMI-1640 supplemented with 20 mM HEPES, pH 7.4, and 1% pasteurized plasma protein solution (5% PPL SRK; Swiss Red Cross Laboratory). Cells were then added to the upper wells and the chemokine (with or without antibodies) to the bottom wells. After 120 min of incubation for primary human B cells and 90 min for the other cell types, the membrane was removed, washed on the upper side with PBS, fixed and stained with DiffQuik. All assays were done in triplicate, and for each well the migrated cells were counted at 100-fold magnification in five randomly selected high-power fields.

Inhibition of chemotaxis by monoclonal antibodies.
Experiments were performed with monoclonal antibodies at a final concentration of 30 μg ml −1 (Extended Data Fig. 8g) or 50 μg ml −1 (Fig. 2e,g and Extended Data Fig. 8j). Baseline migration was determined in the absence of chemoattractant (buffer control).

ELISA
To evaluate the antibodies' binding to chemokine peptides, 96-well plates (ThermoFisher, 442404) or 384-well plates (ThermoFisher, 464718) were coated with 50 μl (or 10 μl for 384-well plates) per well of a 2 μg ml −1 Neutravidin (Life Technologies, 31000) solution in PBS, overnight at room temperature. Plates were washed four times with washing buffer (PBS + 0.05% Tween-20 (Sigma-Aldrich)) and incubated with 50 μl (or 10 μl for 384-well plates) per well of a 50 nM biotinylated peptide solution in PBS for 1 h at room temperature. After washing four times with washing buffer, plates were incubated with 200 μl (or 50 μl for 384-well plates) per well of blocking buffer (PBS + 2% BSA + 0.05% Tween-20) for 2 h at room temperature. Plates were then washed four times with washing buffer, and serial dilutions of monoclonal antibodies or plasma were added in PBS + 0.05% Tween-20 and incubated for 1 h at room temperature. To screen for the presence of chemokine IgGs, plasma samples were assayed (unless otherwise stated) at 1:50 starting dilution followed by three fourfold serial dilutions (1:200, 1:800, 1:3,200). Monoclonal antibodies were tested at 5 μg ml −1 starting concentration followed by 11 threefold serial dilutions. Plates were subsequently washed four times with washing buffer and incubated with an antibody against human IgG secondary antibody conjugated to horseradish peroxidase (HRP) (GE Healthcare, NA933) at a 1:5,000 dilution in PBS + 0.05% Tween-20. Finally, after washing four times with washing buffer, plates were developed by the addition of 50 μl (or 10 μl for 384-well plates) per well of the HRP substrate tetramethylbenzidine (TMB) (ThermoFisher, 34021) for 10 min. The developing reaction was stopped with 50 μl (or 10 μl for 384-well plates) per well of a 1 M H 2 SO 4 solution, and absorbance was measured at 450 nm with an ELISA microplate reader (BioTek) with Gen5 3.12 software. A positive control (broadly reactive plasma from donor CLM70) and negative control (uninfected participant) samples were included in each experiment. Since the basal average optical density likely also depends on intrinsic features of each peptide that is used to coat the ELISA plate, the presented values should be interpreted as relative rather than absolute. The area under the curve (AUC) was obtained from two independent experiments and plotted with GraphPad Prism v.9.0.2. The main findings were further confirmed by assaying subsets of samples belonging to the different groups, side-by-side on the same plates.
Reactivity at month 6 versus 12. Experiments were performed with plasma samples from different time points side-by-side on the same plate. In Extended Data Fig. 4b, plasma was assayed at a 1:50 starting dilution, followed by four additional fivefold dilutions. RBD IgG antibody levels were measured in COVID-19 convalescents who had not received a COVID-19 mRNA vaccine between first and second visits (no vaccination) or in individuals with at least one dose of vaccine at least 10 d before blood sampling at the second visit (Fig. 1d, Extended Data  Fig. 4b,c and Supplementary Table 2).
Kinetics of signature chemokine IgG antibodies. Experiments were performed with plasma samples from different time points assayed at 1:50 dilution side-by-side on the same plate, and the average optical density at 450 nm obtained from two independent experiments was plotted with GraphPad Prism v.9.0.2 (Extended Data Fig. 4e).

IgG antibodies binding to SARS-CoV-2 RBD.
Experiments were performed with 96-well plates coated with 50 μl per well of a 5 μg ml −1 protein solution in PBS overnight at room temperature, and subsequently blocked and treated as described above. In this case, plasma samples were assayed at a 1:50 starting dilution either followed by seven additional threefold serial dilutions (Figs. 1a and 2a and Extended Data Figs. 3a,c,d and 6a) or followed by three additional fivefold serial dilutions (Extended Data Fig. 4b,f).

Antibody sequencing, cloning, production and purification
Antibody genes were sequenced, cloned and expressed as previously reported [50][51][52] . Briefly, reverse transcription of RNA from FACS-sorted single cells was performed to obtain complementary DNA, which was then used for amplification of the immunoglobulin IGH, IGK and IGL genes by nested PCR. Amplicons from this first PCR reaction served as templates for sequence-and ligation-independent cloning into human IgG1 antibody expression vectors. Monoclonal antibodies were produced by transiently transfecting Expi293F cells cultured in Freestyle-293 Expression Medium (ThermoFisher) with equal amounts of immunoglobulin heavy and light chain expression vectors using polyethyleneimine Max (PEI-MAX, Polysciences) as a transfection reagent. After 6-7 d of culture, cell supernatants were filtered through 0.22-μm Millex-GP filters (Merck Millipore), and antibodies were purified using Protein G Sepharose 4 Fast Flow (Cytiva) according to the manufacturer's instructions and buffer-exchanged and concentrated in PBS by Amicon Ultra-4 centrifugal filters (30-kDa cutoff, Millipore). Where indicated, the monoclonal antibody against Zika virus Z02150 was used as an isotype control.

Computational analysis of antibody sequences
Antibody sequences were analyzed using a collection of Perl and R scripts provided by IgPipeline and publicly available on GitHub (https:// github.com/stratust/igpipeline) 27 . In brief, sequences were annotated using IgBlast 53 v.1.14.0 with IMGT domain delineation system and the Change-O toolkit v.0.4.5 (ref. 54 ). CDR3 sequences were determined by aligning the IGHV and IGLV nucleotide sequences against their closest germlines using the blastn function of IgBlast.

SARS-CoV-2 pseudotyped reporter virus and neutralization assay
To generate (HIV-1/NanoLuc2AEGFP)-SARS-CoV-2 particles, HEK293T cells were co-transfected with the three plasmids pHIV NL Gag-Pol, pCCNanoLuc2AEGFP and SARS-CoV-2 S as described elsewhere 27,55 . Supernatants containing virions were collected 48 h after transfection, and virion infectivity was determined by titration on 293T ACE2 cells. The plasma neutralizing activity was measured as previously reported 27,55 . Briefly, threefold serially diluted plasma samples (from 1:50 to 1:328,050) were incubated with SARS-CoV-2 pseudotyped virus for 1 h at 37 °C, and the virus-plasma mixture was subsequently incubated with 293T ACE2 cells for 48 h. Cells were then washed with PBS and lysed with Luciferase Cell Culture Lysis 5X reagent (Promega). Nanoluc Luciferase activity in cell lysates was measured using the Nano-Glo Luciferase Assay System (Promega) with Modulus II Microplate Reader user interface (TURNER BioSystems). The obtained relative luminescence units were normalized to those derived from cells infected with SARS-CoV-2 pseudotyped virus in the absence of plasma. The NT 50 values were determined using four-parameter nonlinear regression with bottom and top constrains equal to 0 and 1, respectively (Graph-Pad Prism v.9.0.2). The dotted line (NT 50 = 5) in the plots represents the lower limit of detection of the assay.

Model interaction between chemokine and chemokine receptor
The illustrative model in Extended Data Fig. 1a was generated from the structure of inactive CCR2 (PDB code: 5T1A) 56 , together with the electron microscopy structures of CCR5 and CCR6 (PDB codes: 6MEO and 6WWZ, respectively) 57,58 , by using the SWISS-MODEL 59 server and the molecular graphics program PyMOL 2.5.0 for modeling the N and C termini of the receptor. The crystal structure of CCL8 (MCP-2) (PDB code: 1ESR) 60 and the electron microscopy structure of CCR6 (ref. 58 ) were used to model the complex. The intracellular residues were removed for clarity.

Statistical analysis
Sample size definition. No statistical methods were used to pre-determine sample sizes but our sample sizes are similar to those reported in previous publications 7 .
Tests for statistical significance. Upon testing of parametric assumptions (Kolmogorov-Smirnov test for normality and Hartley's F max test for homoskedasticity), statistical significance between two groups was determined using parametric paired two-tailed Student's t-test, or nonparametric two-tailed Mann-Whitney U-tests (unpaired samples), or Wilcoxon signed-rank test (paired samples). Statistical significance between more than two groups was evaluated using Kruskal-Wallis test (followed by Dunn multiple comparisons), one-way analysis of variance (ANOVA) (followed by Tukey multiple comparisons) or two-way repeated measures ANOVA (followed by Šídák multiple comparisons), as described in the figure legends. Statistical significance of the signature chemokines (CCL19, CCL22, CXCL17, CXCL8, CCL25, CXCL5, CCL21, CXCL13 and CXCL16) was also confirmed when applying the Bonferroni criterion to guarantee a familywise level of significance equal to 0.05. Statistical significance from a 2 × 2 contingency table was determined with Fisher's exact test. Correlations were assessed using Pearson correlation analysis. A P value of less than 0.05 was considered statistically significant. Data and statistical analyses were performed with GraphPad Prism v.9.0.2. Data collection and analysis were not performed blind to the conditions of the experiments.

t-SNE. t-SNE
analysis was performed using the Rtsne R package v.0.15 (https://CRAN.R-project.org/package=Rtsne) using the AUC values for all chemokines. The theta parameter for the accuracy of the mapping was set to zero in all cases for exact t-SNE.
Clustering. Hierarchical clustering was created using the hclust R function v.4.1.1. Clustering analysis was performed using the correlation as distance and Ward's method for agglomeration. Heatmaps were created with either GraphPad Prism v.9.0.2 ( Fig. 1a and Extended Data Fig.  2d) or the Pretty Heatmaps (pheatmap) R package v.1.0.12 (Extended Data Fig. 10b). In Extended Data Fig. 10b, each column containing a distinct chemokine was scaled with the scaling function provided by R, which sets the mean and the standard deviation to 0 and 1, respectively.

Logistic regression and additional analyses.
Logistic regression was performed using the GLM (Generalized Linear Models) function provided by the R package v.4.1.1. To identify which variables to include in the analysis, AUCs were ranked according to the P value obtained with a Mann-Whitney-Wilcoxon nonparametric test on the Lugano cohort.
Letter https://doi.org/10.1038/s41590-023-01445-w The first N variables minimizing the Akaike information criterion were then used in the fitting. Furthermore, the same set of variables was used to perform the fitting with the Milan and Zurich cohorts. In each plot, values from 0 to 0.5 and from 0.5 to 1 on the y axis represent the assignment of individuals to the A and B groups (of a Prediction A versus B; see gray backgrounds), respectively. On the x axis, samples are divided into the two groups and subsequently ordered according to sample identity as shown in Supplementary Table 2. Dots in the gray area represent individuals that are assigned to the correct group. We additionally performed chi-tests considering covariates that are known to influence COVID-19 severity (demographics (sex and age) and comorbidities (diabetes and cardiovascular diseases)) and found that none of them was significantly different between groups. Race/ ethnicity was not analyzed because the cohort is nearly 100% White; similarly, immune deficiency was rare and was not considered. Logistic regression analysis using the combination of these covariates (age, sex, diabetes and cardiovascular diseases) allowed proper assignment with accuracies of 74.6% (COVID-19 severity; outpatient versus hospitalized) and 68.3% (long COVID; no long COVID versus long COVID). Notably, the accuracy using chemokine antibody values is even better (77.5% (COVID-19 severity) and 77.8% (long COVID)). These analyses are shown in Supplementary Table 9.

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability
All data analyzed during the present study are included in this article and its supporting information files. Source data are provided with this paper. PDB accession codes are 5T1A, 6MEO, 6WWZ and 1ESR.