Limited cross-variant immune response from SARS-CoV-2 Omicron BA.2 in naïve but not previously infected outpatients
====================================================================================================================

* Hye Kyung Lee
* Ludwig Knabl
* Mary Walter
* Priscilla A. Furth
* Lothar Hennighausen

## ABSTRACT

Omicron is currently the dominant SARS-CoV-2 variant and several sublineages have emerged. Questions remain about the impact of previous SARS-CoV-2 exposure on cross-variant immune responses elicited by BA.2 infection compared to BA.1. Here we show that without previous history of COVID-19, BA.2 infection induces a reduced immune response against all variants of concern (VOC) compared to BA.1 infection. The absence of ACE2 binding in sera of previously naïve BA.1 and BA.2 patients indicates a lack of meaningful neutralization. In contrast, anti-spike antibody levels and neutralizing activity greatly increased in the BA.1 and BA.2 patients with a previous history of COVID-19. Transcriptome analyses of peripheral immune cells showed significant differences in immune response and specific antibody generation between BA.1 and BA.2 patients as well as significant differences in expression of specific immune genes. In summary, prior infection status significantly impacts the innate and adaptive immune response against VOC following BA.2 infection.

## INTRODUCTION

The highly transmissible SARS-CoV-2 Omicron (B.1.1.529) variant has replaced previous variants and is less susceptible to neutralizing antibodies elicited by vaccination or infection (Cele et al., 2022; Edara et al., 2022; Schmidt et al., 2022). Currently, BA.2 is the dominant Omicron sublineage and has replaced BA.1 (Viana et al., 2022). While BA.1 and BA.2 share 32 mutations, they differ in 28 and demonstrate different antibody neutralization profiles (Arora et al. 2022a). Omicron spike mutations negatively impact neutralizing activity in sera from previously vaccinated and previously infected persons as well as limited efficacy against monoclonal antibody treatment (Dejnirattisai et al., 2022; Hoffmann et al., 2022; Mannar et al., 2022; VanBlargan et al., 2022). Of note, vaccinated individuals with BA.1 infection develop measurable neutralizing antibody titers against BA.1 and BA.2 (Yu et al., 2022). Prior vaccination also broadens the immunological response against variants of concern (VOC) following Omicron infection whereas unvaccinated individuals show an antibody response more limited to the infecting strain (Lee et al., 2022b; Suryawanshi et al., 2022).

The ability of BA.2 infection to induce neutralizing antibodies in unvaccinated individuals, either without or with previous SARS-CoV-2 infection, is pending definition. Here we conducted a comparative investigation of the innate and humoral immune response elicited in 74 Omicron outpatients infected with BA.1 (Lee et al., 2022b) or BA.2 that were either naïve prior to Omicron infection or had prior exposure to other SARS-CoV-2 variants. We measured antibody titers and ACE2 binding inhibition as a proxy to neutralization. In addition, we investigated the immune transcriptome and the repertoire of antibody germline alleles in these cohorts.

## RESULTS

### Antibody response after Omicron subvariant infection

We measured antibody titers and neutralizing antibody responses in 74 outpatients infected with either Omicron BA.1 or BA.2 sublineages. Demographic and clinical information are provided in Table 1. While none of the patients had been vaccinated, 26 had been infected by an earlier variant. More than 20% of the BA.1 infected, but none of the BA.2 infected patients, had moderate or severe disease. Humoral immune response was measured between 7 and 21 days (median: days 13-14) post documented Omicron BA.1 or BA.2 infection. In previously uninfected and unvaccinated individuals, anti-Omicron BA.1 spike IgG titers were approximately 40-fold lower following BA.2 infection compared to BA.1 infection (Figure 1A; Table 1). IgG titers against the ancestral spike were almost 70-fold lower following BA.2 as compared to BA.1 infection. IgG antibody titers in individuals previously infected with earlier SARS-CoV-2 variants were at least 5-fold higher with no difference in titers between those infected with the BA.1 versus BA.2 variants (Figure 1A; Table 1). This highlights that previous infection with variants has a particular strong impact on the humoral immune response in BA.2 infected patients.

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/05/26/2022.04.07.22273565/T1)

Table 1. Characteristics of Omicron study population.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/26/2022.04.07.22273565/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/05/26/2022.04.07.22273565/F1)

Figure 1. Antibody analysis.
(A-B) Plasma IgG antibody binding the SARS-CoV-2 RBD (spike) from the ancestral and Omicron strains (A) as well as the SARS-CoV-2 RBD (spike) from different strains (B) in the Beta patients as a reference, the unvaccinated Omicron patients without or with prior infection experience. WA, authentic SARS-CoV-2 strain. (C-D) Neutralization response to virus spike protein of the ancestral and Omicron variants (A) as well as the SARS-CoV-2 RBD (spike) variants (B). *p*-value between two groups is from one-tailed Mann-Whitney t-test. Percentage of samples with zero neutralization activity for each group are indicated over the X axis. Asterisks indicate significance between groups. **p* < 0.05, ***p* < 0.01, \***|*p* < 0.001, \**\*|\**P* < 0.0001. Line at median, dotted line at 90%.

For reference, titers from individuals infected with the Omicron Beta variant (Knabl et al., 2022) are shown. These approximate those obtained in individuals previously infected with an earlier variant who then experienced BA.1 or BA.2 infection. The differential response in antibody generation following BA.1 and BA.2 infection was consistent across all SARS-CoV-2 variants tested (Figure 1B).

We next measured neutralization capacity using the angiotensin-converting enzyme 2 (ACE2) binding inhibition assay (Figure 1C). Neutralizing activities against the ancestral strain and Omicron BA.1 were similarly low in antigen naïve individuals infected with either BA.1 or BA.2. Similar to overall antibody titers, neutralizing activity was higher in those previously infected. Activity from BA.1 and BA.2 infected individuals with previous SARS-CoV-2 infection approximated that found for individuals infected with the Omicron Beta variant and was similar across all SARS-CoV-2 variants tested against (Figure 1D).

### Transcriptome response

Next, we investigated the impact of BA.1 and BA.2 infection on the immune transcriptome in Omicron patients with or without a previous infection history (Figure 2; Tables S2-3). Buffy coats were isolated within the first three days after validated Omicron infection and bulk RNA-seq was conducted with an average sequencing depth of 200 million reads per sample. First, we directly compared the transcriptomes of the BA.1 and BA.2 infected Omicron cohorts without prior infection. In a direct comparison of the two cohorts, expression of 38 and 15 genes was induced significantly in patients infected with BA.1 and BA.2, respectively (Figure 2A; Table S2). While genes activated upon BA.1 infection are enriched in neutrophil degranulation and innate immune system, the ones in activated upon BA.2 infection are categorized in complement pathway (Figure 2B; Table S2). Genes induced in BA.1 patients, such as C4BPA (Figure 2C), have also been reported to be regulated in COVID-19 patients (Gutmann et al., 2021; Jabeen et al., 2022). Of note, the complement genes C1QA and C1QB that show elevated expression in nonclassical monocytes from COVID-19 patients (Stephenson et al., 2021) are preferentially expressed in BA.2 patients.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/26/2022.04.07.22273565/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/05/26/2022.04.07.22273565/F2)

Figure 2. Immune transcriptomes following Omicron infection.
(A) Volcano plot of DEGs comparing BA.1 and BA.2 infection in unvaccinated Omicron patients with no prior infection at days 2-3 after infection. (B) Gene categories expressed at significantly higher levels in BA.1 or BA.2 infection group compared to the other. X-axis denotes statistical significance as measured by minus logarithm of FDR q-values. Y-axis ranked the terms by q values. (C) Bar plots with the fold change to mRNA levels between BA.1 and BA.2 groups. Red bars show genes that are higher in the BA.1 group compared to BA.2 and enriched in neutrophil degranulation and innate immune system, and blue bars indicate genes that are higher in the BA.2 group and enriched in complement pathways. (D) Volcano plot of DEGs comparing BA.2 infection in unvaccinated Omicron patients without or with no prior infection at days 2-3 after infection. (E) Gene categories expressed at significantly higher levels in BA.2 infection groups of no prior infection and prior infection. X-axis denotes statistical significance as measured by minus logarithm of FDR q-values. Y-axis ranked the terms by q values. (F) Bar plots with the fold change to mRNA levels in BA.2 groups of no prior infection and prior infection. Blue bars present genes that are higher in the no prior infection group compared to the prior infection group and enriched in innate immune responses.

Next, we compared the immune transcriptomes from BA.2 patients with and without prior infection and identified 342 and 145 significantly induced genes, respectively (Figure 2D; Table S3). GSEA analyses linked the induced genes in the ‘no prior infection’ cohort to innate immune responses, including interferon and inflammatory pathways (Figure 2E). Genes activated in the ‘no prior infection’ cohort (Figure 2F), such as CCL2 and USP18, are part of the interferon response in severe COVID-19 (Cao, 2021; Munnur et al., 2021).

### Antibody germline repertoire

The higher level of anti-spike antibodies and neutralization activity in BA.1 patients of the ‘no prior infection’ group compared to BA.2 infection, led us to dig deeper and interrogate the germline antibody gene usage. Specifically, we determined the range of immunoglobulin heavy chain (IGHV) gene usage elicited by BA.1 and BA.2 sublineage infection in the ‘no previous infection’ cohort (Figure 3). RNA-seq was conducted 13-16 days after the positive PCR test. The number of total immunoglobulin genes (heavy chain variable (IGHV), light chain variable (IGKV and IGLV) and T cell receptor alpha and beta variable (TRAV and TRBV)) revealed the use of a broad range of germlines in all cohorts without statistically significance (Figure 3A). However, we observed increased frequency of transcription of eight VH genes, including IGHV3-11, IGHV3-66, IGHV1-2 and IGHV4-61, in the BA.1 infection group compared to the BA.2 cohort (Figure 3B). These particular heavy chain variants have been identified in SARS-CoV-2 infected patients and COVID-19 vaccinated individuals that developed strongly neutralizing antibodies (Andreano et al., 2021; Andreano and Rappuoli, 2021; Collier et al., 2021; He et al., 2021; Lee et al., 2022a; Zhang et al., 2021). Expression of some of these genes was also activated in octogenarians which had received a single dose of the BNT162b mRNA vaccine 15 months after recovery from COVID-19 (Lee et al., 2022a).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/26/2022.04.07.22273565/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2022/05/26/2022.04.07.22273565/F3)

Figure 3. SARS-CoV-2-RBD-specific B-cell memory in Omicron patients.
(A) The number of total alleles that are generated at days 13-14. (B) Box plots show the frequency (%) of variable (V) gene usage for all paired heavy chains in BA.1 and BA.2 infection of the no prior infection groups. Ones identified in neutralizing antibodies to COVID-19 infection are marked in red. *p*-value between two groups is from one-tailed Mann-Whitney t-test. **p* < 0.05. Median, middle bar inside the box; IQR, 50% of the data; whiskers, 1.5 times the IQR.

## DISCUSSION

While the Omicron BA.1 and BA.2 sublineages share mutations they also carry unique mutations, 28 in BA.2, which can explain their differential response to vaccine-induced and therapeutic antibodies (Ai et al., 2022; Arora et al., 2022a; Arora et al., 2022b). Previous studies have also demonstrated that BA.1 infection of naïve persons yields a very limited antibody response to SARS-CoV-2 VOCs other than Omicron itself (Lee et al., 2022b; Rössler et al., 2022). However, prior vaccination or prior infection with older SARS-CoV-2 variants augment the immune response upon BA.1 breakthrough infection (Lee et al., 2022b; Rössler et al., 2022). Our study now reveals that in antigen naïve individuals the immunologic response following BA.2 infection is even lower than that elicited by BA.1 infection as indicated by significantly lower antibody titers against all VOCs. This also highlights that the set of 28 mutations unique to BA.2 contribute to the profound immune evasion.

While other studies have focused on the humoral response including neutralizing activity induced by Omicron infection (Suryawanshi et al., 2022), we have extended the investigations and explored the presence of germline immunoglobulin variants induced by BA.1 and BA.2 infection. RBD-targeting antibodies can block SARS-CoV-2 binding to ACE2 and neutralizing immunoglobulin G heavy-chain variable (IGHV) gene usage has been reported (Yuan et al., 2020). Our study identified several IGHV germline alleles in Omicron patients that have been identified in neutralizing antibodies in COVID-19 patients (Yan et al., 2021; Yuan et al., 2020).

In conclusion, our study demonstrates that Omicron BA.2 infection of unvaccinated individuals without prior COVID-19 does not elicit effective antibodies against BA.1 and non-Omicron variants. In previously infected individuals, however, BA.2 infection successfully induces immunity against other variants.

### Limitations of the study

There are several limitations to the current study. First, the study was conducted on volunteers from a specific geographical area, Tyrol (Austria). Secondly, limitation is the confinement of the study to a timeframe of three weeks following Omicron infection. Third, the number of the BA.2 patients with the prior infection experiment is limited. Finally, some of the data collected on breakthrough infections was reliant on retrospective chart review and not collected as part of a prospective study.

## Data Availability

RNA-seq data generated from this study will be deposited in the Gene Expression Omnibus (GEO) before publishing the manuscript. RNA-seq data of Omicron patients were obtained under GSE190747. 

## AUTHOR CONTRIBUTIONS

H.K.L., L.K. and L.H. designed the study. L.K recruited patients and collected material. H.K.L. analyzed RNA-seq data. M.W. conducted antibody assays. H.K.L., L.K., P.A.F. and L.H. analyzed data. L.H. supervised project. H.K.L., L.K., P.A.F. and L.H. wrote the paper. All authors read and approved the manuscript.

## DECLARATION OF INTERESTS

The authors declare not competing interests.

## INCLUSION AND DIVERSITY

We worked to ensure ethnic or other types of diversity in the recruitment of human subjects. We worked to ensure that the study questionnaires were prepared in an inclusive way. While citing references scientifically relevant for this work, we also actively worked to promote gender balance in our reference list. The author list of this paper includes contributors from the location where the research was conducted who participated in the data collection, design, analysis, and/or interpretation of the work.

## STAR METHODS

### KEY RESOURCES TABLE

View this table:
[Table2](http://medrxiv.org/content/early/2022/05/26/2022.04.07.22273565/T2)

## RESOURCE AVAILABILITY

### Lead Contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact Lothar Hennighausen (lotharh{at}nih.gov).

### Materials Availability

This study did not generate new unique reagents.

### Data and Code Availability

*   RNA-seq data generated from this study will be deposited in the Gene Expression Omnibus (GEO) before publishing the manuscript.

*   RNA-seq data of Omicron patients were obtained under GSE190747.

*   This paper does not report original code.

*   Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

## EXPERIMENTAL MODEL AND SUBJECT DETAILS

### Study population, study design and recruitment

Data from a total of 74 patients infected with either Omicron BA.1 or BA.2 were analyzed. 13 BA.2 patients and 10 new BA.1 patients were recruited for the study under informed consent 50 with no history of prior vaccination and prior infection with other SARS-CoV-2 variants and 24 unvaccinated patients with prior infections (Table 1). Recruitment and blood sample collection took place between December 2021 and March 2022. Day numbers for samples refers to the number of days following initial positive SARS-CoV-2 RT-PCR test. An additional 53 BA.1 patients from our previous study (Lee et al., 2022b) were included. Patient recruitment was performed by a medic who assessed clinical status including performance of an oxygen saturation (SpO2) test. The clinical spectrum of the patients’ SARS-CoV-2 symptoms were classified based on the National Institutes of Health (NIH) treatment guidelines ([https://files.covid19treatmentguidelines.nih.gov/guidelines/section/section_43.pdf](https://files.covid19treatmentguidelines.nih.gov/guidelines/section/section_43.pdf)). This study was approved (EK Nr: 1064/2021) by the Institutional Review Board (IRB) of the Office of Research Oversight/Regulatory Affairs, Medical University of Innsbruck, Austria, which is responsible for all human research studies conducted in the State of Tyrol (Austria). The investigators do not need to have an affiliation with the University of Innsbruck. Written informed consent was obtained from all subjects. Participant information was coded and anonymized. This study was determined to impose minimal risk on participants. All methods were carried out in accordance with relevant guidelines and regulations. All research has been have been performed in accordance with the [Declaration of Helsinki](https://www.wma.net/policies-post/wma-declaration-of-helsinki-ethical-principles-for-medical-research-involving-human-subjects/) ([https://www.wma.net/policies-post/wma-declaration-of-helsinki-ethical-principles-for-medical-research-involving-human-subjects/](https://www.wma.net/policies-post/wma-declaration-of-helsinki-ethical-principles-for-medical-research-involving-human-subjects/)). In addition, we followed the ‘[Sex and Gender Equity in Research – SAGER – guidelines](https://researchintegrityjournal.biomedcentral.com/articles/10.1186/s41073-016-0007-6)’ and included sex and gender considerations where relevant.

## METHOD DETAILS

### Antibody assay

Whole blood was collected by medical personnel after subjects had tested positive for SARS-CoV-2. Antibody containing sera were obtained by centrifuging EDTA blood samples for 10 min at 4,000g. End-point binding IgG antibody titers to various SARS-CoV-2–derived antigens were measured using the Meso Scale Discovery (MSD) platform. SARS-CoV-2 spike, nucleocapsid, Alpha, Beta, Gamma, Delta, and Omicron spike subdomains were assayed using the V-plex multispot COVID-19 serology kits (Panel 23 (IgG) Kit, K15567U). Plates were coated with the specific antigen on spots in the 96 well plate and the bound antibodies in the samples (1:50000 dilution) were then detected by anti-human IgG antibodies conjugated with the MSD SULPHO-TAG which is then read on the MSD instrument which measures the light emitted from the tag.

### ACE2 binding inhibition (Neutralization) ELISA

The V-PLEX COVID-19 ACE2 Neutralization kit (Meso Scale Discovery, Panel 23 (ACE2) Kit, K15570U) was used to quantitatively measure antibodies that block the binding of ACE2 to its cognate ligands (SARS-CoV-2 and variant spike subdomains). Plates were coated with the specific antigen on spots in the 96 well plate and the bound antibodies in the samples (1:10 dilution) were then detected by Human ACE2 protein conjugated with the MSD SULPHO-TAG which is then read on the MSD instrument which measures the light emitted from the tag.

### Extraction of the buffy coat and purification of RNA

Whole blood was collected, and total RNA was extracted from the buffy coat and purified using the Maxwell RSC simply RNA Blood Kit (Promega) according to the manufacturer’s instructions. The concentration and quality of RNA were assessed by an Agilent Bioanalyzer 2100 (Agilent Technologies, CA).

### mRNA sequencing (mRNA-seq) and data analysis

The Poly-A containing mRNA was purified by poly-T oligo hybridization from 1 mg of total RNA and cDNA was synthesized using SuperScript III (Invitrogen, MA). Libraries for sequencing were prepared according to the manufacturer’s instructions with TruSeq Stranded mRNA Library Prep Kit (Illumina, CA, RS-20020595) and paired-end sequencing was done with a NovaSeq 6000 instrument (Illumina) yielding 200-350 million reads per sample.

The raw data were subjected to QC analyses using the FastQC tool (version 0.11.9) ([https://www.bioinformatics.babraham.ac.uk/projects/fastqc/](https://www.bioinformatics.babraham.ac.uk/projects/fastqc/)). mRNA-seq read quality control was done using Trimmomatic (Bolger et al., 2014) (version 0.36) and STAR RNA-seq (Dobin et al., 2013) (version STAR 2.5.4a) using 150 bp paired-end mode was used to align the reads (hg19). HTSeq (Anders et al., 2015) (version 0.9.1) was to retrieve the raw counts and subsequently, Bioconductor package DESeq2 (Love et al., 2014) in R ([https://www.R-project.org/](https://www.R-project.org/)) was used to normalize the counts across samples (Zhao et al., 2021) and perform differential expression gene analysis. Additionally, the RUVSeq (Risso et al., 2014) package was applied to remove confounding factors. The data were pre-filtered keeping only genes with at least ten reads in total. The visualization was done using dplyr ([https://CRAN.R-project.org/package=dplyr](https://CRAN.R-project.org/package=dplyr)) and ggplot2 (Wickham, 2009). The genes with log2 fold change >1 or <-1 and adjusted p-value (pAdj) <0.05 corrected for multiple testing using the Benjamini-Hochberg method were considered significant and then conducted gene enrichment analysis (GSEA, [https://www.gsea-msigdb.org/gsea/msigdb](https://www.gsea-msigdb.org/gsea/msigdb)).

For T- or B-cell receptor repertoire sequencing analysis, trimmed fastq files from bulk RNA-seq were aligned against human V, D and J gene sequences using the default settings with MiXCR (Bolotin et al., 2017; Bolotin et al., 2015).

### Quantification and statistical analysis

Differential expression gene (DEG) identification used Bioconductor package DESeq2 in R. P-values were calculated using a paired, two-side Wilcoxon test and adjusted p-value (pAdj) corrected using the Benjamini–Hochberg method. The cut-off value for the false discovery rate was pAdj > 0.05. Genes with log2 fold change >1 or <-1, pAdj <0.05 and without 0 value from all sample were considered significant. For significance of each GSEA category, significantly regulated gene sets were evaluated with the Kolmogorov-Smirnov statistic. P-values of antibody between two groups were calculated using one-tailed Mann-Whitney t-test on GraphPad Prism software (version 9.0.0). A value of **P* < 0.05, ***P* < 0.01, \***|*P* < 0.001, \**\*|\**P* < 0.0001 was considered statistically significant.

## ACKNOWLEDGMENTS

This work was supported by the Intramural Research Program (IRP) of the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK).

Our gratitude goes to the participants who contributed to this study to advance our understanding of SARS-CoV-2 Omicron infection. We thank Yuhai Dai from the NIDDK clinical core for helping antibody assay. This work was utilized the computational resources of the NIH HPC Biowulf cluster ([http://hpc.nih.gov](http://hpc.nih.gov)). RNA-sequencing and single cell RNA-sequencing was conducted in the NIH Intramural Sequencing Center, NISC ([https://www.nisc.nih.gov/contact.htm](https://www.nisc.nih.gov/contact.htm))

## Footnotes

*   Lead author: LH: lotharh{at}nih.gov

*   Received April 7, 2022.
*   Revision received May 25, 2022.
*   Accepted May 26, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license

## References

1.  Ai, J., Wang, X., He, X., Zhao, X., Zhang, Y., Jiang, Y., Li, M., Cui, Y., Chen, Y., Qiao, R., et al. (2022). Antibody evasion of SARS-CoV-2 Omicron BA.1, BA.1.1, BA.2, and BA.3 sub-lineages. Cell Host Microbe 30, 1–7. doi:10.1016/j.chom.2022.05.001.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2022.05.001&link_type=DOI) 

2.  Anders, S., Pyl, P.T., and Huber, W. (2015). HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166–169. doi:10.1093/bioinformatics/btu638.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu638&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25260700&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000347832300003&link_type=ISI) 

3.  Andreano, E., Paciello, I., Piccini, G., Manganaro, N., Pileri, P., Hyseni, I., Leonardi, M., Pantano, E., Abbiento, V., Benincasa, L., et al. (2021). Hybrid immunity improves B cells and antibodies against SARS-CoV-2 variants. Nature 600, 530–535 doi:10.1038/s41586-021-04117-7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-04117-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34670266&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

4.  Andreano, E., and Rappuoli, R. (2021). Immunodominant antibody germlines in COVID-19. J Exp Med 218. doi:10.1084/jem.20210281.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1084/jem.20210281&link_type=DOI) 

5.  Arora, P., Zhang, L., Kruger, N., Rocha, C., Sidarovich, A., Schulz, S., Kempf, A., Graichen, L., Moldenhauer, A.S., Cossmann, A., et al. (2022a). SARS-CoV-2 Omicron sublineages show comparable cell entry but differential neutralization by therapeutic antibodies. Cell Host Microbe 30, 1–9. doi:10.1016/j.chom.2022.04.017.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2022.05.001&link_type=DOI) 

6.  Arora, P., Zhang, L., Rocha, C., Sidarovich, A., Kempf, A., Schulz, S., Cossmann, A., Manger, B., Baier, E., Tampe, B., et al. (2022b). Comparable neutralisation evasion of SARS-CoV-2 omicron subvariants BA.1, BA.2, and BA.3. Lancet Infect Dis. doi:10.1016/S1473-3099(22)00224-9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(22)00224-9&link_type=DOI) 

7.  Bolger, A.M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 doi:10.1093/bioinformatics/btu170.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu170&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24695404&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340049100004&link_type=ISI) 

8.  Bolotin, D.A., Poslavsky, S., Davydov, A.N., Frenkel, F.E., Fanchi, L., Zolotareva, O.I., Hemmers, S., Putintseva, E.V., Obraztsova, A.S., Shugay, M., et al. (2017). Antigen receptor repertoire profiling from RNA-seq data. Nat Biotechnol 35, 908–911. doi:10.1038/nbt.3979.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.3979&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

9.  Bolotin, D.A., Poslavsky, S., Mitrophanov, I., Shugay, M., Mamedov, I.Z., Putintseva, E.V., and Chudakov, D.M. (2015). MiXCR: software for comprehensive adaptive immunity profiling. Nat Methods 12, 380–381. doi:10.1038/nmeth.3364.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.3364&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25924071&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

10. Cao, X. (2021). ISG15 secretion exacerbates inflammation in SARS-CoV-2 infection. Nat Immunol 22, 1360–1362. doi:10.1038/s41590-021-01056-3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41590-021-01056-3&link_type=DOI) 

11. Cele, S., Jackson, L., Khoury, D.S., Khan, K., Moyo-Gwete, T., Tegally, H., San, J.E., Cromer, D., Scheepers, C., Amoako, D.G., et al. (2022). Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization. Nature 602, 654–656. doi:10.1038/s41586-021-04387-1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-04387-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35016196&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

12. Collier, D.A., Ferreira, I., Kotagiri, P., Datir, R.P., Lim, E.Y., Touizer, E., Meng, B., Abdullahi, A., Collaboration, C.-N.B.C.-., Elmer, A., et al. (2021). Age-related immune response heterogeneity to SARS-CoV-2 vaccine BNT162b2. Nature 596, 417–422. doi:10.1038/s41586-021-03739-1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-03739-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34192737&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

13. Dejnirattisai, W., Huo, J., Zhou, D., Zahradnik, J., Supasa, P., Liu, C., Duyvesteyn, H.M.E., Ginn, H.M., Mentzer, A.J., Tuekprakhon, A., et al. (2022). SARS-CoV-2 Omicron-B.1.1.529 leads to widespread escape from neutralizing antibody responses. Cell 185, 467–484 e415. doi:10.1016/j.cell.2021.12.046.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2021.12.046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35081335&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

14. Dobin, A., Davis, C.A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., Batut, P., Chaisson, M., and Gingeras, T.R. (2013). STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi:10.1093/bioinformatics/bts635.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts635&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23104886&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312654600003&link_type=ISI) 

15. Edara, V.V., Manning, K.E., Ellis, M., Lai, L., Moore, K.M., Foster, S.L., Floyd, K., Davis-Gardner, M.E., Mantus, G., Nyhoff, L.E., et al. (2022). mRNA-1273 and BNT162b2 mRNA vaccines have reduced neutralizing activity against the SARS-CoV-2 omicron variant. Cell Rep Med 3, 100529. doi:10.1016/j.xcrm.2022.100529.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.xcrm.2022.100529&link_type=DOI) 

16. Gutmann, C., Takov, K., Burnap, S.A., Singh, B., Ali, H., Theofilatos, K., Reed, E., Hasman, M., Nabeebaccus, A., Fish, M., et al. (2021). SARS-CoV-2 RNAemia and proteomic trajectories inform prognostication in COVID-19 patients admitted to intensive care. Nat Commun 12, 3406. doi:10.1038/s41467-021-23494-1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-23494-1&link_type=DOI) 

17. He, B., Liu, S., Wang, Y., Xu, M., Cai, W., Liu, J., Bai, W., Ye, S., Ma, Y., Hu, H., et al. (2021). Rapid isolation and immune profiling of SARS-CoV-2 specific memory B cell in convalescent COVID-19 patients via LIBRA-seq. Signal Transduct Target Ther 6, 195. doi:10.1038/s41392-021-00610-7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41392-021-00610-7&link_type=DOI) 

18. Hoffmann, M., Krüger, N., Schulz, S., Cossmann, A., Rocha, C., Kempf, A., Nehlmeier, I., Graichen, L., Moldenhauer, A.S., Winkler, M.S., et al. (2022). The Omicron variant is highly resistant against antibody-mediated neutralization: Implications for control of the COVID-19 pandemic. Cell 185, 447–456.e411. doi:10.1016/j.cell.2021.12.032.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2021.12.032&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35026151&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

19. Jabeen, A., Ahmad, N., and Raza, K. (2022). Global Gene Expression and Docking Profiling of COVID-19 Infection. Front Genet 13, 870836. doi:10.3389/fgene.2022.870836.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2022.870836&link_type=DOI) 

20. Knabl, L., Lee, H.K., Wieser, M., Mur, A., Zabernigg, A., Knabl, L., Rauch, S., Bock, M., Schumacher, J., Kaiser, N., et al. (2022). BNT162b2 vaccination enhances interferon-JAK-STAT-regulated antiviral programs in COVID-19 patients infected with the SARS-CoV-2 Beta variant. Communications Medicine 2, 17. doi:10.1038/s43856-022-00083-x.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s43856-022-00083-x&link_type=DOI) 

21. Lee, H.K., Knabl, L., Moliva, J.I., Knabl, L., Sr.., Werner, A.P., Boyoglu-Barnum, S., Kapferer, S., Pateter, B., Walter, M., Sullivan, N.J., et al. (2022a). mRNA vaccination in octogenarians 15 and 20 months after recovery from COVID-19 elicits robust immune and antibody responses that include Omicron. Cell Rep 39, 110680. doi:10.1016/j.celrep.2022.110680.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2022.110680&link_type=DOI) 

22. Lee, H.K., Knabl, L., Walter, M., Knabl, L.S., Dai, Y., Füßl, M., Caf, Y., Jeller, C., Knabl, P., Obermoser, M., et al. (2022b). Prior vaccination exceeds prior infection in eliciting innate and humoral immune responses in Omicron infected outpatients. Frontiers in immunology in press.
    
    

23. Love, M.I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15, 550. doi:10.1186/s13059-014-0550-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0550-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

24. Mannar, D., Saville, J.W., Zhu, X., Srivastava, S.S., Berezuk, A.M., Tuttle, K.S., Marquez, A.C., Sekirov, I., and Subramaniam, S. (2022). SARS-CoV-2 Omicron variant: Antibody evasion and cryo-EM structure of spike protein-ACE2 complex. Science 375, 760–764. doi:10.1126/science.abn7760.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abn7760&link_type=DOI) 

25. Munnur, D., Teo, Q., Eggermont, D., Lee, H.H.Y., Thery, F., Ho, J., van Leur, S.W., Ng, W.W.S., Siu, L.Y.L., Beling, A., et al. (2021). Altered ISGylation drives aberrant macrophage-dependent immune responses during SARS-CoV-2 infection. Nat Immunol 22, 1416–1427. doi:10.1038/s41590-021-01035-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41590-021-01035-8&link_type=DOI) 

26. Risso, D., Ngai, J., Speed, T.P., and Dudoit, S. (2014). Normalization of RNA-seq data using factor analysis of control genes or samples. Nat Biotechnol 32, 896–902. doi:10.1038/nbt.2931.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.2931&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25150836&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

27. Rössler, A., Knabl, L., von Laer, D., and Kimpel, J. (2022). Neutralization Profile after Recovery from SARS-CoV-2 Omicron Infection. N Engl J Med 386, 1764–1766. doi:10.1056/NEJMc2201607.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2201607&link_type=DOI) 

28. Schmidt, F., Muecksch, F., Weisblum, Y., Da Silva, J., Bednarski, E., Cho, A., Wang, Z., Gaebler, C., Caskey, M., Nussenzweig, M.C., et al. (2022). Plasma Neutralization of the SARS-CoV-2 Omicron Variant. N Engl J Med 386, 599–601. doi:10.1056/NEJMc2119641.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2119641&link_type=DOI) 

29. Stephenson, E., Reynolds, G., Botting, R.A., Calero-Nieto, F.J., Morgan, M.D., Tuong, Z.K., Bach, K., Sungnak, W., Worlock, K.B., Yoshida, M., et al. (2021). Single-cell multiomics analysis of the immune response in COVID-19. Nat Med 27, 904–916. doi:10.1038/s41591-021-01329-2.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-021-01329-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33879890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

30. Suryawanshi, R.K., Chen, I.P., Ma, T., Syed, A.M., Brazer, N., Saldhi, P., Simoneau, C.R., Ciling, A., Khalid, M.M., Sreekumar, B., et al. (2022). Limited cross-variant immunity from SARS-CoV-2 Omicron without vaccination. Nature. doi:10.1038/s41586-022-04865-0.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-022-04865-0&link_type=DOI) 

31. VanBlargan, L.A., Errico, J.M., Halfmann, P.J., Zost, S.J., Crowe, J.E., Jr.., Purcell, L.A., Kawaoka, Y., Corti, D., Fremont, D.H., and Diamond, M.S. (2022). An infectious SARS-CoV-2 B.1.1.529 Omicron virus escapes neutralization by therapeutic monoclonal antibodies. Nat Med 28, 490–495. doi:10.1038/s41591-021-01678-y.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-021-01678-y&link_type=DOI) 

32. Viana, R., Moyo, S., Amoako, D.G., Tegally, H., Scheepers, C., Althaus, C.L., Anyaneji, U.J., Bester, P.A., Boni, M.F., Chand, M., et al. (2022). Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa. Nature 603, 679–686. doi:10.1038/s41586-022-04411-y.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-022-04411-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35042229&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

33. Wickham, H. (2009). Ggplot2 : elegant graphics for data analysis (Springer).
    
    

34. Yan, Q., He, P., Huang, X., Luo, K., Zhang, Y., Yi, H., Wang, Q., Li, F., Hou, R., Fan, X., et al. (2021). Germline IGHV3-53-encoded RBD-targeting neutralizing antibodies are commonly present in the antibody repertoires of COVID-19 patients. Emerg Microbes Infect 10, 1097–1111. doi:10.1080/22221751.2021.1925594.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/22221751.2021.1925594&link_type=DOI) 

35. Yu, J., Collier, A.Y., Rowe, M., Mardas, F., Ventura, J.D., Wan, H., Miller, J., Powers, O., Chung, B., Siamatu, M., et al. (2022). Neutralization of the SARS-CoV-2 Omicron BA.1 and BA.2 Variants. N Engl J Med 386, 1579–1580. doi:10.1056/NEJMc2201849.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2201849&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35294809&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

36. Yuan, M., Liu, H., Wu, N.C., Lee, C.D., Zhu, X., Zhao, F., Huang, D., Yu, W., Hua, Y., Tien, H., et al. (2020). Structural basis of a shared antibody response to SARS-CoV-2. Science 369, 1119–1123. doi:10.1126/science.abd2321.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwNy8xMTE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDUvMjYvMjAyMi4wNC4wNy4yMjI3MzU2NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

37. Zhang, Q., Ju, B., Ge, J., Chan, J.F., Cheng, L., Wang, R., Huang, W., Fang, M., Chen, P., Zhou, B., et al. (2021). Potent and protective IGHV3-53/3-66 public antibodies and their shared escape mutant on the spike of SARS-CoV-2. Nat Commun 12, 4210. doi:10.1038/s41467-021-24514-w.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-24514-w&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34244522&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F26%2F2022.04.07.22273565.atom) 

38. Zhao, Y., Li, M.C., Konate, M.M., Chen, L., Das, B., Karlovich, C., Williams, P.M., Evrard, Y.A., Doroshow, J.H., and McShane, L.M. (2021). TPM, FPKM, or Normalized Counts? A Comparative Study of Quantification Measures for the Analysis of RNA-seq Data from the NCI Patient-Derived Models Repository. J Transl Med 19, 269. doi:10.1186/s12967-021-02936-w.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12967-021-02936-w&link_type=DOI)