Transient anti-interferon autoantibodies in the airways are associated with efficient recovery from COVID-19

FlowBEAT reveals distinct anti-SARS-CoV-2 and anti-type I IFN autoantibody responses across tissues and disease state

To quantify the full breadth of antibody responses, including isotype usage and antigen specificity across the COVID-19 disease state, we developed a multiparameter assay called “flowBEAT” (flow cytometry-based Bead assay to detect Antigen-specific antibody isoTypes). FlowBEAT is a modular flow cytometry-based technology that measures 176 antibody parameters per sample, including eight human antibody isotypes (IgG1, IgG2, IgG3, IgG4, IgA1, IgA2, IgE, IgM) against a panel of 22 host and viral antigens (Supplemental Table 1), including host type-I interferons (anti-IFN-α2a and anti-IFN-ω) and SARS-CoV-2 structural (anti-spike, -nucleocapsid, -membrane, and - envelope) and non-structural (anti-NSPs and -ORFs) proteins. We applied FlowBEAT to paired airway and blood samples from COVID-19 patients (Fig. 1 and Supplemental Table 1).

• Download figure
• Open in new tab

Figure 1. Multi-tissue validation of the flowBEAT assay and cohort summary

(A) Graphical overview of flowBEAT sampling from different tissues to reveal antibody isotype (color scheme 1) and antigen specificity (color scheme 2). (B) Dot plot showing the linear decrease of anti-RBD signal by flowBEAT serial dilution of seropositive NIH standards and monoclonal human IgG1 anti-RBD. (C) Graphical summary of study cohort spanning mild, moderate, and severe COVID-19 with multi-tissue collection and longitudinal follow-up. See cohort details in Tables S2-5. (D) Sunburst plots of flowBEAT showing the breadth of antibody response against SARS-CoV-2 and autoantigens in the airway and systemic samples, grouped by disease severity.

We first determined the specificity, sensitivity, and reproducibility of the flowBEAT antigen-coated beads. We established the linear range of the assay by generating replicate serial dilutions of mouse and human monoclonal antibodies, NIH COVID-19 human serology standard (24) (Fig. 1B), and pre-pandemic serum as negative controls (Fig. 1B and Supplemental Figs. 1-2). We determined the background signal of non-specific antibody binding using BSA-coated control beads, which were added to each tube during data collection and whose signals were used to establish the lower limits of detection of the assay (Methods).

• Download figure
• Open in new tab

Figure 2. Longitudinal kinetics of anti-SARS-CoV-2 antibody isotype response associated with disease severity and vaccination in the airway mucosa and blood

(A) Longitudinal plots showing isotype-specific anti-Spike signal (Y-axis) throughout infection and up to four doses of mRNA vaccination (X-axis) in both the nasal mucosa and blood of mild and moderate COVID-19 patients. Solid lines connect longitudinal samples from individual patients, and a Loess regression showing isotype induction and wane is plotted on top with standard error as a shaded area. (B) Bar plots show the maximum anti-Spike signal by individual patients (points, one per patient) and are grouped by severity class (fill texture) and antibody isotype (color). Statistics by Wilcoxon rank-sum test. (C) Tissue-specific anti-NSP signature. 31/36 donors show nasal IgA1 anti-NSP (NSP1, NSP2, NSP12, ORF3a, ORF3b, and ORF8) within the first two weeks. Of these, 16/36 progress to develop systemic anti-NSP IgG1 after the first two weeks. Pie charts detail positive-detected donors, with shading corresponding to peak antibody level. Dashed lines represent gates to select donors who generate anti-NSP IgG1. Bar charts of CRP and peak viral load, with donors grouped by anti-NSP IgG1 status, points colored by severity group. Statistics by Wilcoxon rank-sum test. (D) Longitudinal plots showing induction and waning of anti-NP antibodies. Grey points highlight suspected viral exposures, identified by an increase in anti-NP signal. Grey Loess regression indicates anti-Spike signal as in (A).

Using flowBEAT, we measured the breadth of antibody responses in the airway mucosa and matching blood of 129 donors spanning from mild (n = 37), moderate (n = 33), and severe life-threatening (n = 23) COVID-19, sampled longitudinally by repeated collections throughout infection and subsequent vaccination. We also collected samples from infection-naïve donors (n = 36), including vaccinated and pre-pandemic donors, for a total of 505 unique samples (205 airway mucosal samples and 300 systemic serum/plasma samples, Fig. 1C and Supplemental Tables 2-4). FlowBEAT revealed that the amount and diversity of antibody response increased with disease severity, represented as peak antibody expression during acute infection (Fig 1D). As the disease severity increases from mild to severe COVID-19, new antibody isotypes and specificities were generated against SARS-CoV-2 non-structural proteins (NSP and ORF) and host IFNs (Fig. 1D). Thus, flowBEAT distinguished the breadth of antibody response against SARS-CoV-2 and host IFNs across tissues and disease state.

Airway-specific antibody isotype and specificity against SARS-CoV-2 proteins distinguish mild, moderate, and severe disease

To better understand the breadth of antibody response associated with disease progression and recovery, we assessed the isotype and specificity of antibodies against SARS-CoV-2 proteins longitudinally and across disease states. We longitudinally sampled the nasal mucosa and blood starting at disease onset, determined as the date of polymerase-chain reaction (PCR)-confirmed infection, and through recovery and subsequent mRNA vaccination. In mild and moderate COVID-19 patients who recovered from the disease, IgG1, IgG3, and IgA1 were the predominant anti-spike (anti-S) antibody isotypes in both nasal mucosa and blood, reaching peak levels before 30 days post-onset (Fig. 2A). While most antibody isotypes were maintained for more than three months post-onset in blood, they quickly waned in the nasal mucosa after disease recovery, except for IgA1 which was still detectable in > 60% of mild/moderate-recovered donors three months post-onset (Fig. 2A). Interestingly, nasal IgA1 and IgE isotypes of anti-S and anti-nucleocapsid protein (anti-NP) pre-dated the detection of IgM by at least one week (Supplemental Fig. 3), likely representing pre-existing immunity to S and NP epitopes conserved between SARS-CoV-2 and common cold coronaviruses (25, 26).

We also found airway- and blood-specific antibody isotype and specificity associated with COVID-19 disease state. In the airway, the production of IgA2, often related to hyperinflammatory responses in the mucosa (27, 28), increased with disease severity and was highest within the endotracheal aspirates (ETA) of acutely-infected severe COVID-19 patients (Fig. 2B). In blood, while IgM was comparable across disease states, a sharp increase in the production of IgG1 and IgG4 anti-S associated with increased COVID-19 severity (Fig. 2B). High IgG3, which can suppress type-I IFN (29), exclusively distinguished patients with moderate to severe disease. While 17/28 of moderate cases of COVID-19 developed high systemic IgG3, only 4/24 of mild cases had detectable systemic IgG3 during acute infection (<100 days post-onset) (Fig. 2B).

Finally, we identified a nasal-specific IgA1 antibody signature against other SARS-CoV-2 structural and non-structural proteins (NSP and ORF), which emerged within two weeks of infection (Fig. 2C and Supplemental Fig. 4). Notably, only a subset of patients, primarily those with moderate but not mild disease, went on to develop an equivalent anti-NSP IgG1 signature in the blood after two weeks of infection. This progression from nasal IgA1 to systemic IgG1 anti-NSP was linked to an increased level of systemic inflammation at the peak of infection (Fig. 2C). Thus, the initial nasal-specific antiviral (anti-NSP/ORF) response is restricted to the airway site of infection and progresses to systemic response only in patients with increased disease severity.

In contrast to the relative stability of isotype usage, particularly IgA1 and IgG1, throughout SARS-CoV-2 infection and recovery, vaccination induced de novo IgG2 and IgG4 after two or more mRNA vaccine doses (Fig. 2A). Furthermore, the isotypes IgG3 and IgA2 anti-S, which we show were associated with disease severity (Fig. 2B), were not boosted by mRNA vaccination (Fig. 2A). Surprisingly, we found that nasal IgA1 and IgG1 anti-S (but not anti-NP) were boosted by distal intramuscular mRNA vaccination (Fig. 2D), suggesting a role for vaccination in supporting mucosal sterilizing immunity and preventing community transmission. Thus, flowBEAT revealed an anti-NSP/ORF signature that is restricted to the airway site of infection in mild-recovered individuals and identified IgA2, IgG3, and IgG4 as potential biomarkers of COVID-19 disease severity, while IgG2 was induced only by repeated vaccination.

Viral-induced anti-IFN-α autoantibodies in the nasal mucosa are associated with lesser symptoms, robust anti-SARS-CoV-2 antibodies, and full recovery

To gain a better insight into the longitudinal dynamics between anti-SARS-CoV-2 antibodies and autoantibodies against type-I IFN associated with full recovery, we longitudinally assessed 36 of the mild and moderate COVID-19 patients for whom we collected matched nasal and blood samples throughout infection and recovery. We profiled multiple samples from disease onset through recovery and subsequent mRNA vaccination and linked the data to symptomology and other clinical features (Supplemental Table 3). We first investigated whether anti-type I IFN autoantibodies (anti-IFN-α2a and anti-IFN-ω) could be detected locally in the airways throughout the COVID-19 disease progression. To confirm the specificity of the anti-IFN autoantibodies in our flowBEAT assay and to compare them to other published studies (16, 19, 30, 31), we coated our assay beads with the same IFN-α2a and IFN-ω proteins previously shown to identify anti-IFN autoantibodies in COVID-19 patients (30). To distinguish background binding from antigen-specific signal, we performed serial dilutions of positive samples and spiked BSA-coated beads in positive samples to calculate non-specific antibody binding to our IFN-coated beads (Supplemental Fig. 2). We were able to identify anti-IFN responders and non-responders in all our COVID-19 cohorts.

Strikingly, we discovered a robust anti-IFN-α2a autoantibody response induced by SARS-CoV-2 soon after infection, which peaked within the first two weeks post-onset in the nasal mucosa (nasal anti-IFN-α2a detected in 72% or 26/36 of mild and moderately-infected patients longitudinally assessed from disease onset) (Fig. 3A). These new-onset nasal anti-IFN-α2a autoantibodies were transient and waned as patients recovered from symptomatic disease (Fig. 3A). Consistent with the anti-SARS-CoV-2 antibody isotype response described above, nasal anti-IFN-α2a was dominated by the IgA1 isotype in mild and moderate cases of COVID-19, while IgA2, often associated with mucosal hyperinflammation (27, 28), was detected at lower levels (Fig. 3A). Notably, the initial anti-IFN-α response was limited to the nasal mucosa and rarely progressed to a systemic response in the blood. Only 13 of the 36 donors with matched blood and nasal swabs developed blood anti-IFN-α, which appeared later (after two weeks post-onset) (Fig. 3C) and was dominated by IgG1 (Fig. 3C and Supplemental Fig. 4). However, when anti-IFN-α progressed to the blood, it persisted for over three months post-disease onset (Fig. 3C) and was more frequently detected in the blood of donors with moderate (10/24) as compared to mild (4/28) disease (Fig. 3D). We observed a similar pattern of anti-IFN-ω in the same donors, but at lower levels (Supplemental Fig. 4). Finally, in contrast to COVID-19 disease, mRNA vaccination did not induce anti-IFN-α autoantibodies (Figs. 3A, C), indicating that anti-IFN-α autoantibodies are a specific response to SARS-CoV-2 infection.

• Download figure
• Open in new tab

Figure 3. Transiently induced nasal anti-IFN-α autoantibodies are associated with less symptomology and efficient recovery from COVID-19

(A) Longitudinal plot showing the post-onset induction of nasal (IgA1, IgA2) anti-IFN-α2a autoantibodies in mild and moderate COVID-19. Days post-onset is shown on the X-axis with longitudinal samples from single donors linked by solid lines with population trends (Loess regression) plotted as solid black lines with shaded standard error. Pie charts show the frequency and maximum signal intensity of total mucosal IgA anti-IFN-α2a. (B) The longitudinal decline of nasal IFN-α2a cytokine plotted as pg/mL of nasal swab supernatant. Longitudinal plots binned by week, showing a decline of nasal IFN-α2a cytokine and an increase in nasal anti-IFN-α2a autoantibodies. (C) Longitudinal plot showing post-onset induction of systemic (IgG1) anti-IFN-α2a autoantibodies in mild and moderate COVID-19. Axes and Loess regression as in (A). Pie charts show the frequency and maximum signal intensity of total systemic IgG1 anti-IFN-α2a. (D) Frequency summary (percent detected) of systemic IgG1 anti-IFN-α2a grouped by COVID-19 severity (numbers of positive donors over total donors). Chi-square test, p < 0.001. (E) Dot plot of nasal IgA1 anti-IFN-α2a (Y-axis) and systemic anti-IFN-α2a IgG1 signals (X-axis). Points represent the patient’s maximum signal within the first six weeks post-onset. Patient populations are gated (dashed lines) as non-producers (black), nasal IgA1-producers (red), and systemic IgG1-producers (green). Heatmap of nasal IgA1 and IgA2, and blood IgG1 anti-IFN-α2a, patient demographic data, and features of local immune response and symptoms. Rows correspond to individual patients, while columns indicate features as labeled. Pie charts show the frequency of “shortness of breath” (SOB) symptoms between red and green gates, as well as the maximum signal for nasal IFN-α2a cytokine and systemic C-reactive protein (CRP) for each patient in the gated populations. Wilcoxon rank-sum test, p < 0.01 or p < 0.05, as indicated. (F) Longitudinal plots showing nasal viral load per mL (qPCR) and tissue-specific anti-IFN-α2a autoantibody kinetics. Groups are divided by red and green gates, as in (C). The bar chart summarizes the peak viral load. (G) Correlation plots between anti-IFN-α2a (X-axis) and composite anti-non-structural proteins (anti-NSP/ORF) (Y-axis). IgA1 anti-SARS-CoV-2 NSP/ORF = anti-ORF3a, -NSP7. IgA2 anti-SARS-CoV-2 NSP/ORF = anti-ORF3a, -ORF8, -NSP1, -NSP2, ORF4 (E), ORF5 (M). IgG1 anti-SARS-CoV-2 NSP/ORF = anti-ORF3b, -NSP1. IgM anti-SARS-CoV-2 NSP/ORF = anti-ORF8, -NSP1, ORF4 (E). Pearson’s product moment correlation coefficient, p < 0.001.

We also measured host IFN-α2a (IFN-α) cytokine secretion longitudinally in the same matched nasal swab and serum samples. We detected nasal IFN-α at least one week before the peak of anti-IFN-α autoantibodies. While nasal IFN-α peaked as early as one week post-onset, anti-IFN-α autoantibodies continued to increase, reaching peak levels around two weeks post-onset when host IFN-α waned (Fig. 3B). We observed a similar pattern of IFN-α and anti-IFN-α dynamics in the blood (Supplemental Fig. 5). Thus, we show a coordinated balance between IFN-α secretion and subsequent anti-IFN-α production in the nasal mucosa and blood of mild to moderate COVID-19.

Next, we explored whether anti-IFN-α autoantibodies in the nasal mucosa or blood could be associated with symptomatology or other clinical features. When distinguishing between nasal and blood anti-IFN-α producers, we found no differences with regard to sex, age, or BMI (Fig. 3E). We did, however, identify a significant association between anti-IFN-α autoantibodies and disease features, including patient symptoms and biomarkers of systemic inflammation. Patients that developed systemic IgG1 anti-IFN-α had elevated C-reactive protein (CRP), a marker for systemic inflammation (Wilcoxon rank-sum test, p < 0.01) and reported persistent shortness of breath (SOB) and other relevant symptoms (Chi-square test, p < 0.01, Fig. 3E). Of note, SOB is used to indicate disease severity by NIH guidelines and may act as an indication of viral spread to the lower respiratory tract (32). In contrast to systemic IgG1 producers, nasal IgA1 anti-IFN-α autoantibody producers showed significantly lower CRP, SOB, and other symptoms related to disease severity (Fig. 3E). Interestingly, patients who produced IgA2 anti-IFN-α in the absence of IgA1 had elevated CRP levels and even worse clinical symptoms (Fig. 3E), suggesting that IgA2 anti-IFN-α may be induced under exacerbated mucosal hyperinflammatory conditions. Thus, nasal IgA1 anti-IFN-α, but not nasal IgA2 or systemic IgG1, may benefit mild to moderate patients when induced at the mucosal site of infection. Finally, patients that produced nasal IgA1 anti-IFN-α, but not systemic IgG1, showed the highest production of host IFN-α2a in the airways (Fig. 3E), suggesting that anti-IFN-α autoantibodies may be induced as a local feedback loop by the host IFN-α response.

To determine whether the viral-induced anti-IFN-α response was correlated with viral load, we assessed the same nasal swabs for SARS-CoV-2 copy numbers. We found no significant differences in viral load (maximum viral copy numbers/mL by qPCR) between patients that produced either nasal or systemic anti-IFN-α (Fig. 3F). However, we found that anti-IFN-α was secreted days after the peak viral load (Fig. 3F), confirming that nasal anti-IFN-α was induced post-viral-exposure and not pre-existing in our mild to moderate patients. Since nasal anti-IFN-α autoantibodies were associated with efficient disease recovery, we assessed their impact on anti-SARS-CoV-2 antibodies. We found a positive linear correlation between anti-IFN-α and anti-SARS-CoV-2 immunity (Fig. 3G), particularly against non-structural proteins (NSP/ORF), which are known to induce host-cell apoptosis and immunomodulatory functions (6, 33). We found nasal-specific correlations between anti-IFN-α and anti-NSP7/ORF3a and an equivalent blood-specific correlation between anti-IFN-α and anti-NSP1/ORF3b (Fig. 3G). Thus, transient but robust production of nasal IgA1 anti-IFN-α was induced after SARS-CoV-2 infection in patients with mild and moderate COVID-19 and was associated with favorable outcomes, including fewer symptoms and robust anti-SARS-CoV2 immunity.

Hospitalized patients show elevated anti-IFN-α in the airways associated with uncontrolled IFN-α, higher viral load, and lower anti-SARS-CoV-2 immunity

To determine whether the coordinated balance between anti-IFN-α and IFN-α we identified in mild/moderate-recovered patients is disrupted in severe COVID-19, we performed flowBEAT on endotracheal aspirates (ETA) and peripheral blood from 23 hospitalized and unvaccinated patients admitted to the intensive care unit (ICU) with life-threatening (severe) COVID-19. We detected anti-IFN-α autoantibodies in the airway mucosa (ETA) of 95% of severe patients, while only 43% (10/23) of them had detectable anti-IFN-α in the blood (Fig. 4A). We observed a similar pattern with anti-IFN-ω in the same donors, but at lower levels (Supplemental Fig. 4). In contrast to mild and moderate outpatients, severe patients in the ICU sustained high levels of IFN-α secretion in the ETA that positively correlated with the amount of anti-IFN-α autoantibodies (Fig. 4A), revealing a dysregulated IFN-α hyperinflammatory response at the airway site of infection. Consistent with sustained hyperinflammation in the airways, hospitalized patients produced higher levels of the pro-inflammatory IgA2 anti-IFN-α, unlike mild and moderate patients, which predominantly produced the IgA1 isotype (Fig. 4B). Similarly, the anti-SARS-CoV-2 response (anti-S) and germline autoreactive VH4-34 antibodies, which we previously showed to be increased in the blood of hospitalized patients (34), also switched from the IgA1 to the pro-inflammatory IgA2 isotype in the airways (Fig. 4B).

• Download figure
• Open in new tab

Fig. 4: Severe COVID-19 patients show elevated airway anti-IFN-α associated with hyperinflammation and elevated viral load

(A) Pie charts of airway (ETA) and systemic anti-IFN-α2a autoantibody frequency and signal intensity in severe hospitalized COVID-19. The dot plot shows a positive correlation between the ETA’s anti-IFN-α2a autoantibodies (IgA1/2) and IFN-α2a cytokine. Pearson’s product moment correlation coefficient, p < 0.01. (B) The bar chart compares the ratio of IgA2:IgA1 isotype usage in the airways of mild/moderate (nasal swabs) and severe (ETA) patients. Wilcoxon rank-sum text, p < 0.001. (C) Dot plot of individual hospitalized patients showing systemic IgG(1-4) and airway IgA(1-2) anti-IFN-α2a. Patient populations are gated (dashed lines) as airway IgA-producers (red) and non-producers (green). The bar chart shows viral copies/mL in the ETA of each gated population. (D) Heatmap of ETA IgA and IgG anti-IFN-α2a, along with ETA IFN cytokine and other hyper-inflammatory markers correlated to anti-IFN-α2a (Pearson correlation, p<0.05). Rows correspond to individual patients, while columns indicate features as labeled. (E) Correlation plots between anti-IFN-α2a (X-axis) and composite anti-non-structural proteins (anti-NSP/ORF) (Y-axis). Colored dots indicate severe COVID-19 and hollow dots indicate mild/moderate COVID-19, as presented in Fig. 3G. Anti-SARS-CoV-2 NSP/ORF programs shown as in Fig. 3G. Pearson’s product moment correlation coefficient, p < 0.05.

Interestingly, airway anti-IFN-α autoantibodies in severe COVID-19 patients were associated with sustained viral load (Wilcoxon rank-sum test, p < 0.01, Fig. 4C), suggesting an ongoing antiviral immune response. Indeed, airway anti-IFN-α correlated positively with increased inflammatory cytokines along with elevated signals of myeloid recruitment and activation, including IFN-γ, IFN-λ1, I-TAC, IP-10, MCP-1, MIP-1β, TRAIL, and SDF-1a (Fig. 4D). However, while anti-IFN-α in mild and moderate patients correlated with enhanced anti-SARS-CoV-2 immunity (i.e., increased anti-NSP, - ORF, -M, and -E in the airway and blood, Fig. 3F), this positive correlation no longer persisted in the airways of severe patients (Fig. 4E). Thus, in hospitalized patients, higher levels of airway IgA2 anti-IFN-α in the absence of IgA1 or the presence of high systemic IgG1 correlated with higher viral load, hyperinflammation, lower anti-SARS-CoV-2 immunity and a more severe disease course.

Study participants and ethical approvals

The study included 129 participants enrolled through the University of California San Francisco (UCSF)’s COVID-19 Host Immune Response Pathogenesis (CHIRP) study, Emory University’s PROATECT study previously described (1), and the pre-pandemic and healthy control samples were provided by the Clinical and Translational Discovery Core at Children’s Healthcare of Atlanta (CHOA) and Emory University under IRB 00089506. All participants were assigned a new, sequential ID in the format “ERxxx.” IDs are known only to the study investigators as per Emory IRB 00003368.

CHIRP study: eligible participants were identified by confirmed SARS-CoV-2 RT-PCR nasal swab test results and had matched plasma and nasal swab biospecimens collected between 4/2020 – 11/2022 (Supplemental Tables 2, 3). All participants provided written informed consent. Longitudinal biospecimens were collected at weeks 0, 1, 3, 10, and 24 from baseline visits. A subset of participants also provided additional samples post-vaccination (2 weeks post-dose #1, 2 weeks post-dose #2, 8 weeks post-dose #1 and 2 weeks post-dose #4). Participants completed a detailed questionnaire at each study visit, including demographic, social, medical history, concomitant medication use, and symptom severity reporting as previously described (43) (Supplemental Table 2). The date of symptom onset (or lack of symptoms, if asymptomatic disease) was recorded at the baseline visit. COVID-19 symptom severity (mild versus moderate) was defined based on symptom burden (total number of reported symptom categories) and duration (weeks of symptoms) (Supplemental Table 2). Clinical labs, including SARS-CoV-2 viral load (nasal swab), complete blood count with differential, comprehensive metabolic panel, erythrocyte sedimentation rate (ESR), and high sensitivity C-reactive protein (hs-CRP) were collected at each visit, as well as whole blood for processing for plasma and serum. This study was approved by the UCSF and Emory University Committee on Human Subjects Research (UCSF IRB 20-30588; Emory IRB 00003368. All participants provided written informed consent, and all biospecimens used for research were de-identified.

Emory PROATECT study: eligible participants were identified by confirmed SARS-CoV-2 RT-PCR test results and enrolled as previously described (1) (Supplemental Table 3). Briefly, 23 patients with severe, life-threatening COVID-19 were recruited from the ICU of Emory University, Emory St. Joseph’s, Emory Decatur, and Emory Midtown Hospitals between 6/2020 – 2/2021. We also recruited 24 outpatients with mild and moderate COVID-19 in the Emory Acute Respiratory Clinic. Outpatients were classified by disease severity (mild or moderate) following NIH guidelines (32) using symptom data collected by a healthcare provider (Supplemental Table 3). Infection-naïve patients were recruited for pre- and post-vaccination antibody assessments in blood and nasal swabs. Informed consent was obtained for all study participants. Emory’s IRB approved all studies under protocol numbers 00003368, 00058507, 00057983, and 00058271.

Sample processing

Plasma and serum: samples were collected for all study participants, regardless of COVID-19 diagnosis or severity. To obtain serum, whole peripheral blood was collected and allowed to clot, and supernatant serum was isolated and cryopreserved. For plasma collection, whole blood was collected in EDTA tubes, and plasma was isolated by centrifugation at 400 × g for 10 minutes at 4°C followed by removal of the supernatant. To further precipitate platelets, contaminating cells, or remaining debris, the isolated plasma and serum were centrifuged at 4000 × g for 10 minutes at 4°C. Samples were distributed in aliquots and cryopreserved at -80C.

Nasal swabs: samples were collected from both nostrils and transferred to 3 mL of viral preservation media. Swabs (in media) were vortexed to liberate collected antibodies into media. Samples were further mechanically dissociated using a syringe and 25 Ga needle to disrupt aggregates of mucins and other respiratory secretions. Samples were centrifuged at 2000 × g to precipitate non-soluble aggregates. Processing of samples from SARS-CoV-2 positive donors was performed within a biosafety cabinet following BSL2+ protocols approved by the Emory University Institutional Biosafety Committee and Biosafety Officer.

Endotracheal aspirates (ETA): samples were collected by healthcare providers in the ICU from either endotracheal intubation or tracheostomy (Supplemental Table 5). ETA samples were transferred to a USDA-approved BSL3 containment facility at Emory for further processing. To disrupt mucin aggregates, ETA was mixed 1:1 with a 50 mM EDTA solution (final concentration 25 mM EDTA) in custom RPMI-1640 media deficient in biotin, L-glutamine, phenol red, riboflavin, and sodium bicarbonate (defRPMI-1640) and dissociated using a syringe as previously described (1). Samples were centrifuged at 250 × g, and supernatants were cryopreserved at -80°C. Before removal from BSL3 containment, ETA supernatants were UV-inactivated following our published protocols (44).

Antigen preparation and conjugation

All SARS-CoV-2 and auto-antigens were obtained from commercial sources, except for SARS-CoV-2 non-structural antigens, which were produced in collaboration with ChemPartner (Supplemental Table 1). All antigens, including a purified bovine serum albumin (BSA) control antigen, were biotin-conjugated using EZ-link NHS-biotin reagents at a 20-fold molar excess. Antigens were loaded onto commercially sourced streptavidin-coated beads at an empirically determined saturation point in 50 µL of PBS + 0.05% BSA + 0.2% Tween-20 (PBS-BSA) and allowed to bind for 20 minutes on ice before the excess (uncaptured) antigen was removed by washing with additional PBS-BSA.

Fluorescent detection antibodies

Fluorescently conjugated monoclonal antibodies were obtained from commercial sources for the specific detection of eight human antibody isotypes and subclasses (IgG1, IgG2, IgG3, IgG4, IgA1, IgA2, IgE, and IgM). IgG1-, IgA2-, IgE- and IgM-specific antibodies were purchased with fluorophores pre-conjugated by Southern Biotech, while IgG2-, IgG3-, IgG4- and IgA1-specific antibodies were conjugated using the Lightning-Link fluorescent conjugation kits (Abcam) following the optimal manufacturer protocol. We used mouse monoclonal anti-SARS-CoV-2 and anti-interferon antibodies to determine assay sensitivity and specificity, and a secondary goat polyclonal anti-mouse (PE or AF488 fluorophore) to reveal monoclonal anti-antigen captured on each bead species.

FlowBEAT calibration

We experimentally determined the antigen-bead saturation point for RBD, S1, S2, IFN-α2a, and IFN-ω antigens. All were found to have an identical saturation point of 2 pmol per 50,000 beads in 50 µL of PBS-BSA. We then confirmed the capacity of saturated beads to capture antigen-specific antibodies by flow cytometry. To confirm the antigen specificity signal for each mouse monoclonal antibody, we performed serial dilutions of antigen-specific antibodies co-incubated with antigen-saturated beads in PBS-BSA and revealed by a secondary stain with fluorescent anti-mouse. Each antigen-specific antibody had a maximum signal output at least 1000-fold greater than the unstained control, which decreases linearly with serial dilution.

FlowBEAT assays

After thawing, samples were centrifuged at 2,000 × g for 10 minutes at 4°C to precipitate any remaining aggregates. Serum or plasma samples were diluted to a final concentration of 1:125 with PBS-BSA in a reaction mix containing beads coated with antigen at an empirically determined saturation point. Airway samples were diluted to a final concentration of 1:2. Samples and beads were incubated on ice in the dark for 30 minutes, mixing at least every 15 minutes to prevent beads from settling. Beads were washed with PBS-BSA and resuspended before fluorescent antibody staining. Antibody staining was performed by mixing equal volumes of resuspended beads and a master mix containing eight fluorescently conjugated monoclonal antibodies specific to IgG1, IgG2, IgG3, IgG4, IgA1, IgA2, IgE and IgM. The optimal concentration of each antibody species was empirically determined by titration. The reaction was incubated on ice in the dark for 30 minutes, mixing at least every 15 minutes. Stained beads were washed with PBS-BSA and resuspended in 1X PBS. Samples were fixed in 1:6 final dilution of FACS Lysing solution (BD Biosciences), which we showed inactivates SARS-CoV-2 (44). While fixing, samples were kept in the dark at room temperature for 20 minutes, mixing at least every 10 minutes. A final wash with PBS-BSA was performed, and beads were resuspended in PBS-BSA. All samples were analyzed within 24 hours of staining using a 5-laser Cytek Aurora flow cytometer. All data collection occurred on the same instrument, and the cytometer was calibrated with a consistent lot of Cytek SpectroFlo QC beads before each data collection. To standardize each flowBEAT batch, we established a daily standard plasma sample with strong reactivity to SARS-CoV-2 spike subunits. The standard positive and a sample-free control were included in all flowBEAT assays to confirm the assay’s performance and control for potential variations.

Cytokine quantification

The levels of human cytokines in serum, nasal swabs, and ETA were measured using a Meso Scale Discovery U-plex custom multiplex assay kit, as we previously described (1). We measured an additional ten analytes (IFN-α2a, IFN-γ, IL-15, IL-29/IFN-λ1, IP-10, MCP-1, MDC, MIP-3β, SDF-1a, and TRAIL) following the manufacturer’s protocol (Meso Scale Discovery) in plasma, nasal swabs, and ETA, which were prepared following our established UVC-inactivation protocol (44).

Computational analysis

Flow Cytometry Standard (FCS) files were gated and analyzed in FlowJo version 10 software (BD Biosciences) to identify median fluorescence intensity (MFI) values for each antigen-specific antibody isotype fluorescence (example gating strategy in Supplementary Fig. 1). Exported MFI values were background-corrected using a linear regression correlating background antigen-specific signal (as determined in pre-pandemic donors) and control BSA signal. Figures were generated using R version 4.2.2 using the packages “ggplot2”, “ggcyto,” and “flowcore.” The code will be available upon publication on the Ghosn Lab’s GitHub page: github.com/Ghosn-Lab.