Abstract

Background

Excessive complement activation has been implicated in the pathogenesis of coronavirus disease 2019 (COVID-19), but the mechanisms leading to this response remain unclear.

Methods

We measured plasma levels of key complement markers, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA and antibodies against SARS-CoV-2 and seasonal human common cold coronaviruses (CCCs) in hospitalized patients with COVID-19 of moderate (n = 18) and critical severity (n = 37) and in healthy controls (n = 10).

Results

We confirmed that complement activation is systemically increased in patients with COVID-19 and is associated with a worse disease outcome. We showed that plasma levels of C1q and circulating immune complexes were markedly increased in patients with severe COVID-19 and correlated with higher immunoglobulin (Ig) G titers, greater complement activation, and higher disease severity score. Additional analyses showed that the classical pathway was the main arm responsible for augmented complement activation in severe patients. In addition, we demonstrated that a rapid IgG response to SARS-CoV-2 and an anamnestic IgG response to the nucleoprotein of the CCCs were strongly correlated with circulating immune complex levels, complement activation, and disease severity.

Conclusions

These findings indicate that early, nonneutralizing IgG responses may play a key role in complement overactivation in severe COVID-19. Our work underscores the urgent need to develop therapeutic strategies to modify complement overactivation in patients with COVID-19.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first detected in December 2019 and quickly spread throughout the globe, causing the coronavirus disease 19 (COVID-19) pandemic. Severe COVID-19 has been correlated with activation of the complement system, as indicated by the presence of low levels of C3 and high levels of the anaphylatoxin C3a, a cleavage product of C3 [1–9]. Recent studies have found deposits of complement factors (C1q, C4, C3, soluble C5-9 [sC5-9], and factor B) in lung tissue [10, 11] and in the kidneys [12] of patients with severe COVID-19. In addition, overactivation of the complement system has been strongly associated with higher mortality risk [4, 7, 8, 13]. Nevertheless, the mechanisms driving overactivation of the complement system in severe COVID-19 are poorly understood.

The complement system comprises a complex network of proteins, forming a highly regulated cascade of reactions. The system is activated through 3 pathways: the classical, the lectin, and the alternative pathway. Activation of the classical pathway is mediated by the binding of C1q to immune complexes or by binding directly to pathogen surfaces. It has been hypothesized that excessive immune complexes could be the cause of overactivation of the complement system leading to severe COVID-19 [14], but this has not yet been tested. It is well established that severe COVID-19 is correlated with faster and more intense binding and neutralizing antibody response; however, the nature and role of these antibodies in activating complement and causing severe disease have not been characterized. Moreover, several reports indicate that a significant fraction of plasmablasts, antibodies and CD4+ T cells generated in response to SARS-CoV-2 infection are cross-reactive with the seasonal human common cold coronaviruses (CCCs) [15–20], but the impact of these anamnestic immune response on protection and disease is not well understood.

Here we report a comprehensive analysis of the complement and antibody response to SARS-CoV-2 infection in a cohort of moderately and critically ill patients hospitalized with COVID-19. We demonstrate that systemic complement activation is highly correlated with disease severity. Our findings reveal that overactivation of the complement system is mediated by the classical pathway in response to increased levels of circulating immunoglobulin (Ig) G immune complexes. Antibody responses to SARS-CoV-2 develop faster and at greater levels in patients with more severe COVID-19. Importantly, we show that antibody responses to the nucleoprotein (NP) of CCCs are boosted in patients with more severe COVID-19 and that high levels of IgG anti-NP against CCCs are associated with complement activation via the classical pathway.

METHODS

Patients and Clinical Samples

The study design and data collection for the patients with COVID-19 included in this study have been described elsewhere [21, 22]. Briefly, we enrolled patients with COVID-19 hospitalized at the University of Pittsburgh Medical Center Presbyterian and Shadyside hospitals (April–September 2020) in an observational study. We included 55 adult patients (≥18 years old) who were hospitalized either in an intensive care unit (ICU) or in dedicated hospital wards (non-ICU) for COVID-19. SARS-CoV-2 infections were documented by positive quantitative polymerase chain reaction (PCR) test results, using nasopharyngeal swab samples. Demographic and clinical variables and plasma samples were collected at day 1 after enrollment for all patients. Prepandemic plasma samples (n = 10) were obtained from healthy adults and included as controls (healthy controls [HCs]) for the complement protein measurements. We also included samples from patients enrolled in the ongoing University of Pittsburgh Acute Lung Injury Registry study [23] (n = 15; Supplementary Table 1) for comparisons of CCC antibody response. These patients were hospitalized in the ICU and met the diagnostic criteria for acute respiratory distress syndrome but tested negative for SARS-CoV-2 infection (by PCR and serology).

Laboratory Assays

Plasma levels of complement proteins C3, C3a, C5, C1q, factor H, and sC5b-9 were measured using in house quantitative enzyme-linked immunosorbent assays (ELISAs). Plasma levels of complement factor D and C5a (R&D Systems), human mannose-binding lectin (Hycult Biotech), and C1 inhibitor, Bb, and CIC-C1q (Quidel) were measured using commercially available ELISA kits, following the manufacturer’s instructions. Binding IgG and IgA to SARS-COV-2, and IgG to CCCs (NL63, 229E, OC43, and HKU1) and severe acute respiratory syndrome coronavirus (SARS-CoV) were measured using indirect ELISAs. Plaque reduction neutralization test for SARS-CoV-2 was performed in the University of Pittsburgh Regional Biocontainment Laboratory biosafety level 3 facility. Laboratory procedures are described in detail in the Supplementary Materials.

Statistical Analysis

Comparisons among HCs and ICU and non-ICU patients were performed using nonparametric tests (Mann-Whitney, Kruskal-Wallis, or Fisher exact tests, as appropriate). Spearman correlation was used to measure the degree of association between 2 variables. We performed a receiver operating characteristic curve analysis for complement activation based on the C3a/C3 ratio as a predictor of COVID-19 severity (non-ICU versus ICU). The length of hospital stay was estimated using unadjusted Kaplan-Meier curve analyses. Statistical analyses were not corrected for confounders owing to the small sample size. Data were analyzed with Prism software version 9.0a (GraphPad). Differences were considered significant at P < .05.

We used linear discriminant analysis (LDA) to design a classifier for COVID-19 severity using complement markers and antibody levels. LDA classifiers were obtained using least-squares fitting, shrinkage, and equal priors. The classification errors were estimated using the error count on the training data and bolstered resubstitution estimates using gaussian kernels [24–27]. Log was applied to all features, and missing values were omitted from the analysis.

Study Approvals

All participants signed written informed consent forms for the study before participation. The Institutional Review Board of the University of Pittsburgh reviewed and approved all study protocols for the COVID-19 observational studies (STUDY19050099 and STUDY20040036) and for the prepandemic samples (STUDY19060213).

RESULTS

Markedly Increased Complement Activation in Patients With Severe COVID-19

To investigate the degree of systemic complement activation, we quantified plasma levels of soluble complement proteins and their cleavage products in 55 patients with COVID-19. Baseline demographic and clinical characteristics of 37 critically ill (ICU) and 18 moderately ill (non-ICU) patients are shown in Table 1. As expected, ICU patients presented with higher World Health Organization 10-point ordinal scale of severity scores at admission and during the period of hospitalization. They also had a significantly higher 90-day mortality rate than non-ICU patients (28.2% vs 0%, respectively; P = .005). There was no significant difference in age (P = .16), sex (P = .15), comorbid conditions, and time from symptom onset to sample collection between non-ICU and ICU patients (P = .24).

Table 1.

Descriptive Characteristics of Patients with Coronavirus Disease 2019, Stratified by Disease Severity

CharacteristicsAll Patients (n = 55)Moderately Ill (Non-ICU) Patients (n = 18)Critically Ill (ICU) Patients (n = 37)P Valuea
Age, median (IQR), y62 (53–69)59 (52–66)63 (57–73).16
Male sex, no. (%)30 (54.5)7 (38.9)23 (62.1).15
White race, no. (%)38 (69.1)15 (83.3)23 (67.6).33
Comorbid conditions, no. (%)
 Diabetes mellitus20 (36.4)5 (27.8)15 (40.5).39
 Active neoplasm3 (5.4)2 (11.1)1 (2.7).25
 Obesity (BMI ≥30)b39 (70.9)13 (72.2)26 (70.3)>.99
 Congestive heart failure4 (7.3)0 (0.0)4 (10.8).29
 Chronic obstructive 
pulmonary disease11 (20.0)1 (5.6)10 (27.0).08
 Chronic renal failure8 (14.5)1 (5.6)7 (18.9).25
Smoking (current or former), no (%)27 (49.0)6 (33.3)21 (56.7).15
WHO score, median (IQR)
 Score at admission5 (5–5)4 (4–5)5 (5–5.5)<.001c
 Worst score during 
hospitalization6 (5–9)5 (5–5)9 (6–10)<.001c
90-d Mortality, no. (%)12 (21.8)0 (0.0)12 (32.4).005c
Time from symptom onset to sampling, median (IQR), d10 (5–12)7 (5–11)10 (6–13).24
CharacteristicsAll Patients (n = 55)Moderately Ill (Non-ICU) Patients (n = 18)Critically Ill (ICU) Patients (n = 37)P Valuea
Age, median (IQR), y62 (53–69)59 (52–66)63 (57–73).16
Male sex, no. (%)30 (54.5)7 (38.9)23 (62.1).15
White race, no. (%)38 (69.1)15 (83.3)23 (67.6).33
Comorbid conditions, no. (%)
 Diabetes mellitus20 (36.4)5 (27.8)15 (40.5).39
 Active neoplasm3 (5.4)2 (11.1)1 (2.7).25
 Obesity (BMI ≥30)b39 (70.9)13 (72.2)26 (70.3)>.99
 Congestive heart failure4 (7.3)0 (0.0)4 (10.8).29
 Chronic obstructive 
pulmonary disease11 (20.0)1 (5.6)10 (27.0).08
 Chronic renal failure8 (14.5)1 (5.6)7 (18.9).25
Smoking (current or former), no (%)27 (49.0)6 (33.3)21 (56.7).15
WHO score, median (IQR)
 Score at admission5 (5–5)4 (4–5)5 (5–5.5)<.001c
 Worst score during 
hospitalization6 (5–9)5 (5–5)9 (6–10)<.001c
90-d Mortality, no. (%)12 (21.8)0 (0.0)12 (32.4).005c
Time from symptom onset to sampling, median (IQR), d10 (5–12)7 (5–11)10 (6–13).24

Abbreviations: BMI, body mass index; ICU, intensive care unit; IQR, interquartile range; WHO, World Health Organization.

Nonparametric test comparisons were performed using Mann-Whitney and Fisher exact tests for continuous and categorical variables, respectively.

BMI calculated as weight in kilograms divided by height in meters squared.

Significant at P < .05.

Table 1.

Descriptive Characteristics of Patients with Coronavirus Disease 2019, Stratified by Disease Severity

CharacteristicsAll Patients (n = 55)Moderately Ill (Non-ICU) Patients (n = 18)Critically Ill (ICU) Patients (n = 37)P Valuea
Age, median (IQR), y62 (53–69)59 (52–66)63 (57–73).16
Male sex, no. (%)30 (54.5)7 (38.9)23 (62.1).15
White race, no. (%)38 (69.1)15 (83.3)23 (67.6).33
Comorbid conditions, no. (%)
 Diabetes mellitus20 (36.4)5 (27.8)15 (40.5).39
 Active neoplasm3 (5.4)2 (11.1)1 (2.7).25
 Obesity (BMI ≥30)b39 (70.9)13 (72.2)26 (70.3)>.99
 Congestive heart failure4 (7.3)0 (0.0)4 (10.8).29
 Chronic obstructive 
pulmonary disease11 (20.0)1 (5.6)10 (27.0).08
 Chronic renal failure8 (14.5)1 (5.6)7 (18.9).25
Smoking (current or former), no (%)27 (49.0)6 (33.3)21 (56.7).15
WHO score, median (IQR)
 Score at admission5 (5–5)4 (4–5)5 (5–5.5)<.001c
 Worst score during 
hospitalization6 (5–9)5 (5–5)9 (6–10)<.001c
90-d Mortality, no. (%)12 (21.8)0 (0.0)12 (32.4).005c
Time from symptom onset to sampling, median (IQR), d10 (5–12)7 (5–11)10 (6–13).24
CharacteristicsAll Patients (n = 55)Moderately Ill (Non-ICU) Patients (n = 18)Critically Ill (ICU) Patients (n = 37)P Valuea
Age, median (IQR), y62 (53–69)59 (52–66)63 (57–73).16
Male sex, no. (%)30 (54.5)7 (38.9)23 (62.1).15
White race, no. (%)38 (69.1)15 (83.3)23 (67.6).33
Comorbid conditions, no. (%)
 Diabetes mellitus20 (36.4)5 (27.8)15 (40.5).39
 Active neoplasm3 (5.4)2 (11.1)1 (2.7).25
 Obesity (BMI ≥30)b39 (70.9)13 (72.2)26 (70.3)>.99
 Congestive heart failure4 (7.3)0 (0.0)4 (10.8).29
 Chronic obstructive 
pulmonary disease11 (20.0)1 (5.6)10 (27.0).08
 Chronic renal failure8 (14.5)1 (5.6)7 (18.9).25
Smoking (current or former), no (%)27 (49.0)6 (33.3)21 (56.7).15
WHO score, median (IQR)
 Score at admission5 (5–5)4 (4–5)5 (5–5.5)<.001c
 Worst score during 
hospitalization6 (5–9)5 (5–5)9 (6–10)<.001c
90-d Mortality, no. (%)12 (21.8)0 (0.0)12 (32.4).005c
Time from symptom onset to sampling, median (IQR), d10 (5–12)7 (5–11)10 (6–13).24

Abbreviations: BMI, body mass index; ICU, intensive care unit; IQR, interquartile range; WHO, World Health Organization.

Nonparametric test comparisons were performed using Mann-Whitney and Fisher exact tests for continuous and categorical variables, respectively.

BMI calculated as weight in kilograms divided by height in meters squared.

Significant at P < .05.

Plasma levels of C3 were markedly reduced in ICU patients compared with non-ICU patients and HCs, indicating increased C3 consumption in patients with severe COVID-19. C3a levels were increased in patients with COVID-19 compared with HCs and approximately 100-fold greater in patients with severe COVID-19 compared with non-ICU patients. Similarly, the ratio between C3a and C3 was markedly increased in ICU patients (Figure 1A). C5 levels were significantly reduced in ICU patients, and the levels of C5a and C5a/C5 ratio were significantly increased in ICU compared with non-ICU patients (Figures 1B). Plasma levels of sC5b-9, the activation product of the terminal complement cascade, were significantly elevated in patients with COVID-19 relative to those in HCs (Figure 1C), though no significant differences were evident based on illness severity.

Markers of complement activation are increased in patients with coronavirus disease 2019 (COVID-19) and are associated with disease severity. Plasma levels of various markers of complement activation were quantified in critically ill patients hospitalized in intensive care units (ICU patients; n = 37) and moderately ill patients hospitalized in dedicated COVID-19 wards (non-ICU patients; n = 18). Plasma samples from healthy controls (HCs; n = 10) were included as a reference population for comparisons. Scatterplots show the differences between HCs and non-ICU and ICU patients in the following plasma levels and ratios: C3, C3a, and C3a/C3 (A); C5, C5a, and C5a/C5 ratio (B); and soluble C5b-9 (sC5b-9) (C). Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Statistical significance was determined using Kruskal-Wallis 1-way analysis of variance (#P < .05; ##P < .01; ###P < .001; ####P < .0001) or the Mann-Whitney U test (*P < .05; **P < .01; ***P < .001; ****P < .0001).
Figure 1.

Markers of complement activation are increased in patients with coronavirus disease 2019 (COVID-19) and are associated with disease severity. Plasma levels of various markers of complement activation were quantified in critically ill patients hospitalized in intensive care units (ICU patients; n = 37) and moderately ill patients hospitalized in dedicated COVID-19 wards (non-ICU patients; n = 18). Plasma samples from healthy controls (HCs; n = 10) were included as a reference population for comparisons. Scatterplots show the differences between HCs and non-ICU and ICU patients in the following plasma levels and ratios: C3, C3a, and C3a/C3 (A); C5, C5a, and C5a/C5 ratio (B); and soluble C5b-9 (sC5b-9) (C). Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Statistical significance was determined using Kruskal-Wallis 1-way analysis of variance (#P < .05; ##P < .01; ###P < .001; ####P < .0001) or the Mann-Whitney U test (*P < .05; **P < .01; ***P < .001; ****P < .0001).

We next performed a receiver operating characteristic analysis for C3a/C3 ratio as a predictor for COVID-19 severity (Figure 2A). A threshold for C3a/C3 ratio of ≥4.5 (log10) classified patients in high and low complement activation and provided a sensitivity and specificity for COVID-19 severity of 94.4% and 70.6%, respectively. We then analyzed the relationship between complement activation, viral RNA levels, and clinical laboratory markers of inflammation. Data on viral load was available for a total of 45 patients (16 non-ICU and 29 ICU patients) [28]. We found that plasma levels of SARS-CoV-2 RNA were highest in patients classified as having high complement activation (Figure 2B). Patients with high C3a/C3 ratios also had significantly higher plasma levels of proinflammatory cytokines IL-6 and IL-8 (Figure 2C) and procalcitonin and suppression of tumorigenicity (Supplementary Figure 1). Levels of fractalkine, pentraxin 3, tumor necrosis factor receptor 1, and receptor for advanced glycation end products were not associated with C3a/C3 ratios (Supplementary Figure 1).

Relationship between complement activation, viral RNA levels, and clinical laboratory markers of inflammation. A, Receiver operating characteristic curve (ROC) analysis for C3a/C3 ratio as a predictor for coronavirus disease 2019 (COVID-19) severity. A C3a/C3 ratio threshold ≥4.5 log10 provided a sensitivity of 94.4% and a specificity of 7.06%. Abbreviations: AUC area under the ROC curve; CI, confidence interval. B, C, Scatterplots of the differences in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plasma RNA (B) and inflammatory cytokine (C) levels between patients classified as having low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Abbreviations: IL-6, IL-8, and IL-10, interleukin 6, 8, and 10. *P < .05 (Mann-Whitney U test). D, Kaplan-Meier curves for 90-day survival and time to discharge from hospital admission, stratified by low versus high complement activation.
Figure 2.

Relationship between complement activation, viral RNA levels, and clinical laboratory markers of inflammation. A, Receiver operating characteristic curve (ROC) analysis for C3a/C3 ratio as a predictor for coronavirus disease 2019 (COVID-19) severity. A C3a/C3 ratio threshold ≥4.5 log10 provided a sensitivity of 94.4% and a specificity of 7.06%. Abbreviations: AUC area under the ROC curve; CI, confidence interval. B, C, Scatterplots of the differences in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plasma RNA (B) and inflammatory cytokine (C) levels between patients classified as having low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Abbreviations: IL-6, IL-8, and IL-10, interleukin 6, 8, and 10. *P < .05 (Mann-Whitney U test). D, Kaplan-Meier curves for 90-day survival and time to discharge from hospital admission, stratified by low versus high complement activation.

We then assessed the relationship between complement activation and clinical parameters of disease severity. High levels of complement activation were significantly associated with worse 90-day survival (log-rank P = .004) and longer time to discharge from the hospital (log-rank P = .002; Figure 2D). When restricted to the group of patients with severe COVID-19, high complement activation was associated with an increased risk of intubation (odds ratio, 7.00 [95% confidence interval, 1.47–29.21; P= .02) and with a greater 90-day mortality risk (7.03 [1.38–34.4]; P= .02). Collectively, these data indicate that the complement cascade is activated in patients hospitalized with COVID-19 and that overactivation of the complement system is correlated with greater viral burden, increased systemic inflammation, and worse disease outcomes.

High Levels of Classical Pathway Activation Correlated With COVID-19 Severity

We next sought to determine which arms of the system contribute to the observed overactivation of complement in patients with severe COVID-19. We first quantified markers that are specific to the initial activation steps of the alternative pathway. We observed no differences in the plasma levels of factors D and Bb, a byproduct of factor B cleavage, between patients with COVID-19 and HCs (Figure 3A). Subsequently, we found that plasma levels of factor H, a regulatory protein of the alternative pathway, were increased approximately 2-fold in ICU compared with non-ICU patients. Factor H–factor D ratios, which act in antagonistic ways to influence the activation/deactivation of the alternative pathway, were higher in both non-ICU and ICU patients compared with HCs, suggesting that the alternative pathway is being suppressed in patients with COVID-19. We also found that plasma levels of functional mannose-binding lectin, the main recognition molecule of the lectin pathway, were similar among HCs and ICU and non-ICU patients (Figure 3B), indicating that this arm of the complement system is not being disproportionally activated in patients with COVID-19.

Markers of complement activation by the classical pathway are significantly increased in patients with severe coronavirus disease 2019 (COVID-19). Plasma levels of markers of complement initiation and regulatory proteins unique to each individual arm (alternative, lectin, and classical pathways) of the complement system were quantified in critically ill patients hospitalized in intensive care units (ICU patients; n = 37) and moderately ill patients hospitalized in dedicated COVID-19 wards (non-ICU patients; n = 18). Plasma samples from healthy controls (HCs; n = 10) were included as a reference population for comparisons. A–C, Scatterplots of the differences in plasma levels of markers of the alternative pathway, including factor D, Bb, factor H, and factor H–factor D ratio (A); functional mannose-binding lectin (MBL) of the lectin pathway (B); and markers of the classical pathway, (C) including total human C1q, circulating immune complexes (CICs) binding to C1q (CIC-C1q), and C1 inhibitor (C). D, Scatterplots of the differences in CIC-C1q between patients classified as having low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Statistical significance was determined using Kruskal-Wallis 1-way analysis of variance (###P < .001; ####P < .0001) or the Mann-Whitney U test (*P < .05; **P < .01; ***P < .001).
Figure 3.

Markers of complement activation by the classical pathway are significantly increased in patients with severe coronavirus disease 2019 (COVID-19). Plasma levels of markers of complement initiation and regulatory proteins unique to each individual arm (alternative, lectin, and classical pathways) of the complement system were quantified in critically ill patients hospitalized in intensive care units (ICU patients; n = 37) and moderately ill patients hospitalized in dedicated COVID-19 wards (non-ICU patients; n = 18). Plasma samples from healthy controls (HCs; n = 10) were included as a reference population for comparisons. A–C, Scatterplots of the differences in plasma levels of markers of the alternative pathway, including factor D, Bb, factor H, and factor H–factor D ratio (A); functional mannose-binding lectin (MBL) of the lectin pathway (B); and markers of the classical pathway, (C) including total human C1q, circulating immune complexes (CICs) binding to C1q (CIC-C1q), and C1 inhibitor (C). D, Scatterplots of the differences in CIC-C1q between patients classified as having low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Statistical significance was determined using Kruskal-Wallis 1-way analysis of variance (###P < .001; ####P < .0001) or the Mann-Whitney U test (*P < .05; **P < .01; ***P < .001).

We next assessed whether complement system initiation through the classical pathway accounted for the increased complement activation observed in patients with COVID-19. We observed an increase in plasma levels of total C1q in patients with COVID-19, which were significantly higher in ICU than in non-ICU patients (Figure 3C). In addition, we found that levels of circulating immune complexes (CICs) with the ability to bind to functional C1q (CIC-C1q) were markedly increased in patients with COVID-19 and were greater in ICU than in non-ICU patients. Furthermore, levels of CIC-C1q were 2-fold greater in patients with high complement activation (Figure 3D). We found no differences in levels of functional C1 inhibitor, a plasma protease that regulates activation of the classical pathway, between patients with COVID-19 and HCs (Figure 3C). Collectively, these findings suggest that the classical pathway is the main factor leading to the overactivation of the complement system in patients with severe COVID-19.

Increased Antibody Titers Against SARS-CoV-2 in Patients With Severe COVID-19 and Association With Greater Complement Activation

We subsequently analyzed the interactions between the antibody response to SARS-CoV-2 and complement activation. We first developed ELISAs to estimate titers of binding IgG and IgA against SARS-CoV-2 structural antigens (Supplementary Table 2 and Supplementary Figures 2 and 3). Similar to previous reports [29–31], we observed that patients with severe COVID-19 mounted a faster and stronger binding IgG and IgA response against the spike protein and NP (Figure 4A) and neutralizing antibody response (Figure 4D) against SARS-CoV-2, compared with patients who had moderate disease. Notably, IgG, IgA, and neutralizing antibody titers (Figure 4B–4D) trended higher in patients with severe COVID-19 even relatively soon after the onset of illness. These findings point to an early and robust antibody response as one determinant of COVID-19 severity.

Binding immunoglobulin (Ig) G and IgA, and neutralizing antibody titers against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are markedly increased in patients with severe coronavirus disease 2019. A–C, Scatterplots of the differences in plasma titers of binding IgG and IgA against the spike protein and the nucleoprotein (NP) of SARS-CoV-2 between intensive care unit (ICU) and non-ICU patients, regardless of time from illness onset (A) and stratified by days after symptom onset (B, C). D, Titers of neutralizing antibodies estimated by plaque reduction neutralization test (PRNT), regardless of time from illness onset and stratified by days after symptom onset. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. *P < .05; **P < .01; ***P < .001 (Mann-Whitney U test).
Figure 4.

Binding immunoglobulin (Ig) G and IgA, and neutralizing antibody titers against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are markedly increased in patients with severe coronavirus disease 2019. A–C, Scatterplots of the differences in plasma titers of binding IgG and IgA against the spike protein and the nucleoprotein (NP) of SARS-CoV-2 between intensive care unit (ICU) and non-ICU patients, regardless of time from illness onset (A) and stratified by days after symptom onset (B, C). D, Titers of neutralizing antibodies estimated by plaque reduction neutralization test (PRNT), regardless of time from illness onset and stratified by days after symptom onset. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. *P < .05; **P < .01; ***P < .001 (Mann-Whitney U test).

We next investigated the degree of association between the antibody response and complement system activation. Levels of CIC-C1q were positively correlated with binding titers of IgG against the spike protein and NP and with neutralizing antibody titers to SARS-CoV-2 (Figure 5A), particularly in the ICU patients. Levels of CIC-C1q were not correlated with binding IgA titers, as this immunoglobulin class does not activate complement via the classical pathway. In line with this, the IgG/IgA binding titer ratio was greater in ICU than in non-ICU patients (Figure 5B). Substantiating our hypotheses of complement overactivation by the classical pathway, we also found that IgG titers were markedly increased in patients with COVID-19 classified as having high levels of complement activation, while levels of IgA did not differ between groups (Figure 5C). Finally, a 2-dimensional LDA classifier using C3a/C3 ratio and CIC-C1q levels distinguished non-ICU and ICU patients with an error <0.2 and with an estimated sensitivity and specificity of 73.5% and 94.4%, respectively (Figure 5D). Collectively, these data provided strong evidence that COVID-19 severity is associated with a predominant IgG response and with high levels of CIC-C1q that ultimately lead to overactivation of the complement system.

Binding immunoglobulin (Ig) G and IgA, and neutralizing antibody titers against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are associated with complement activation in patients with severe coronavirus disease 2019 (COVID-19). A, Degree of correlation between plasma levels of circulating immune complexes (CIC-C1q) and binding IgG against the spike and nucleoprotein (NP) and neutralizing antibody titers, assessed using the Spearman rank correlation coefficient test. B, Scatterplots of the differences in IgG/IgA ratios between intensive care unit (ICU) and non-ICU patients with COVID-19. Abbreviation: PRNT, plaque reduction neutralization test. C, Scatterplots of binding titers of IgG against spike protein and NP, neutralizing antibody titers, and binding titers of IgA against spike protein and NP, stratified by the C3a/C3 ratio as low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. *P < .05; **P < .01 (Mann-Whitney U test). D, Linear discriminant analysis using all measurements for complement and SARS-CoV-2 antibody, showing a 2-dimensional discriminator using CIC-C1q and the C3a/C3 ratio.
Figure 5.

Binding immunoglobulin (Ig) G and IgA, and neutralizing antibody titers against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are associated with complement activation in patients with severe coronavirus disease 2019 (COVID-19). A, Degree of correlation between plasma levels of circulating immune complexes (CIC-C1q) and binding IgG against the spike and nucleoprotein (NP) and neutralizing antibody titers, assessed using the Spearman rank correlation coefficient test. B, Scatterplots of the differences in IgG/IgA ratios between intensive care unit (ICU) and non-ICU patients with COVID-19. Abbreviation: PRNT, plaque reduction neutralization test. C, Scatterplots of binding titers of IgG against spike protein and NP, neutralizing antibody titers, and binding titers of IgA against spike protein and NP, stratified by the C3a/C3 ratio as low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. *P < .05; **P < .01 (Mann-Whitney U test). D, Linear discriminant analysis using all measurements for complement and SARS-CoV-2 antibody, showing a 2-dimensional discriminator using CIC-C1q and the C3a/C3 ratio.

High IgG Titers Against the NP of CCCs and Complement Overactivation in Patients With COVID-19

To obtain insights into the apparent rapid increase in antibody response to SARS-CoV-2 in patients with severe COVID-19, we analyzed the IgG responses to the antigenically related CCCs. We developed binding assays to estimate IgG titers toward a range of spike and NP antigens of CCCs and the spike protein of SARS-CoV. We validate our binding assays using a panel of remnants samples from a pediatric seroprevalence study (in patients aged 1–16 years) living in the same region as the adult patients with COVID-19 [32] (Supplementary Tables 3 and 4 and Supplementary Figure 4). The frequency of binding IgG against all CCCs increased with age (Supplementary Table 5), which is consistent with endemic circulation of these viruses in the United States.

We then quantified IgG titers against CCC antigens in our patients with COVID-19. All patients with COVID-19 exhibited substantial IgG reactivity to ≥1 of the CCC antigens (Supplementary Table 6). Notably, the increase in IgG titer to CCC antigens was more prominent in ICU than in non-ICU patients and most striking against NP antigens (Figure 6A). We also noted an increase in IgG titers against SARS-CoV spike protein in patients with COVID-19, especially in the ICU patients. This likely reflects cross-reactive antibody responses, given the structural similarities between SARS-CoV and SARS-CoV-2 [33].

Binding immunoglobulin (Ig) G titers against various antigens of the common cold coronaviruses (CCCs) are higher in patients with severe coronavirus disease 2019 (COVID-19) and are also associated with complement activation. A, Scatterplots of the differences in plasma titers of binding IgG against the spike protein and the nucleoprotein (NP) of various CCCs (NL63, 229E, OC43, and HKU1) and the spike protein of severe acute respiratory syndrome coronavirus (SARS-CoV) between intensive care unit (ICU) and non-ICU patients. Plasma samples from patients who were hospitalized in ICU wards with respiratory illness not related to severe acute respiratory syndrome coronavirus (SARS-CoV-2) infection were included as a group for comparisons (non–COVID 19 ICU; n = 20). Scatterplots of the differences in plasma titers of binding IgG against the spike protein of CCCs and of SARS-CoV (B) and against the NP of the CCCs (C), stratified by the C3a/C3 ratio as low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Statistical significance was determined using Kruskal-Wallis 1-way analysis of variance (#P < .05; ##P < .01; ###P < .001) and the Mann-Whitney U test (*P < .05; **P < .01; ***P < .001). D, Degree of correlation between plasma levels of circulating immune complexes (CIC-C1q) and IgG titers to CCC NP, assessed using the Spearman rank correlation coefficient test.
Figure 6.

Binding immunoglobulin (Ig) G titers against various antigens of the common cold coronaviruses (CCCs) are higher in patients with severe coronavirus disease 2019 (COVID-19) and are also associated with complement activation. A, Scatterplots of the differences in plasma titers of binding IgG against the spike protein and the nucleoprotein (NP) of various CCCs (NL63, 229E, OC43, and HKU1) and the spike protein of severe acute respiratory syndrome coronavirus (SARS-CoV) between intensive care unit (ICU) and non-ICU patients. Plasma samples from patients who were hospitalized in ICU wards with respiratory illness not related to severe acute respiratory syndrome coronavirus (SARS-CoV-2) infection were included as a group for comparisons (non–COVID 19 ICU; n = 20). Scatterplots of the differences in plasma titers of binding IgG against the spike protein of CCCs and of SARS-CoV (B) and against the NP of the CCCs (C), stratified by the C3a/C3 ratio as low (C3a/C3 ratio, <4.5 log10) or high (C3a/C3 ratio, ≥4.5 log10) complement activation. Each symbol in the scatterplots represents an individual patient, and horizontal lines represent medians. Statistical significance was determined using Kruskal-Wallis 1-way analysis of variance (#P < .05; ##P < .01; ###P < .001) and the Mann-Whitney U test (*P < .05; **P < .01; ***P < .001). D, Degree of correlation between plasma levels of circulating immune complexes (CIC-C1q) and IgG titers to CCC NP, assessed using the Spearman rank correlation coefficient test.

We next evaluated the relationship between antibody binding to CCCs and complement activation and found that IgG titers to NP antigens (Figure 6C) but not to the spike (Figure 6B) were greater among patients classified as having high complement activation. In addition, patients with COVID-19 displayed a significant positive correlation between levels of CIC-C1q and IgG titers against NP antigens for all CCCs (Figure 6D). This correlation was not significant when stratifying patients by disease severity. Collectively, these findings reveal a cross-reactive anamnestic response to CCCs in patients with COVID-19 that is correlated with increased complement activation.

DISCUSSION

We show that IgG responses against SARS-CoV-2 and CCCs are associated with complement activation in patients with severe COVID-19. In line with earlier reports [1, 6–9], we confirm that markers of systemic complement activation are significantly increased in patients with severe COVID-19 and distinguish those with worse disease outcomes. More importantly, our in-depth analysis of pathway-specific components delineate the role of the classical pathway in the enhanced complement activation observed in patients with severe COVID-19 and suggests a previously unrecognized mechanism of antibody-mediated activation driven by an overexuberant IgG response against SARS-CoV-2 and CCC antigens. These data indicate that an immune dysregulation of both the innate and adaptive immune systems is a primary feature of COVID-19 severity.

Our data showed that high complement activation is associated with some established markers of inflammation commonly elevated in patients with severe COVID-19 [34, 35]. This is consistent with a hyperinflammatory milieu characteristic of SARS-CoV-2 pathophysiology, which is likely influenced by the enhanced levels of the potent anaphylatoxins C3a and C5a, particularly in ICU patients [6, 8]. In contrast with a prior report [1], our findings indicate that complement overactivation is correlated with levels of SARS-CoV-2 RNA in the plasma. This discrepancy between studies might be related to the sensitivity of the assays used for detection of viral RNA. SARS-CoV-2 RNA in our patients was measured using an ultrasensitive quantitative reverse-transcription PCR method, resulting in a greater proportion of RNAemia-positive ICU patients than reported elsewhere [28].

While earlier reports indicate that overactivation of the complement system plays a critical role in COVID-19 immunopathology [1, 5–10], our detailed analysis of several pathway-specific markers extend these findings by revealing that the classical pathway and the humoral response are associated with unfavorable COVID-19 outcomes. The high levels of CIC-C1q seen in our patients are likely due to the presence of greater amounts of antigens and antibodies present in the circulation, reduced CIC clearance, or both. Earlier reports have shown extensive deposition of C1q and C4d, a marker of complement activation by CICs, in postmortem lung sections of patients with COVID-19 [10, 36]. Complement activation by CICs has also been associated with the pathophysiology of other respiratory viruses [37, 38]. Interestingly, a prior in vitro study reported no C1q binding when mixing SARS-CoV-2 proteins with nonimmune human serum, indicating that the classical pathway is not activated in the absence of specific antibodies [39]. Of note, the assay used in our study measures CICs with the ability to bind to functional C1q and thus activate the complement system—a feature that is further confirmed by the association between levels of CIC-C1q and C3a/C3 ratio seen in our patients. In addition, it is important to mention that CICs can also induce a complement-independent dysregulated inflammatory milieu via CIC engagement to FcγR, leading to excessive activation of alveolar macrophages that likely contributes to amplification of the proinflammatory cytokine release seen in patients with severe COVID-19 [40, 41].

Consistent with earlier reports [29–31], our analysis of the magnitude of the antibody response to SARS-CoV-2 suggest that ICU patients seroconvert earlier and developed higher titers of binding IgG and IgA and neutralizing antibodies than patients with moderate disease. Interestingly, patients with higher IgG/IgA ratios had more severe disease, which is consistent with the greater capacity of IgG to mediate complement activation. Patients with severe COVID-19 also displayed a strong rise in the levels of IgG against the NP of CCCs, in agreement with 2021 reports showing that IgG antibodies and B cell clones targeting structural antigens of CCCs are boosted on SARS-CoV-2 infection [20, 42, 43]. Extending on these findings, our data reveal a marked correlation between enhanced complement activation by CIC-C1q and IgG levels to SARS-CoV-2 and CCCs. These findings suggest that overactivation of the complement system by CICs at the pulmonary and systemic levels may contribute to the secretion of proinflammatory cytokines, immune cell recruitment, and complement deposition within the lung, in a mechanism similar to the antibody-mediated enhanced respiratory disease observed in other respiratory viruses [38, 44].

Because activation markers of the alternative and the lectin pathways were present at similar levels, we were unable to ascribe their role in the observed complement overactivation seen in our patients with COVID-19. However, this does not exclude the possibility that these arms of the complement system also contribute to COVID-19 severity. Earlier reports have shown the ability of mannose-binding lectin to bind SARS-CoV-2 proteins and activate complement [39, 45], suggesting that the lectin pathway is possibly involved in complement activation in COVID-19. Activation of the alternative pathway has also been implicated in COVID-19 pathogenesis using in vitro systems and in clinical studies [9, 46, 47]. Interestingly, recent reports have suggested that activation of the alternative pathway in SARS-CoV-2–infected patients is factor D independent [9, 46]. This is to some extent consistent with our data showing no differences in factor D levels between patients with COVID-19 and HCs. Notably, negative regulatory mediators of the alternative pathway (eg, factor H) were markedly up-regulated in our patients, suggesting that this pathway is being suppressed rather than synergistically amplifying complement activation in COVID-19. This variability among studies may result from differences in timing of sample collection, type of sample tested (plasma EDTA, citrate, and serum), and sensitivity of the methods used for complement measurements. However, based on available data, it is reasonable to hypothesize that SARS-CoV-2 infection triggers activation of the complement cascade via all 3 pathways and that the classical pathway likely contributes the most to overactivation of the system and hence disease severity.

We acknowledge that the number of patients included in our analysis was relatively small, which might underpower some of the conclusions. In addition, all patients enrolled in our cohort were hospitalized, and complement activation and antibody response in asymptomatic or mild COVID-19 cases could not be assessed in this study. More extensive studies are needed to clearly elucidate the mechanistic basis by which SARS-CoV-2 mediates excessive activation of the complement system, ideally with longitudinal sampling throughout the course of the disease, including samples from patients with mild COVID-19 and from those who develop moderate symptoms and later progress to severe disease.

In summary, our findings highlight the role of the classical pathway in mediating complement overactivation in patients with severe COVID-19 and reveal that the interplay between IgG responses to SARS-CoV-2 and to CCCs in SARS-CoV-2–infected patients may be associated with adverse outcomes. Although the viral and immunoglobulin content of the CICs needs to be determined, our findings provide new insights into the immunopathogenesis of SARS-CoV-2 infection and indicate that IgG-mediated effector functions may contribute to COVID-19 severity. Importantly, such identification of immune parameters might help in accurately selecting interventional strategies, including inhibition of the C4bC2a and C3bBb convertases, blockers of C3a and C5a receptors, and a variety of others reviewed elsewhere [48].

Supplementary Data

Supplementary materials are available at The Journal of Infectious Diseases online. Supplementary materials consist of data provided by the author that are published to benefit the reader. The posted materials are not copyedited. The contents of all supplementary data are the sole responsibility of the authors. Questions or messages regarding errors should be addressed to the author.

Notes

Acknowledgments. We thank all the clinical staff at the University of Pittsburgh Medical Center Presbyterian and Shadyside Hospital units for coordinating patient enrollment and specimen collection. We also thank Willian Bain, MD for assistance with collection of patient samples.

Author contributions. Conceptualization: P. M. S. C., S. M. B. B., D. B., and E. T. A. M. Investigation: P. M. S. C., D. J. T., G. D. K., J. L. J., U. B. N., M. D., S. R., M. M. M., S. W., A. N., B. S., J. M., C. A. F., S. M. B. B., D. B., and E. T. A. M. Funding acquisition: P. M. S. C., G. D. K., S. W., J. M., B. J. M., B. J. C. M., D. B., and E. T. A. M. Clinical sample acquisition: G. D. K., A. M., B. J. M., F. S., C. S., B. J. C. M., and B. M. Writing—original draft: P. M. S. C., S. M. B. B., D. B., and E. T. A. M. Writing—review and editing: P. M. S. C., G. D. K., J. L. J., S. W., A. M., B. J. M., F. S., S. M. B. B., D. B., and E. T. A. M.

Financial support. This study was supported in part by the University of Pittsburgh Clinical and Translational Science Institute (pilot coronavirus disease 2019 awards to P. M. S. C., G. D. K., S. W., J. M., D. B., and E. T. A. M.); the DSF Charitable Foundation (G. D. K. and S. W.); the National Heart, Lung, and Blood Institute, National Institutes of Health (K23 grant 139987 to G. D. K. and grant to B. J. M.); the University of Pittsburgh Medical Center (internal grant to B. J. C. M.); and the National Institute of Allergy and Infectious Diseases (Pittsburgh AIDS Research Training Fellowship–National Research Service Award to 5T32AI065380-15 to D. J. T.).

References

1.

Holter
JC
,
Pischke
SE
,
Boer de
E
, et al.
Systemic complement activation is associated with respiratory failure in COVID-19 hospitalized patients.
Proc Natl Acad Sci U S A
2020
;
117
:
25018
25
.

2.

Fang
S
,
Wang
H
,
Lu
L
,
Jia
Y
,
Xia
Z.
Decreased complement C3 levels are associated with poor prognosis in patients with COVID-19: a retrospective cohort study.
Int Immunopharmacol
2020
;
89
:
107070
.

3.

Ramlall
V
,
Thangaraj
PM
,
Meydan
C
, et al.
Immune complement and coagulation dysfunction in adverse outcomes of SARS-CoV-2 infection.
Nat Med
2020
;
26
:
1609
15
.

4.

Zhang
J
,
Wang
Z
,
Wang
X
,
Hu
Z
,
Yang
C
,
Lei
P.
Risk factors for mortality of COVID-19 patient based on clinical course: a single center retrospective case-control study.
Front Immunol
2021
;
12
:1–11.

5.

Henry
BM
,
Szergyuk
I
,
Oliveira
MHS
, et al.
Complement levels at admission as a reflection of coronavirus disease 2019 (COVID-19) severity state.
J Med Virol
2021
;
93
:
5515
22
.

6.

Nooijer
AH
de,
Grondman
I
,
Janssen
NAF
, et al.
Complement activation in the disease course of coronavirus disease 2019 and its effects on clinical outcomes.
J Infect Dis
2021
;
223
:
214
24
.

7.

Sinkovits
G
,
Mező
B
,
Réti
M
, et al.
Complement overactivation and consumption predicts in-hospital mortality in SARS-CoV-2 infection.
Front Immunol
2021
;
12
:1–13.

8.

Alosaimi
B
,
Mubarak
A
,
Hamed
ME
, et al.
Complement anaphylatoxins and inflammatory cytokines as prognostic markers for COVID-19 severity and in-hospital mortality.
Front Immunol
2021
;
12
:1–13.

9.

Ma
L
,
Sahu
SK
,
Cano
M
, et al.
Increased complement activation is a distinctive feature of severe SARS-CoV-2 infection.
Sci Immunol
2021
;
6
:1–12.

10.

Macor
P
,
Durigutto
P
,
Mangogna
A
, et al.
Multiple-organ complement deposition on vascular endothelium in COVID-19 patients.
Biomedicines
2021
;
9
:
1003
.

11.

Stenmark
KR
,
Frid
MG
,
Gerasimovskaya
E
, et al.
Mechanisms of SARS-CoV-2-induced lung vascular disease: potential role of complement.
Pulm Circ
2021
;
11
:
1204589402110157
.

12.

Sethi
S
,
D’Costa
MR
,
Hermann
SM
,
Nasr
SH
,
Fervenza
FC.
Immune-complex glomerulonephritis after COVID-19 infection.
Kidney Int Reports
2021
;
6
:
1170
3
.

13.

Zhao
Y
,
Nie
HX
,
Hu
K
, et al.
Abnormal immunity of non-survivors with COVID-19: predictors for mortality.
Infect Dis Poverty
2020
;
9
:
108
.

14.

Manzo
G.
COVID-19 as an immune complex hypersensitivity in antigen excess conditions: theoretical pathogenetic process and suggestions for potential therapeutic interventions.
Front Immunol
2020
;
11
:1–6.

15.

Gong
F
,
Dai
Y
,
Zheng
T
, et al.
Peripheral CD4+ T cell subsets and antibody response in COVID-19 convalescent individuals.
J Clin Invest
2020
;
130
:
6588
99
.

16.

Mateus
J
,
Grifoni
A
,
Tarke
A
, et al.
Selective and cross-reactive SARS-CoV-2 T cell epitopes in unexposed humans.
Science
2020
;
370
:
89
94
.

17.

Kuri-Cervantes
L
,
Pampena
MB
,
Meng
W
, et al.
Comprehensive mapping of immune perturbations associated with severe COVID-19.
Sci Immunol
2020
;
5
:1–15.

18.

Henss
L
,
Scholz
T
,
Rhein von
C
, et al.
Analysis of humoral immune responses in patients with severe acute respiratory syndrome coronavirus 2 infection.
J Infect Dis
2021
;
223
:
56
61
.

19.

Ng
KW
,
Faulkner
N
,
Cornish
GH
, et al.
Preexisting and de novo humoral immunity to SARS-CoV-2 in humans.
Science
2020
;
370
:
1339
43
.

20.

Aguilar-Bretones
M
,
Westerhuis
BM
,
Raadsen
MP
, et al.
Seasonal coronavirus–specific B cells with limited SARS-CoV-2 cross-reactivity dominate the IgG response in severe COVID-19.
J Clin Invest
2021
;
131
:1–13.

21.

Bain
W
,
Yang
H
,
Shah
FA
, et al.
COVID-19 versus Non–COVID-19 acute respiratory distress syndrome: comparison of demographics, physiologic parameters, inflammatory biomarkers, and clinical outcomes.
Ann Am Thorac Soc
2021
;
18
:
1202
10
.

22.

Drohan
CM
,
Nouraie
SM
,
Bain
W
, et al.
Biomarker-based classification of patients with acute respiratory failure into inflammatory subphenotypes: a single-center exploratory study.
Crit Care Explor
2021
;
3
:
e0518
.

23.

Bain
W
,
Li
H
,
Geest van der
R
, et al.
Increased alternative complement pathway function and improved survival during critical illness.
Am J Respir Crit Care Med
2020
;
202
:
230
40
.

24.

Virtanen
P
,
Gommers
R
,
Oliphant
TE
, et al.
SciPy 1.0: fundamental algorithms for scientific computing in Python.
Nat Methods
2020
;
17
:
261
72
.

25.

Pedregosa
F
,
Varoquaux
G
,
Gramfort
A
, et al.
Scikit-Learn: machine learning in Python.
J Mach Learn Res
2011
;
12
:
2825
30
.

26.

Reback
J
,
jbrockmendel
,
McKinney
W
, et al.
pandas-dev/pandas: Pandas 1.3.4.
Indexed in OpenAIRE,
2021
. https://zenodo.org/record/5574486. Accessed 19 October 2021.

27.

Braga-Neto
UM
,
Dougherty
E.
Error estimation for pattern recognition
. Hoboken, NJ: Wiley-IEEE Press,
2004
.

28.

Jacobs
JL
,
Bain
W
,
Naqvi
A
, et al.
Severe acute respiratory syndrome coronavirus 2 viremia is associated with coronavirus disease 2019 severity and predicts clinical outcomes.
Clin Infect Dis
2022
;
74
:
1525
33
.

29.

Zhao
J
,
Yuan
Q
,
Wang
H
, et al.
Antibody responses to SARS-CoV-2 in patients with novel coronavirus disease 2019.
Clin Infect Dis
2020
;
71
:
2027
34
.

30.

Marklund
E
,
Leach
S
,
Axelsson
H
, et al.
Serum-IgG responses to SARS-CoV-2 after mild and severe COVID-19 infection and analysis of IgG non-responders.
PLoS One
2020
;
15
:
e0241104
.

31.

Woodruff
MC
,
Ramonell
RP
,
Nguyen
DC
, et al.
Extrafollicular B cell responses correlate with neutralizing antibodies and morbidity in COVID-19.
Nat Immunol
2020
;
21
:
1506
16
.

32.

Rapsinski
GJ
,
Freeman
MC
,
Haidar
G
,
Belle
SH
,
Hasskamp
JH
,
Wheeler
SE.
Pediatric SARS-CoV-2 seroprevalence during mitigation procedures in Southwestern Pennsylvania.
J Clin Virol Plus
2021
;
1
:
100026
.

33.

Ou
X
,
Liu
Y
,
Lei
X
, et al.
Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV.
Nat Commun
2020
;
11
:
1620
.

34.

Han
H
,
Ma
Q
,
Li
C
, et al.
Profiling serum cytokines in COVID-19 patients reveals IL-6 and IL-10 are disease severity predictors.
Emerg Microbes Infect
2020
;
9
:
1123
30
.

35.

Bermejo-Martin
JF
,
González-Rivera
M
,
Almansa
R
, et al.
Viral RNA load in plasma is associated with critical illness and a dysregulated host response in COVID-19.
Crit Care
2020
;
24
:
691
.

36.

Magro
C
,
Mulvey
JJ
,
Berlin
D
, et al.
Complement associated microvascular injury and thrombosis in the pathogenesis of severe COVID-19 infection: a report of five cases.
Transl Res
2020
;
220
:
1
13
.

37.

Monsalvo
AC
,
Batalle
JP
,
Lopez
MF
, et al.
Severe pandemic 2009 H1N1 influenza disease due to pathogenic immune complexes.
Nat Med
2011
;
17
:
195
9
.

38.

Polack
FP
,
Teng
MN
,
Collins
PL
, et al.
A role for immune complexes in enhanced respiratory syncytial virus disease.
J Exp Med
2002
;
196
:
859
65
.

39.

Ali
YM
,
Ferrari
M
,
Lynch
NJ
, et al.
Lectin pathway mediates complement activation by SARS-CoV-2 proteins.
Front Immunol
2021
;
12
:1–8.

40.

Ankerhold
J
,
Giese
S
,
Kolb
P
, et al.
Circulating multimeric immune complexes drive immunopathology in COVID-19.
bioRxiv
[Preprint: not peer reviewed]. 9 September 2021. Available from: https://www.biorxiv.org/content/10.1101/2021.06.25.449893v4.

41.

Larsen
MD
,
Graaf de
EL
,
Sonneveld
ME
, et al.
Afucosylated IgG characterizes enveloped viral responses and correlates with COVID-19 severity.
Science
2021
;
371
:1–9.

42.

Aydillo
T
,
Rombauts
A
,
Stadlbauer
D
, et al.
Immunological imprinting of the antibody response in COVID-19 patients.
Nat Commun
2021
;
12
:
3781
.

43.

Anderson
EM
,
Goodwin
EC
,
Verma
A
, et al.
Seasonal human coronavirus antibodies are boosted upon SARS-CoV-2 infection but not associated with protection.
Cell
2021
;
184
:
1858
64.e10
.

44.

Lee
WS
,
Wheatley
AK
,
Kent
SJ
,
DeKosky
BJ.
Antibody-dependent enhancement and SARS-CoV-2 vaccines and therapies.
Nat Microbiol
2020
;
5
:
1185
91
.

45.

Gao
T
,
Hu
M
,
Zhang
X
, et al.
Highly pathogenic coronavirus N protein aggravates lung injury by MASP-2-mediated complement over-activation.
medRxiv
[Preprint: not peer reviewed]. 18 June 2020. Available from: https://www.medrxiv.org/content/10.1101/2020.03.29.20041962v3.

46.

Yan
B
,
Freiwald
T
,
Chauss
D
, et al.
SARS-CoV-2 drives JAK1/2-dependent local complement hyperactivation.
Sci Immunol
2021
;
6
:1–14.

47.

Yu
J
,
Yuan
X
,
Chen
H
,
Chaturvedi
S
,
Braunstein
EM
,
Brodsky
RA.
Direct activation of the alternative complement pathway by SARS-CoV-2 spike proteins is blocked by factor D inhibition.
Blood
.
2020
;
136
:
2080
9
.

48.

Harris
CL.
Expanding horizons in complement drug discovery: challenges and emerging strategies.
Semin Immunopathol
2018
;
40
:
125
40
.

Author notes

Potential conflicts of interest. U. B. N. reports grants from the National Science Foundation, and Texas A&M Institute of Data Science, and book royalties from Springer, un­related to the current work. J. M. reports grants from the National Institutes of Health (NIH), USAID, Gilead Sciences, and Janssen Pharmaceuticals; serves or has served as a consultant for Gilead Sciences, Accelevir Diagnostics, and Xi’an Yufan Biotechnologies; owns share options in Co-Crystal Pharmaceuticals and Infectious Diseases Connect; and is a shareholder of Abound Bio. His holdings in Co-Crystal Pharmaceuticals, Infectious Diseases Connect, and Abound Bio are unrelated to the current work. A. M. reports grants from the National Heart, Lung, and Blood Institute and the National Institute of Allergy and Infectious Diseases, NIH during the conduct of the study; B. J. M. reports grants from Bayer Pharmaceuticals, the Translational Breast Cancer Research Consortium, and the UPMC Learning While Doing Program, during the conduct of the study; consulting fees from Boehringer Ingelheim, the VeraMedica Institute; and the Patient-Centered Outcomes Research Institute for the ACCOMPLISH Trial, outside the submitted work. F. S. reports grants from the National Institute of General Medical Sciences, NIH, during the conduct of the study. B. J. C. M. reports grants from NIH and Gilead Sciences during the conduct of the study. D. B. reports leadership roles as president and board chair of Epistemix, director of the Magee Women’s Research Institute, and member of the Allegheny County Board of Health, unrelated to the current work, and is also a shareholder of Epistemix, unrelated to the current work. All other authors report no potential conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs licence (https://creativecommons.org/licenses/by-nc-nd/4.0/), which permits non-commercial reproduction and distribution of the work, in any medium, provided the original work is not altered or transformed in any way, and that the work is properly cited. For commercial re-use, please contact [email protected]