Combinatorial analysis reveals highly coordinated early-stage immune reactions that predict later antiviral immunity in mild COVID-19 patients

Summary While immunopathology has been widely studied in patients with severe COVID-19, immune responses in non-hospitalized patients have remained largely elusive. We systematically analyze 484 peripheral cellular or soluble immune features in a longitudinal cohort of 63 mild and 15 hospitalized patients versus 14 asymptomatic and 26 household controls. We observe a transient increase of IP10/CXCL10 and interferon-β levels, coordinated responses of dominant SARS-CoV-2-specific CD4 and fewer CD8 T cells, and various antigen-presenting and antibody-secreting cells in mild patients within 3 days of PCR diagnosis. The frequency of key innate immune cells and their functional marker expression are impaired in hospitalized patients at day 1 of inclusion. T cell and dendritic cell responses at day 1 are highly predictive for SARS-CoV-2-specific antibody responses after 3 weeks in mild but not hospitalized patients. Our systematic analysis reveals a combinatorial picture and trajectory of various arms of the highly coordinated early-stage immune responses in mild COVID-19 patients.


In brief
In a parallel, longitudinal, prospective cohort based on a general population, using a systems-immunology analysis approach, Capelle et al. provide a comprehensive resource of multi-layered early immunological responses in mild COVID-19 patients. They reveal that a combination of highly coordinated earlystage immune responses is unique to mild patients.

INTRODUCTION
The immunopathology underlying severe COVID-19 has been thoroughly studied over the last 2 years, including antibody (Ab) responses, cellular immune subsets, cytokines, and chemokines that were linked to characteristics and outcome of the disease. [1][2][3][4][5][6] However, with few exceptions, 7 relatively little is known about the details of the immune response in mild patients (MP) and asymptomatic COVID-19 patients (ASP). Using profiling analyses of immune cell subsets, several studies have identified crucial alterations in hospitalized patients (HP) with severe symptoms versus HP with moderate disease, convalescent patients, and healthy controls. These studies, which mainly utilized flow cytometry or single-cell mRNA sequencing, have demonstrated a wide spectrum of abnormal immune responses to SARS-CoV-2. 3,8,9,10 However, ASP and/ or MP have only rarely been focused upon in these studies to draw conclusions. Therefore, it is unclear whether early protective immune signatures are identifiable in ASP or MP, and how such immune signatures might compare with HP and CoV-2 S inhibition % healthy people. Importantly, as reported cohort studies often lack a simultaneous analysis of the different aspects of multifaceted immune responses, it still remains unknown whether early-stage immune responses have any consequence for the evolution of other later immune reactions in MP. Several studies have included ASP and MP in cross-sectional analyses. For example, Ab responses and several cytokines have been analyzed in ASP versus symptomatic patients. 11 Also, SARS-CoV-2-specific and functional memory T cells have been detected in recovered ASP and MP 12 or in recovered patients with undefined disease severity. 13 Such cross-sectional studies were critical to identifying dysregulated immune factors associated with severe COVID-19. However, the isolated analysis of specific cellular immune subsets or cytokines and Ab responses alone will only allow for a partial understanding of the complexity of the early immune trajectories following infection. Furthermore, owing to different kinetics of immune responses among various patient groups, only a head-to-head comparison in a longitudinal, prospective study design can guarantee the comparability of observations. Because of partial cohort aggregation and non-harmonized sampling time points, this was insufficiently addressed in another longitudinal project. 7 In our longitudinal cohort with $220 samples characterized by a parallel and prospective study design, we sought to address the aforementioned open questions. To this end, we simultaneously analyzed 484 immune combinations resulting from 36 lineage and functional markers in three multicolor flowcytometry panels, 24 serological cytokine markers, Ab titers to SARS-CoV-2 spike (S), receptor-binding domain (RBD), N-terminal domain (NTD), and nucleocapsid (N), ACE2-binding inhibition to S and RBD as surrogate for Ab neutralization capacity, and SARS-CoV-2-specific T cells using both ex vivo T cell receptor b (TCR-b) repertoire sequencing and in vitro activation-induced marker (AIM) assays. Such a comprehensive, simultaneous, and integrated analysis in a longitudinal cohort using a systems-immunology strategy, as we and others have suggested in other contexts, 14,15 helps to draw a full picture of immune responses among MP.

RESULTS
Serological and whole blood count analysis distinguishing hospitalized from mild and asymptomatic COVID-19 patients We established the longitudinal Predi-COVID cohort 16 in Luxembourg during the first wave of the pandemic with the aim of gaining a systematic understanding of the early antiviral immune response across the full spectrum of COVID-19 disease phenotypes ( Figure 1A). All patients were included in Predi-COVID with a maximum delay of 3 days after a positive PCR test. According to available clinical data, patients were stratified into ASP (n = 14), mild to moderate (n = 63; referred to as MP in Figure 1A; for details, refer to STAR Methods), and HP (n = 15) subgroups for further analyses. All patients (n = 92) were sampled on the day of inclusion (D1) and 3 weeks after inclusion (D21). We also included control individuals (n = 26) from patients' households (HC), who were sampled on D1 and day 14 (D14). While the age was not different for ASP and MP as compared with HC (Table S1 and Figure 1B), HP were older than MP (median, $57 versus $38 years) and HC ( Figure 1B). Body mass index was not different among any of the analyzed groups (Table S1). In general, more males were included in the patient subgroups (between 57% and 69%), while only around 30% of HC were male (Table S1). No comorbidity information was available for HP and HC. For the other patient groups, the prevalence of comorbidities (asthma, chronic hematologic disease, obesity, and uncomplicated diabetes) was higher among MP than ASP (5%-8% in MP versus none in ASP) (Table S1).
A whole blood count analysis was used to further characterize COVID-19 patients on D1. We found no significant difference between ASP and MP in any of the tested 17 general blood count parameters (Table S1 and Figure S1). However, HP showed a remarkable difference compared with MP or ASP as demonstrated in the principal component analysis (PCA) plot ( Figure S1A). In line with other reports, 17 the frequency of lymphocytes was substantially decreased in HP versus the two other groups ( Figure S1B), while C-reactive protein (CRP) was elevated in HP only ( Figure S1C). Furthermore, compared with other groups, HP showed a significant but modest increase in the number or frequency of white blood cells, monocytes, granulocytes, and platelets, whereas a significant but modest decrease in red blood count and hematocrit was observed (Figures S1D-S1K). Considering the fact that the time from the onset of first symptoms to diagnosis might take longer in severe cases, i.e., in HP, we further analyzed only those HP with a shorter prodromal phase that was comparable with MP. To this end, we performed a sub-cohort analysis by selecting HP (n = 4) with at most a 7-day delay from the appearance of the first symptoms to inclusion to match the average time interval of all the MP (for details, refer to STAR Methods). Encouragingly, both the CRP levels and the percentages of lymphocytes were still significantly higher or lower, respectively, among HP versus MP on D1 (Figures S1L and S1M).  As expected, on D1 we did not observe a significant increase in immunoglobulin G (IgG) Ab titers to SARS-CoV-2 S, RBD, N, and NTD in the non-hospitalized groups versus HC ( Figures 1C-1F). In contrast, while only 48% and 43% of MP showed slightly increased IgG levels, 93% and 73% of HP already displayed significantly enhanced IgG titers to CoV-2 S and N antigens, respectively. Three weeks later, we observed a significant increase in IgG levels to all four antigens also in MP compared with HC together with a further enhancement of IgG levels in HP. The positivity rate for IgG Abs to CoV-2 S, RBD, NTD, or N reached up to 89% among MP and 100% among HP at D21. ASP had a lower positivity rate for IgG against S/RBD (67%) and N (75%) than MP at D21 (Figures 1C-1E). IgG Abs against NTD were in general much lower in ASP and MP, except for HP at D21 ( Figure 1F). Next, we tested the functional capacity of the induced Abs in a surrogate virus neutralization assay. In line with the serology findings on D1 ( Figures 1C and  1D), HP already showed blocking Abs that interfered with ACE2 binding to CoV-2 S or RBD at this early stage ( Figures 1G and 1H). On D21, MP also had developed a significant inhibitory capacity to block ACE2 binding to CoV-2 S or RBD versus HC. Similar to the correlation reported by others, 18 IgG Ab titers against both CoV-2 S and RBD were highly correlated (Spearman's r = 0.89 and 0.82 for S and RBD, respectively) with the inhibitory capacity of sera from all the patients (Figures 1I and 1J).
Early-stage highly coordinated innate and adaptive immune responses in mild COVID-19 As shown above, routine laboratory data as well as deep serologic profiling of SARS-CoV-2-specific Ab responses already distinguished HP from MP or ASP and HC. However, these analyses were not sufficient to further differentiate non-hospitalized clinical phenotypes. A meta-analysis has found that only 3% of the mean convalescent neutralizing Ab levels are necessary to predict protection from severe COVID-19, 19 but the same work has also indicated that those Abs are not sufficient to protect from severe cases.
Thus, we aimed to explore the full complexity of innate and adaptive cellular immune signatures that orchestrate the response to SARS-CoV-2 infection across the full spectrum of COVID-19 disease phenotypes. We systematically investigated 484 cellular immune subsets or combinations of various lineage and functional markers ( Figure 1A) by three staining panels using 18-color flow cytometry (for general gating strategy, see Figure S2; for cellular markers analyzed, refer to key resources table) in our longitudinal cohort. When compared with HC, ASP displayed no obvious change in all of the analyzed 484 immune profiles at D1 ( Figure S3A). Interestingly, at D21, ICOS + cells were the only significantly changed immune subset with a decrease in the frequency among CD8 T cells from ASP ( Figures S3B and S3C).
We next used PCA to show that 484 immune features were only able to partition HP at D1 from all other groups, but not MP at D1 from any HC (both D1 and D14) (Figure 2A). We then asked whether specific immune subsets were differentially present in MP versus HC at D1 ( Figure 2B). We observed differences in the frequency of several CD8 T cell subsets, such as Ki67 + , CD38 + , and HLA-DR + CD38 + , representing proliferating, activated, and antigen-specific responsive CD8 T cells, respectively that were significantly enhanced in MP ( Figures 2C-2E). The profile included an increase of both T-bet-dependent (Tbet + Ki67 + ) and -independent (EOMES + Ki67 + ) responsive CD8 T cells ( Figures 2F-2H). Also, the fraction of proliferating CD4 T cells, especially Th1-responsive (T-bet + Ki67 + ) cells, was already enhanced early on in MP on D1 ( Figures 2I, 2J, and S4A). In parallel, the frequency of antigen-presenting cells (APCs) and Ab-secreting cells, such as mature dendritic cells (HLA-DR + CD38 high DCs) and short-lived plasmablasts (CD27 + CD38 high ), was increased in MP ( Figures 2K-2M, S4B, and S4C). Notably, the frequency of activated CD38 + CD8 T cells, HLA-DR + CD38 + CD4 T cells, and mature DCs measured at D1 was highly predictive for the degree of the serological titers of anti-SARS-CoV-2 N IgG and other CoV-2 antigens at D21 among ASP and MP ( Figures 2N-2P). Since the magnitude of the responses of CD4 T cells, CD8 T cells, and mature DCs at D1 was all highly correlated with the same parameter (i.e., Ab levels) at D21, the response levels of these three subsets at D1 should correlate with each other. Such early responses of those immune subsets were highly coordinated only in ASP and MP. On the contrary, it is noteworthy that neither the frequency of activated subsets among CD4 and CD8 T cells nor of mature   DCs at D1 was significantly correlated with the even higher titers of CoV-2 N IgG ( Figures 2Q-2S) and other CoV-2 antigens in HP at D21. This finding indicates that the progression and deterioration of COVID-19 is averted only in the presence of a highly coordinated interplay of early innate and adaptive immune responses, which are strongly correlated with the subsequent production of Ab titers. On D21, MP were characterized by an enhanced cytotoxic CD8 T cell (GZMB + ) response, especially of terminally differentiated responsive CD8 T cells (CD45RO À Ki67 + ) ( Figures S3D-S3F). CD4 T cells also showed similar changes. The frequencies of CD4 T cells expressing CD57 or GZMB as well as of CD45RO and CD57 double-positive cells were also significantly enhanced in MP versus HC ( Figures S3G-S3I). Notably, the frequency of GZMB + CD4 cytotoxic T cells showed a trend to be elevated (p = 0.053, Kruskal-Wallis test including multiple-group correction) already on D1 ( Figure S3H). The CD57 expressing CD4 T cells detected on D21 appeared to be mainly cytotoxic effector cells, since the percentage of GZMB + CD57 + cells was also significantly enhanced among CD4 T cells ( Figure S3J). Although CD57 is known as a T cell senescence marker, CD57 + T cells, similar to the scenario of PD-1 + T cells, 20 were apparently still functional during the acute phase of COVID-19, thus likely contributing to sufficient control of the infection in MP.
Early-stage impaired innate immunity in hospitalized but not mild patients Until now, we parsed primarily early immune signatures in MP in relation to HC on D1 and D21. Yet the analysis of early cellular responses and later Ab responses showed a positive correlation only in MP and ASP, but not HP ( Figures 2N-2S). This finding prompted us to further compare between HP and other groups. Thus, we asked whether any additional early immune signatures observed in MP were significantly different from HP and determined the immune signatures that were significantly upregulated or downregulated in MP versus HP on D1 ( Figure 3A). As shown in the volcano plot, major differences were present primarily among innate immune cells, such as monocytes, DCs, and natural killer (NK) cells, and to a lesser extent also among B and T cells. Compared with HP, MP showed a much higher frequency ($40% in MP versus $10% in HP) of non-classical monocytes (ncMono, HLA-DR + CD38 À ). 21 The ncMono were not only higher in frequency among MP but also expressed higher levels of critical functional markers, such as CD86/CD80 double positivity ( Figures 3B-3D), PD-L1, and CD13 ( Figures S4D, S5A, and S5B). Similar to monocytes, the frequency of APCs such as plasmacytoid DCs (pDCs) and myeloid DCs (mDCs, also known as cDCs) was significantly higher in MP versus HP (Figures 3A and 3D-3G). It is noteworthy that the frequency of ncMono and mDCs ( Figures 3B and 3F) was also slightly lower in MP than in HC, indicating a disease severityrelated effect and further supporting the involvement of both cell types in early protective immune responses of COVID-19. Although mature DCs were higher in both MP and HP ( Figure 2L), the frequency of CD86 À CD80 + cells ( Figures 3G and 3H) and of CD13 + cells ( Figure S5C) among total DCs was decreased in HP only, thus indicating a reduction in phagocytic and antigen-presenting capacity of individual DCs. 22 In line with the notion of reduced APC functions, the downstream events of APC activation, the frequency of activated CD4 T cells (CD27 + ICOS + ), and the ICOS mean fluorescence intensity (MFI) among CD8 T cells were decreased only in HP but not MP ( Figures S5D and S5E). Furthermore, the frequency of NK cells was also significantly decreased only in HP, but not in MP versus HC (Figures 3I and 3J). In line with the overall compromised innate immune cell profile, critical senescence and exhaustion markers such as KLRG1 and PD-1 were enhanced among several NK subsets in HP only ( Figures 3K, S4E, S5F, and S5G).
In contrast to the compromised innate immune compartment, the expression levels of CD86 and the percentages of PD-L1 + cells among class-switched memory B cells were even substantially enhanced in HP versus both HC and MP at D1 (Figures S5H and S5I). Considering these results together with the high SARS-CoV-2-specific IgG titers, the ACE2-blocking capacity of patient serum, and the high frequency of plasmablasts (even higher than MP at both D1 and D21, Figure 2M), we concluded that Absecreting cells were not impaired in both MP and HP at D1, thus leading to a robust Ab response after 3 weeks.
We next sought to understand whether the impaired innate immune response in HP was paralleled by early deviated CD8 T cell profiles. Although the MFI of ICOS on total CD8 T cells was decreased ( Figure S5E), the frequency of ICOS + cells was unchanged ( Figure S3C) and the frequency of CD40L + and PD-1 + GZMB + cells among CD8 T cells was even significantly enhanced in HP on D1 ( Figures S5J and S5K). Furthermore, since the frequency of CD8 T cells expressing other key functional markers was not decreased in HP ( Figures 2C-2G), the functional antiviral capacity of CD8 T cells was most likely not compromised in HP versus MP. Overall, these data indicate that the major deficiencies observed in HP on D1 were impaired innate immune cells and their APC functions rather than adaptive immunity.
Due to the difference between the onset of first symptoms and the inclusion day among different groups (refer to the rationale described in the first part of this section regarding whole blood count analysis), we further performed a sub-cohort analysis by comparing HP with a shorter prodromal phase with all MP. We  Figure S6). Although one could not expect a statistical significance in all the comparisons with a lower case number, the change trend between the four HP and all MP ( Figures S6A-S6K) clearly matched the results displayed in the entire cohort. Notably, even with only four selected HP, the frequency of the ncMono expressing the key functional marker PD-L1 was still significantly lower than in MP at D1 ( Figure S6C). This also held true for an increase of the inhibitory marker KLRG1 in one of the NK subsets ( Figure S6K). To further substantiate our conclusions, we also compared innate immune responses of MP on D21 with those of HP on D1. All MP at D21 were in a later stage of their infection course than all HP at D1. Highly encouragingly, almost all the impaired critical immune subsets or functional markers (Figures 3 and  S5) in HP at D1, except for the frequency of one NK subset ( Figure 3J), remained significantly higher (or lower for inhibitory markers) in MP even at D21 versus HP at D1. In summary, our data firmly support that it is the early-stage response of innate immunity (including NK cells) that differentiates MP and HP during natural SARS-CoV-2 infection.

Early-stage temporary elevation of IP10 and IFN-b in mild COVID-19 patients
To gain further insight into the coordinated early immune response of COVID-19, we analyzed 24 different cytokines (refer to STAR Methods) at both D1 and D21 in sera of all groups. Interestingly, at D1, we observed increased levels of interferon-g (IFN-g)-inducible protein 10 (IP10/CXCL10) in HP and MP ( Figures 4A and 4B). Unexpectedly, a similar regulation was found for the type I interferon IFN-b, which was previously reported to be undetectable in severe COVID-19 patients at around 10 days after symptom onset. 23 Both MP and HP showed a significant increase of IFN-b compared with HC at D1 ( Figure 4C). While the levels of IP10 and IFN-b showed only a temporary increase among MP, declining to normal levels at D21, both IP10 and IFN-b levels remained elevated in HP after 3 weeks ( Figure 4C). These results point to a crucial and dynamic role of IP10 and IFN-b, which is tightly regulated during the early stage of protective immune responses in COVID-19 patients. This notion is also supported by the fact that levels of IP10 and IFN-b were significantly correlated with the frequency of mature DCs among all the analyzed patients at D1 (Figures 4D and 4E). Interestingly, the serological IP10 and IFN-b levels were significantly correlated with the magnitude of the CD8 T cell response among ASP and MP at D1 (Figures 4F and 4G). IFN-b levels were also correlated with the frequency of HLA-DR + CD38 + cells among CD4 T cells ( Figure 4H), indicating that the temporary IFN-b surge was indeed correlated with T cell responses. Such a positive correlation did not exist in HP at D1 (Figures S7A-S7C), indicating again the importance of early-stage orchestrated immune responses preventing severe illness. The CD8 responses and IFN-b levels in HP even showed a trend to be negatively correlated (r = À0.48, p = 0.08, Figure S7B). Furthermore, both IP10 and IFN-b levels among MP significantly correlated with the SARS-CoV-2 viral load at D1 (Figures 4I and 4J), as demonstrated by a negative correlation with the PCR Cq values. On D21, none of the 24 circulating immune analytes showed a significant change in MP versus HC at D14 ( Figure S7D). We also could not observe any significant change among ASP versus HC at either D1 or D21 (Figures S7E and S7F). The increase of IP10 in MP and HP appears to be independent of IFN-g, which was only elevated in HP at D1 (Figures 4B, 4K, and S7G) and remained high at D21, but still within the normal range for most of the patients ( Figure S7G).
Regarding other circulating soluble factors, we found a significant and substantial increase in plasma levels of eosinophil chemotactic protein (eotaxin-1/CCL11) and vascular endothelial growth factor A (VEGFA) only in HP at D1 (Figures 4L and 4M). The enhanced levels of interleukin-6 (IL-6) and the regulatory cytokine IL-10 in HP versus MP on D1 were still mostly seen within the normal range ( Figures S7H and S7I). The elevation of IL-10 levels was in line with an increased percentage of FOXP3 + regulatory T cells in both MP and HP at D1 ( Figure S7J). On D21, only IFN-b and IP10/CXCL10 remained significantly and substantially higher in HP versus MP ( Figure S7K). The T helper 2 cytokine IL-5, which was not different on D1, was modestly decreased in HP versus MP on D21 ( Figure S7L).

Early-stage dominant expansion of CD4 + SARS-CoV-2specific T cells in mild patients
To identify the T cell response on a broader scale, we performed TCR-b sequencing analysis among 45 MP versus 8 ASP and 21 Significantly increased or decreased cytokines (p % 0.05 and fold change R2) are marked in red and green, respectively. Due to its highly significant p value, IFN-b is also marked although its fold change was slightly less than 2. The gray dot represents the analytes showing a significant change but with values less than the reported normal levels even in HC. HC on D1 and D21. Aging has a strong impact on the TCR repertoire. 24 As expected, sample clonality, the inverted normalized diversity index, was significantly correlated with age ( Figure 5A). Since a decrease in TCR diversity was previously associated with aging and impaired immunity against influenza virus infection and other diseases, 25,26 we sought to compare the TCR diversity between different groups. The productive clonality of the sequenced TCR-b repertoire was increased in MP at D21 versus HC ( Figure 5B). At D21, only the use of one specific V gene (TCRBV06-07) was significantly under-represented in MP compared with HC ( Figure 5C). Notably, the SARS-CoV-2-specific T cell clonotypes were substantially expanded (about six times higher than in HC) among MP already at D1, as reflected by both clonal breadth and depth, 27 and maintained at D21 (Figures 5D and 5E). These results indicate a key functional role of early-responsive SARS-CoV-2-specific T cells in MP. Unexpectedly, inferred expanded CD4 SARS-CoV-2-specific T cells showed an average frequency about six times higher than that of inferred CD8 SARS-CoV-2-specific TCR clonotypes among MP at D1 (Figures 5F-5H). This finding was in line with a trend for increased frequency of GZMB + cells among CD4 Tconv cells, but not among CD8 T cells, in MP versus HC on D1 (Figures S3D and S3H). At D21, inferred CD4 SARS-CoV-2-specific TCR clonotypes continued to dominate over CD8 clonotypes in MP (Figures 5G and 5H). The expansion of inferred CD4 or CD8 SARS-CoV-2-specific T cells was highly correlated with the frequency of responsive ICOS + Ki67 + cells among total CD4 or CD8 T cells in both ASP and MP at D1 (Figures 5I and 5J). These data highlight a crucial role of early-expanding SARS-CoV-2-specific T cells, especially of the CD4 + cells, for subsequent coordinated antiviral immune responses. An earlier study has already reported SARS-CoV-2specific T cells early after diagnosis. 28 However, the authors have not correlated their findings with a severity stratification of COVID-19 patients nor have they performed further sub-analysis of CD4 and CD8 T cells. Therefore, our findings advance the understanding of the early antiviral responses of T cell subsets, particularly in MP.
To further consolidate our observations on inferred SARS-CoV-2-specific T cells based on sequencing approaches, 27 we used another more direct experimental approach to analyze virus-specific T cells in a cytokine-independent but viral-peptide-specific way. To this end, the AIM assay was selected to independently quantify SARS-CoV-2-specific T cell responses. The AIM assay has been successfully used by others to quantify SARS-CoV-2-specific CD4 T cells by detecting 4-1BB and OX40 double-positive 29,30 or CD40L and 4-1BB double-positive cells. 31 Encouragingly, following stimulation with a peptide pool selected from the wild-type SARS-CoV-2 whole-virus proteome, the detected SARS-CoV-2-specific responses among CD4 T cells, as reflected by the stimulation index (SI) (i.e., after excluding background effects; for details refer to STAR Methods) of both AIM combinations, were indeed higher in D1 samples of MP versus HC (Figures 6A-6D). In line with our other T cell datasets (Figures 2C-2J) and humoral immune results ( Figure 2M), the frequency of SARS-CoV-2-specific CD4 T cells, as quantified by both AIM combinations, was even higher in HP than in MP already at D1 (Figures 6B-6D).
SARS-CoV-2-specific CD8 responses were detected by calculating the portion of AIM (CD69 + 4-1BB/CD137 + ) 29,32 among total CD8 T cells. In line with our TCR sequencing data ( Figures 5F-5H), the overall restimulation response (SI) for CD8 T cells was much lower than for CD4 T cells in both groups ( Figures 6C-6E). Nevertheless, the SI for CD8 responses was still significantly higher in HP versus MP at D1 (Figures 6E and 6F). Contrary to CD4 T cells, we could not identify a significant difference between MP and HC for CD8 T cell responses ( Figures 6E and 6F). As shown above using the sequencing method, a significant difference was observed between MP and HC ( Figure 5H), and the CD8 T cell activation/proliferation response was highly correlated with the frequency of sequencing-derived SARS-CoV-2-specific T cells ( Figure 5J). Notably, the frequency ($5 3 10 À5 ) of SARS-CoV-2-specific cells among total CD8 T cells detected using the robust tetramer approach 12 was in a range similar to that inferred by the TCR sequencing approach here. This indicates that the ex vivo TCR sequencing-based approach might be even more sensitive than the in vitro AIM assay, at least in detecting antigen-specific CD8 T cells.
Excitingly, the SARS-CoV-2-specific CD4 T cells quantified by the AIM assays were highly correlated with the clonal breadth determined by the TCR-b repertoire sequencing method ( Figures 6G and 6H). This correlation did not exist for the quantified SARS-CoV-2-specific CD8 T cells ( Figure 6I). It is worth noting that our positive controls (those stimulated by the peptide pool of cytomegalovirus [CMV]) did not show any difference for any of the CD4 or CD8 AIM combinations between HC, MP, and HP ( Figures S8A and S8B). No difference in CMV peptide-stimulated samples between different participant groups demonstrated that the observed SARS-CoV-2-specific T cell responses were not due to the unspecific immune responses of the participants within a particular group to a type of common pre-exposed antigen, e.g., CMV. 33 In short, the direct antigen-dependent AIM assays have consolidated our TCR-sequence-inferred results, firmly revealing a dominant early-stage SARS-CoV-2-specific CD4 T cell response among MP. The AIM assays also showed that HP had no deficiency in generating both CD4 and CD8 SARS-CoV-2-specific T cell responses.

DISCUSSION
In May 2020, a mass PCR screening program was implemented on a population-wide level in Luxembourg, 34 which allowed us to obtain unique access to PCR-positive ASP and MP and to prospectively recruit them into our longitudinal study that was initiated simultaneously. 16 This endorsed us to access a rich resource to fully explore and understand all essential facets of the early-stage and dynamic immunological changes following recent SARS-CoV-2 infection in MP, using an unbiased, combinatorial, and prospective approach. So far, the immune response in COVID-19 patients has only been investigated in a few longitudinal cohort studies. 17,[35][36][37][38] These studies concentrated on different time windows and usually put their major focus on one or two selected immunological aspects, which makes it challenging to directly compare them with our multi-faceted analysis. In line with another recent longitudinal study from the United Kingdom (UK), 7 we observed enhanced early-stage CD8 T cells and plasmablast responses in MP. In contrast to that report, our current study provides not only information on the number and frequency of a wide spectrum of peripheral immune subsets but also on the functional status of individual immune cell types. We found both early-stage T-bet-dependent and -independent CD8 T cell responses among MP. Also distinct from the UK cohort, 7 we observed robust early-stage responses of CD4 T cells with a profoundly enhanced frequency of type I IFN-dependent Tbet + Ki67 + CD4 T cells. As another unique result of our study, many other innate immune cells expressing all key functional markers were observed to be intact early on in MP versus HC, but substantially compromised in HP. The notion that early coordinated DC and CD4/CD8 T cell responses have indeed a functional consequence in SARS-CoV-2-specific immunity was further supported by the correlation between mature DCs, CD4, and CD8 T cell responses at early stage (D1) with serum Ab responses 3 weeks later (D21), only among ASP/MP but not HP. Consistent with our own observation of early-stage overall CD4 and CD8 T cell responses, we found a substantial expansion of SARS-CoV-2-specific T cells, predominantly of CD4 T cells, as demonstrated by both TCR repertoire sequencing and AIM assays, in MP already at D1. A more dominant SARS-CoV-2-specific CD4 T cell response in MP might be attributable to pre-existing CD4 memory T cells that cross-reacted with SARS-CoV-2. 39,40 Since we observed enhanced DC and coordinated CD4 and CD8 T cell responses very early on among MP, at most 3 days after PCR diagnosis, the concept of bystander CD8 T cell responses 7 might need to be adapted. Our findings suggest that appropriate and highly coordinated early-stage DC and antigen-specific CD4 and CD8 T cell responses are predictive of the evolution of appropriate immune responses in MP individuals at a later stage. In addition, our longitudinal analysis firmly demonstrates that adaptive immunity was not compromised in HP. Since a highly coordinated early-stage adaptive and innate immune response was only observed in MP but not HP, such a connection between the two arms of immunity appears to be one of the key conditions for a long-term favorable outcome in MP. In this context, it is worth highlighting the strength of our simultaneous and combinatorial analysis strategy of the different immunological facets. By comparing various individual immune features in an isolated and one-dimensional way, various groups in our cohort could not be fully distinguished. For example, both MP and HP displayed enhanced T cell and Ab responses, indicating that one cannot separate them based only on those responses. When analyzing innate immunity, although HP showed impaired responses compared with MP, the innate immune parameters were comparable between MP and HC. As a consequence, none of the individual immune factors alone can fully distinguish the three groups in our cohort. Only a two-or even multi-dimensional view that integrates all layers of innate and adaptive immunity allows for such a differentiation of the three groups in one snapshot. There are other established examples in medicine using multi-dimensional analysis, such as the identification of the different naive and memory T cell subsets that can only be fully distinguished by simultaneously analyzing CD45RA/CD45RO and CCR7 in two dimensions. 41 Another important observation that is very much in line with the coordinated early-stage DC and antigen-specific CD4 and CD8 T cell responses was the strong early induction of the type I interferon IFN-b in MP. Notably, the early rise of IFN-b levels in MP was followed by a decline to normal levels 3 weeks later. However, such a contraction of enhanced IFN-b levels was not seen in HP. Our current results on early temporary induction of IFN-b levels in MP were further confirmed by a strong correlation between the frequency of mature DC, one of the main producers of type I IFN, 42 and the circulating IFN-b levels on D1. More excitingly, the frequency of activated and proliferating CD4 and CD8 subsets among total CD4 or CD8 T cells was highly correlated with IFN-b/IP10 levels in ASP and MP, but not in HP, on D1. The highly synchronized early-stage immune Article ll OPEN ACCESS responses further correlated with viral load in MP. Such a multilayered early-stage coordination of responses including multiple factors such as SARS-CoV-2 viral load, DC activation, CD4 and CD8 T cell stimulation, and relevant circulating cytokine (e.g., IFN-b/IP10) levels might decisively build up the basis for a later beneficial outcome in MP. In contrast to our results, a different study 43 has shown that another antiviral IFN, IFN-l1 (type III), was negatively correlated with viral load in severe patients, which, however, was apparently not sufficient to guarantee a beneficial clinical outcome. The inverted correlation observed by us and others suggests that the regulatory roles of antiviral IFNs might be very different, if not opposite, in MP versus HP. Our observations of an early time-dependent induction of IFN-b in MP are at least partially in line with studies on genetic and autoimmune defects leading to impaired type I IFN responses that were correlated with severe COVID-19. 23,44,45 Those findings have prompted the proposition of early and transient intervention with recombinant type I IFN as a treatment option in severe patients. 46 Parallel to IFN-b, we observed substantially enhanced early IP10/CXCL10 levels without other signs of systemic inflammation in MP. IP10 showed a dynamic change very similar to that of IFN-b in MP. The key difference in the cytokine and chemokine responses (IFN-b/IP10) observed between MP and HP is their elevation duration rather than the magnitude. IP10, previously known to exclusively bind to CXCR3, has recently been identified as a high-affinity agonist for the anti-inflammatory atypical chemokine scavenger receptor ACKR2/D6. 47,48 Thus, during the coordinated early anti-SARS-CoV-2 immune response in MP, IP10 might play a role in resolving inflammatory responses. The longer-lasting high levels of IP10 are not unique to severe COVID-19 but occur also in other infectious diseases, 49 including SARS. 50 Our observation in HP is in line with a cross-sectional comparison study, 11 where the IP10 level was enhanced in HP versus ASP. Our results of HP were also similar to a longitudinal cytokine-focused observation in HP, 51 although no MP were included in these analyses. Furthermore, a sustained elevation of IP10 was already reported in severe cases, 52 but only a few MP were included. Similar observations have also been reported 38,53 or reviewed elsewhere. 54,55 However, such a temporal early upregulation of IP10 in MP has never been convincingly reported so far.
Another crucial observation of our study was a very-earlystage (i.e., on D1) signature of reduced frequency and functional impairment of innate immune cells, such as ncMono, pDCs, mDCs, and NK cells, in HP versus MP. Both the frequency of ncMono and the expression of their key functional markers were significantly reduced in HP. Notably, although the frequency of mature DCs was induced early on in both MP and HP, their functional subsets among total DCs were substantially reduced only in HP. In addition, similar findings were also made in HP regarding the reduced frequency of several NK subsets, paralleled by a substantial enhancement in the expression of the inhibitory and terminal differentiation marker KLRG1 on NK cells. Thus, our data of impaired innate immune signatures in HP confirm previous findings of impaired innate immunity in severe or critically ill patients. 38,56,57 However, without including MP into the analyses, none of those previous studies has highlighted that critical differences in early-stage innate immune responses actually exist between non-HP (i.e., ASP and MP) and HP. Thus, our results provide strong support for an early-stage impairment of innate immunity being a unique clinical immunological feature in HP.
Since our analysis was only performed in peripheral blood, our observations might be affected by a potential redistribution of immune cell types between blood and inflammatory tissues. According to the work deposited by Kedzierska and colleagues, 58 the respiratory tract, relative to blood, was dominated by infiltrating neutrophils and a higher frequency of intermediate monocytes/macrophages and effector T cells in severe patients. In our study, we observed a lower frequency of ncMono, pDC, mDC, and NK subsets and substantial changes in their critical functional markers in blood of HP versus MP on D1. Therefore, most of our results on the frequency of these immune cell subsets should not be caused by a potential redistribution between blood and infected tissues. Our observation on the expression levels of functional markers is obviously independent of a potential redistribution of immune cells. Our findings about the frequency of ncMono in blood might have to be evaluated with more caution, as monocyte infiltration into the airways has been identified as an important driver of severe COVID-19. 36 In summary, most of the differences we observed between HP and MP are likely not due to a potential redistribution between blood and infected tissues.
Our data demonstrate that only the combination of different immunological facets can simultaneously distinguish MP, HP, and HC. Notably, based on a sufficiently powered sample size (63 MP) and on the prospective longitudinal nature of our study, we discovered the frequency of specific activated/proliferating CD4 and CD8 T cell subsets and mature DCs at an early disease stage as even better predictors of ensuing humoral responses 3 weeks later than early plasmablast responses in MP. This points to a critical role of DC activation in coordinating very early antigen-specific T cell and later Ab responses. The immune signatures identified in our study bear the potential to be extrapolated to predict protective immune responses in vaccinated people early on. In fact, our discoveries are in line with a recent report showing that mRNA vaccination induces rapid abundant antigen-specific CD4 T cell responses in SARS-CoV-2-naive participants following the first dose, 59 which phenocopies our findings regarding natural infection in MP, thus indicating that the discoveries of our study will have a more general impact on understanding and further dissecting SARS-CoV-2-related immunity.

Limitations of the study
We would like to point out that the unique combination of various early immune response features were only observed in MP rather than in HP and HC. However, so far we only have data supporting that these immunological features were correlated with mild symptoms in those patients. More mechanistic investigations are still required to figure out exactly which of the early immunological factors play a causal protective role in MP. Although this work critically focused on MP with a plausible number of participants (n = 63), the sample size of ASP and HP was moderate (n = 14 and 15). Recruitment of more ASP and HP might enhance the analysis power of this study in the conclusions related to those groups of patients. Furthermore, we have only analyzed two time points. Although with all the restrictions for obtaining ethical approval for a cohort of mostly non-critically ill patients in the early pandemic (April 2020), we should have sampled at more time points following the infection, thus allowing us to provide a more detailed dynamic analysis of evolving immune responses. Furthermore, for practical and logistical reasons, we could not perform repetitive sampling in the same individuals immediately after infection for a size of $120 participants recruited through a population-based program. In other human cohorts, this might be possible if the focus would be on a group of selected people who are PCR tested regularly for professional reasons (e.g., sports athletes), which was not within the scope of our current study. Another limitation is the partial comparability of time from symptom onset to inclusion between MP and HP.
In our cohort, we could not observe any distinguishable immune signature in the peripheral blood of ASP throughout all adaptive and innate immunity analyses. Although the sample size was relatively small in our ASP, the most plausible explanation is a more prominent role of a tissue-resident rather than a systemic response, as recently shown in a pediatric cohort. 60

STAR+METHODS
Detailed methods are provided in the online version of this paper and include the following: d This work does not report original code. d Any additional information required to re-analyze the data reported in this work is available from the lead contact upon request.

EXPERIMENTAL MODEL AND SUBJECT DETAILS
Human subjects and cohort design The Luxembourg National Research Ethics Committee (CNER) has given approval to this study (Predi-COVID) (202003/07). The trial has been registered at ClinicalTrials.gov (NCT04380987). All collections were performed with approval from relevant ethic organizations. Informed consent was obtained from each participant prior to collection. The blood/swab sampling was performed by nurses from Clinical and Epidemiological Investigation Centre (CIEC) of LIH.
Predi-COVID is a prospective longitudinal cohort study composed of individuals older than 18 years of age with a positive PCR test for SARS-CoV-2 in Luxembourg. Blood or swab samples were collected by a nurse at the latest 3 days post clinical PCR diagnosis (baseline, as day 1) at home for asymptomatic and mild participants. The mean time lag from the onset of first symptoms to the inclusion for mild patients was 6 days. For hospitalized patients, except for two of them (sampled 5 or 6 days post hospital arrival), the remaining 13 patients were all sampled at the latest 3 days after hospitalization. The mean time lag from the onset of first symptoms to the inclusion for mild patients was 12 days. A follow-up visit was organized 3 weeks (day 21, D21) later. Please note that the number of participants in the same category might be slightly different between baseline and the follow-up due to either practical or technical issues (refer to Figure 1A). The indicated number in each figure referred to the larger one among the two analyzed time points for each category. The group of mild COVID-19 patients also contained 11 patients with self-reported shortness-of-breath symptoms that could not be confirmed by a physician and therefore could not be classified as moderate patients following the NIH guideline.
The Predi-COVID-H sub-study is a prospective longitudinal cohort study composed of household members of a Predi-COVID participant as controls. Biological samples were collected at the same time as for the Predi-COVID participant sharing the house (baseline, as day 1, D1) and 2 weeks later, at day 14 (D14). The swab samples were collected at day 14 for the household controls. More details on the study design have been described elsewhere. 16 The time point at day 21 was chosen as the ideal follow-up time point in all COVID-19 patients recruited. Although it would have been ideal to do the follow-up of the household controls at day 21 as well, the ethical committee requested us to perform the follow-up earlier, i.e. after 14 days. The reason for that decision was to capture a potential COVID-19 infection in the household control members early enough and, in such a case, to refer them early to the medical healthcare system.

METHOD DETAILS
Blood sampling and PBMC isolation Samples were collected from confirmed SARS-COV-2 positive patients and household controls by trained nurses from the LIH-CIEC. Blood samples were collected in CAT, K 2 EDTA and CPT (all from BD, Erembodegem, Belgium) by the standard phlebotomy procedure. Blood samples were transported daily to centralized processing laboratory (IBBL) at ambient temperature.
CAT tubes were centrifuged for 10 min at 2000 x g, room temperature (RT). Serum upper layer was sterile aliquoted and stored at À80 C. Prior to centrifugation, 200 mL from the K 2 EDTA was transferred into 0.5 mL microcentrifuge tube for complete blood count (CBC) on ABX Micros CRP200 (Horiba, Japan). The K 2 EDTA tube was centrifuged for 20 min at 2000xg, 4 C. Plasma upper layer and buffy-coat were aliquoted and stored at À80 C. The CPT tubes were centrifuged for 20 min at 1800 x g, RT. The collected PBMCs were washed twice in Ca 2+ free PBS and counted using a Cellometer (Nexcelom, UK). Fresh PBMCs were partly used for direct flow cytometry and partly cryopreserved in CryoStor CS10 (Biolife solutions, USA) by controlled-rate freezing using Mr. Forsty (Nalgene, USA), followed by a long-term storage in liquid nitrogen.
Ex vivo multicolour flow-cytometry-based deep immunophenotyping analysis For each panel staining, 1 3 10 6 isolated fresh PBMCs, rather than frozen PBMCs, were used since cryopreservation affects several relevant markers as we demonstrated elsewhere. 61 The cells were resuspended in 50 ul of Brilliant stain buffer (BD, 563794) containing 2.5 ul of Fc blocking antibodies (BD, 564765) and incubated for 15 min. The suspension was then mixed with 50 ul of the respective 2x concentrated mastermix for the surface staining. The fluorochromes associated to the different markers are specified in key resources table. After 30 min of incubation in the dark at 4 C, the cells were washed three times with FCM buffer (flow cytometry [FCM] staining buffer, Ca 2+ -and Mg 2+ -free PBS + 2% heat-inactivated FBS) (5 min, 300xg). Following the final washing step, the stained PBMCs were fixed in 200 ul of 4% PFA (ThermoFisher Scientific, 28906) and incubated at RT for 30 min in the dark. After the PFA fixation, the PBMCs were washed once in FCM buffer (5 min, 400 x g) and resuspended in 200 ul fixation buffer from the True-Nuclear Transcription Factor Buffer Set (Biolegend, 424401). After 1 h of incubation in the dark at RT, the cells were centrifuged down (5 min, 400xg) and resuspened in 200 ul FCM buffer and left at 4 C overnight.
In the next morning the cells were resuspended in permeabilisation buffer (Biolegend, 424401) containing 2.5 ul of Fc blocking antibodies (BD, 564765) and incubated for 15 min at RT. After the intracellular blocking step, the cells were resuspended in 100 ul  (key resources table). Of note, different fluorochromes were used for some markers among different panels. After 30 min of incubation at RT in the dark, the cells were washed three times with permeabilisation buffer (5 min, 400 x g) and resuspended in 100 ul of FCM buffer to proceed to the acquisition on a BD LSRFortessa TM analyzer. To ensure a consistent acquisition of all the markers over the whole duration of the study, the application settings of the instrument were saved during the first acquisition and applied to all the following samples of the cohort as a part of our clinical research standard. The data was analyzed using FlowJo v10.5.6.

Activation induced marker (AIM) assay
The cryopreserved PBMCs were recovered in 10 mL of pre-warmed RPMI medium (Thermo Fisher Scientific, 21870084) [supplemented with 5% heat-inactivated human serum (Merck, H5667), 2 mM GlutaMAX TM (Thermo Fisher Scientific, 35050061), 50 U/mL penicillin and 50 mg streptomycin (Thermo Fisher Scientific, 15070063)] and centrifuged for 10 min, 200xg, RT. They were then resuspended in 2 mL of pre-warmed medium and rested overnight in a 24-well plate at 37 C, 5% CO2. The next day 1 3 10 6 of PBMCs per well were seeded in a round-bottom 96-well plate, first blocked for 15 min at 37 C, 5% CO2 with 0.5 ug/mL anti-CD40 antibodies (Miltenyi Biotec, 130-094-133) and then stimulated with 1 ug/mL of the specific peptide pool or H20 for 24 h at 37 C, 5% CO2 following manufacture's recommendation. The peptide pool we used was either PepTivatorâ SARS-CoV-2 (derived from the whole-virus proteome of wildtype SARS-CoV-2, premium grade, Miltenyi Biotec, 130-127-309) or PepTivatorâ CMV pp65 (premium grade, Miltenyi Biotec, 130-093-438). Of note, we used H2O rather than DMSO as a negative control because the peptide pools were reconstituted in H2O according to the manufacturer's recommendation. Following stimulation, the cells were stained using the master mixture of surface Abs provided in key resources table. The staining procedure was identical to that described in immunophenotyping section above. Considering the extremely-low frequency of antigen-specific T cells, we acquired up to 1.5E5 of gated lymphocytes for the AIM assay samples on a BD LSRFortessa TM analyzer followed by the analysis using FlowJo v10.5.6. The statistical analysis and visualization was then performed using Graphpad prism 9.0. The stimulation index was calculated as the ratio between the original percentage of the corresponding AIM combination stimulated by the peptide pool and background percentage in the negative control of the same sample. If the background percentage was zero, we used the minimal background values in the given group (household controls, mild or hospitalized patients at the given day) as the corresponding denominator. Serological detection of IgG against SARS-CoV-2 by MSD assay V-Plex COVID-19 Coronavirus Panel 1 serology kits from MSD (reference K15362U) were used to detect the presence of IgG antibodies to SARS-CoV-2-Spike (S), SARS-CoV-2 Nucleocapsid (N), SARS-CoV-2-S N-terminal domain (NTD) and SARS-CoV-2-S receptor binding domain (RBD) in diluted sera (1/500) according to the manufacturer's instructions. The plate was read on an MSD instrument, which measures the light emitted from the MSD SULFO-TAG. To determine the cutoff values for positivity for SARS-CoV-2, we measured 35 patients in another cohort from Central Hospital of Luxembourg (PCR positive and >15 days symptom onset) and negative sera before the pandemic from 2019 stored in Luxembourg National Laboratory (LNS). Of note, the hospitalized samples (both at day 1 and day 21 of inclusion) in this cohort were measured using the plates with the lot number different from the other groups and the corresponding positive thresholds of those plates were calculated accordingly. To guarantee the comparability of positive percentages between the two batches, the choice of the cutoffs aimed for a similar sensitivity and specificity between the two batches of plates (the AUC analysis was done by GraphPad Prism 9.0). In Figure 1, the statistical test between non-hospitalized groups was based on the signals. Only the positive percentages, rather than the signals of hospitalized samples were directly compared with that of other groups because abs from different lots were used.

Determination of neutralization antibody capacity by MSD assay
The assay analyzes the capacity of antibodies to inhibit the binding of labelled recombinant ACE2, the human receptor for SARS-CoV-2, to CoV-2 S or CoV-2 RBD in a multiplex high-throughput format. Multiplex assays for the detection of neutralizing antibodies against SARS-CoV-2 (SARS-CoV-2-Spike and SARS-CoV-2 S RBD) were done on patient sera using the MSD COVID-19 ACE2 Neutralization Kits from MSD (Panel 1 reference K15375U) according to the manufacturer's instructions. The samples were diluted 50 times for the neutralization assay. Data were recorded on a MESO QuickPlex SQ 120 instrument, which measures the light emitted from the MSD SULFO-TAG. Results were reported as inhibition percentage calculated using the equation below. % Inhibition was calculated using the following equation: (1-Average Sample ECL Signal/Average ECL signal of calibrator) *100. Determination of SARS-CoV-2 viral load by realtime-reverse transcription PCR (rRT-PCR) Oropharyngeal and nasopharyngeal swabs were collected from SARS-CoV-2 positive patients on day 1 and from household controls on day 14. Both swabs were discharged in a single tube containing Universal Transport Medium (Copan, Italy). Swabs were transported at 4 C to the processing laboratory, where they were homogenized by vortexing for 30 sec. After centrifugation for 1 min at 1000xg, 4 C, supernatant was aliquoted and stored at À80 C until further analysis.
Viral RNA was extracted using the QIAamp Viral RNA Mini Kit (52906, Qiagen, Venlo, the Netherlands) from 140 mL of swab medium according to the manufacturer protocol and eluted in 60 mL of elution buffer. Presence of SARS-CoV-2 viral RNA was assessed independently of initial clinical PCR diagnosis by analyzing the N gene using rRT-PCR (CDC, N1 target). The reaction was carried using TaqPath TM one-Step RT-qPCR Master Mix, CG (A15299, Life Technologies, Merelbeke, Belgium) with a 20 mL of final volume, 0.5 mM final concentration of both primers (Fwd: 5 0 -GAC CCC AAA ATC AGC GAA AT-3 0 ; Rev: 5 0 -TCT GGT TAC TGC CAG TTG AAT CTG-3 0 ) and 0.125 mM probe (5 0 -FAM-ACC CCG CAT TAC GTT TGG TGG ACC-BHQ1-3 0 ) and 5 mL of RNA. Real-time PCR cycling conditions were used as follows: reverse transcription for 15 min at 50 C and 2 min at 95 C followed by 45 PCR cycles at 95 C for 3 sec; 55 C for 30 sec. Human RNase P RNA levels were assessed using similar conditions (annealing/elongation step at 58 C; Fwd primer 5 0 -AGA TTT GGA CCT GCG AGC G-3 0 ; Rev primer: 5 0 -GAG CGG CTG TCT CCA CAA GT-3 0 ; Probe: 5 0 -FAM-TTC TGA CCT GAA GGC TCT GCG CG-BHQ1-3 0 , CDC) to monitor sample quality and the RNA extraction process. All PCRs were carried on CFX96 Touch real time instruments and results analyzed with the CFX Maestro software (BioRad, Temse, Belgium).
To quantify the viral load in swabs, RNA extracts were tested in duplicates with the rRT-PCR protocol described above, together with a 3-fold dilution curve of the EDX SARS-CoV-2 standard (COV019, BioRad) tested in triplicate on each PCR plate. Samples for which viral RNA concentration exceeded the upper range of the dilution curve were further retested in duplicates after dilutions of RNA extracts.

TCR repertoire analysis and SARS-CoV-2-specific T-cell inference
Immunosequencing of the CDR3 regions of human TCRb chains was performed using the ImmunoSEQâ Assay (Adaptive Biotechnologies, Seattle, WA). Extracted genomic DNA from $5 3 e6 cryopreserved PBMCs was amplified in a bias-controlled multiplex PCR, followed by high-throughput sequencing. Sequences were collapsed and filtered in order to identify and quantitate the absolute abundance of each unique TCRb CDR3 region for further analysis as previously described. 62 Clonality was defined as 1-normalized Shannon's Entropy and was calculated on productive rearrangements. Clonality values approaching 0 indicate a very even distribution of frequencies, whereas values approaching 1 indicate an increasingly asymmetric distribution in which a few clones are present at high frequencies. Clonal breadth and depth of SARS-CoV2-associated TCRb sequences were calculated as previously described, 27 using a set of sequences described elsewhere. 28 Briefly, breadth is calculated as the proportion of unique annotated SARS-CoV-2 specific rearrangements out of the total number of unique productive rearrangements, while depth accounts for the extent of the expansion of those clonal lineages in the repertoire.
We inferred whether TCR sequences are CD4 or CD8 T-cells by associating each sequence to a Class II or Class I HLA using a logistic regression classifier with L1 regularization in an independent dataset. The HLA class associations were validated using data from independent multiplexed antigen-specific T-cell receptor stimulation assay (MIRA) experiments 63 in which the sequence was directly observed in a Class I or Class II experiment. We find $90% consistency between associations based on MIRA and our inferred sequence-HLA associations. Results from these models were used to infer CD4 or CD8 labels for the sequences from the current study with high confidence. Clonal breadth and depth for these labeled sequences were estimated as described above.

Statistical analysis
Both PCA and volcano plots were visualized using GraphPad Prism 9.0. Correlation analysis was based on either Spearman or Pearson correlation as indicated in the corresponding figures. We only displayed the top-ranked highly correlated results if the correlation coefficient was ranked in the top positions among 484 correlations, calculated between antibody levels, cytokines or TCR breadth and any of the 484 subsets. The corresponding P-values from correlation analysis were based on a two-tailed analysis using GraphPad Prism 9.0. P-value from each scatter dot plots was determined by the Kruskal-Wallis (nonparametric) test and corrected using the Dunn's multiple comparisons test from GraphPad Prism 9.0. P-value in the AIM assay was calculated by nonpaired Student t test. In addition to each individual value measured from each participant, data from each group were presented as mean ± standard deviation (S.D.). ns, not significant, *p<=0.05, **p<=0.01 and ***p<=0.001.