Impact of BNT162b first vaccination on the immune transcriptome of elderly patients infected with the B.1.351 SARS-CoV-2 variant

Fast-spreading variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) energize the COVID-19 pandemic. The B.1.351 variant carrying the escape mutation E484K in the receptor binding domain is of particular concern due to reduced immunological protection following vaccination. Protection can manifest as early as 10 days following immunization with full protection two weeks following the second dose, but the course is not well-characterized for variants. Here, we investigated the immune transcriptome of six elderly individuals (average age 82 yr.) from an old people’s home, who contracted B.1.351, with four having received the first dose of BNT162b eight to 11 days prior to the onset of COVID-19 symptoms. The patients were hospitalized and received dexamethasone treatment. Immune transcriptomes were established from PBMCs approximately 10 and 35 days after the onset of COVID-19 symptomology. RNA-seq revealed a more intensive immune response in vaccinated patients as compared to unvaccinated ones. Specifically, transcription factors linked to the JAK/STAT pathway, interferon stimulated genes, and genes associated with innate antiviral immunity and COVID-19-SARS-CoV-2 infection were highly enriched in vaccinated patients. This rendered the transcriptomes of the older vaccinated group significantly different than older unvaccinated individuals infected at the same institution and more similar to the immune response of younger unvaccinated individuals (age range 48–62) following B.1.351 infection. All individuals in this study whether vaccinated or not were hospitalized due to B.1.351 infection and one vaccinated patient died illustrating that although an enhanced immune response was documented infection it was insufficient to protect from disease. This highlights the need for maintaining physical distancing and prevention measures throughout the time course of vaccination in older adults.


Introduction
The SARS-CoV-2 variant B.1.351 1 has multiple changes in the immunodominant spike protein that facilitates viral cell entry via the angiotensin-converting enzyme-2 (ACE2) receptor. The receptor-binding domain (RBD) mutation E484K provides tighter ACE2 binding and widespread escape from monoclonal antibody neutralization 2-5 , a concern for a limited protection by some vaccines 6 , and recent data from Qatar demonstrate a BNT162b vaccine 7 effectiveness of 75% with the B.1.351 variant 7 . Vaccine response as well as the immunological response to SARS-CoV-2 infection are time-dependent, with first a rise and then a fall to normal levels of immune-related gene transcription 8 .

Differences in sampling times post vaccination and/or infection strongly influence results
derived. Limited data is available on the comparative response of the immune transcriptome between infected individuals who are vaccinated versus unvaccinated following the first dose. Since the immune response declines with age, it is most critical to explore this question in the older age group that is most affected by the COVID-19 pandemic. Here, we conducted a study on elderly patients in an old-people's home that experienced an outbreak of the SARS-CoV-2 B.1.351 variant. A unique aspect of this single setting study is that the infected individuals were all exposed within the same for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 14, 2021. ; https://doi.org/10.1101/2021.05.11.21256862 doi: medRxiv preprint timeframe and vaccinated individuals were all immunized within the same timeframe, controlling for environmental and kinetic variables.

Study design
To better understand the extent of protection of elderly individuals from COVID-19 after receiving the first dose of BNT162b vaccine, we investigated a unique cohort of hospitalized patients from an old people's home that experienced an outbreak of the

Immune transcriptome
Blood was drawn at ~10 and 30 days after the development of COVID-19 symptoms, RNA was prepared from buffy coats followed by RNA-seq analyses. The principal component analysis (PCA) showed a separation between RNA-seq samples from non-COVID controls 8 and the first and second time points from COVID-19 patients on the first principal component (PC1) (Fig. 1b) This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Elevated immune response in patients infected after the first vaccination
Next, we focused on the immune transcriptome response in six hospitalized elderly patients (average age 82 yr.) from the same old people's home (Supplementary Table 1 Table 1). The PCA plot shows a separation between the vaccinated and unvaccinated groups on PC1 (Fig. 2a), which is further strengthened by hierarchical clustering of the 181 significantly differentially expressed genes between the six patients (  Table 3). Specifically, the highly ranked ETV7, STAT1 and STAT2 ( Fig. 2f) are key transcription factors executing interferon responses and are top predictors for ACE2 expression in human airway epithelium 14,[16][17][18][19][20] . Notably, expression of IFITM3, a gene with a polymorphism that has been linked to the severity of COVID-19 21 , is highly induced in vaccinated patients ( Supplementary Fig. 1a). While three out of the four vaccinated patients had a highly activated immune transcriptome, patient S-02 had a blunted immune response, more similar to that seen in the non-vaccinated patients for use under a CC0 license.
This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Temporal progression of the immune response
From our cohort of the four vaccinated patients, one died of COVID-19 and three were discharged from the hospital after a single stay of 10-14 days (Fig. 1a). From the two nonvaccinated patients, one was discharged after a single ten day stay and one was readmitted. To dig deeper into the temporal immune response after discharge from the  Table 4). Expression of specific ISGs (e.g., IFI35, IFIT5, IFITM3, STAT1, STAT2, OAS2) was sharply reduced (Fig. 3b), approaching levels seen in non-COVID samples.

Immune responses at different ages
Having analyzed the immune transcriptome from the elderly hospitalized population (average age 82 yr.) we addressed potential age differences and included RNA-seq transcriptomes from three hospitalized patients all below 62 years of age (average age 56 yr.) (Fig. 1a). The two groups were separated in the PCA plot (  Table 5).
These results support the concept that age might be a defining factor in the immune response in COVID-19 patients.

Discussion
COVID-19 is a rapidly evolving field in which real-world data, especially those concerning the now prevalent variants of concern (VOC), is critical for everyday management of the COVID-19 disease by public health officials, government agencies and institutions. There for use under a CC0 license.
This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  Our study illustrates that whole transcriptome study from buffy coat is feasible and yields actionable data when performed in well-controlled clinical settings. It demonstrates a technique that is feasible to expand and reproduce across settings. While our study was limited to infection ~11 days following first vaccination, extensions could compare the immune transcriptome after the second vaccination and investigate individuals who may have waning antibody response ten or more months following either infection and/or full immunization. As a caveat, it can be challenging to reliably compare RNA-seq data generated on different platforms and the disease severity at the time of the analysis greatly impacts the transcriptome 10 .
While there was a measurable immune transcriptome response to first vaccination, it was insufficient to prevent significant clinical disease and protective measures against for use under a CC0 license.
This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 14, 2021. ; https://doi.org/10.1101/2021.05.11.21256862 doi: medRxiv preprint spread should be fully employed in individuals who have received only one dose.
Moreover, SARS-CoV-2 variants of concern (VOC) such as B.1.351 studied here demonstrate resistance to neutralization by antibodies generated from some current vaccines 11 . Further research in vaccinology will continue to address mechanisms to improve the protective immunological response to SARS-CoV-2 and VOC through boosters and vaccines engineered specifically for VOC. Serial immune transcriptome studies should be included in addition to antibody tests for a fuller understanding of the spectrum of immune response in real-world situations as more immunized individuals encounter SARS-CoV-2 VOC.

Limitations
The findings in this report are subject to at least three limitations. We have investigated the transcriptome from patients infected with B.1.351 but not with other variants. Our study focused on elderly (average age 82 yr.) vaccinated and unvaccinated patients with limited comparison to younger individuals (48-62 yr). We have investigated the effect of the BNT162b vaccine but not of any other vaccine. The study group was small; however, the parallel timing of vaccination, infection and treatment in the elderly group enabled the study to be well controlled for environmental factors.

SARS-CoV-2 virus sequencing
RNA was extracted from patient's blood using a Maxwell RSC simply RNA Blood purification kit according to the manufacturer's instructions (Promega, USA). Library preparation and sequencing was performed as described 26 . In short, cDNA was obtained by using reverse transcriptase with random priming. Following cDNA synthesis, primers based on sequences from the ARTICnetwork were used to generate 400 bp amplicons in two different PCR pools. After merging of pools and amplification, libraries were constructed using QIASeq FX DNA Library UDI Kit following the manufacturer's instructions (Qiagen GmbH, North Rhine-Westphalia, Germany).
Sequencing was done with Illumina NextSeq® 500/550 using 149-bp paired-end reads with 10-bp indices (Illumina, California, USA). Obtained viral sequences were for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Extraction of the buffy coat and purification of RNA
Whole blood was collected, and total RNA was extracted from the buffy coat and purified The raw data were subjected to QC analyses using the FastQC tool (version 0.11.9) (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). mRNA-seq read quality control was done using Trimmomatic 27 (version 0.36) and STAR RNA-seq 28 (version STAR 2.5.4a) using 150 bp paired-end mode was used to align the reads (hg19).
HTSeq 29 (version 0.9.1) was to retrieve the raw counts and subsequently, R (https://www.R-project.org/), Bioconductor 30 and DESeq2 31 were used. Additionally, the RUVSeq 32 package was applied to remove confounding factors. The data were prefiltered keeping only genes with at least ten reads in total. The visualization was done using dplyr (https://CRAN.R-project.org/package=dplyr) and ggplot2 33 . Genes were categorized as significant differentially expressed with an adjusted p-value (pAdj) below 0.05 and a fold change > 2 for up-regulated genes and a fold change of < -2 for downregulated ones and then conducted gene enrichment analysis (https://www.gseamsigdb.org/gsea/msigdb).
for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Statistical analysis
Data were presented as the means ± s.e.m. (standard error of the mean) of all experiments with n = number of biological replicates. For comparison of RNA expression levels between two groups, data were presented as standard deviation in each group and were evaluated with a two-way ANOVA followed by Tukey's multiple comparisons test or a two-tailed unpaired t-test with Welch's correction using GraphPad PRISM (version 9.0).

This study was approved by the Institutional Review Board (IRB) of the Office of Research
Oversight / Regulatory Affairs, Medical University of Innsbruck, Austria (EC numbers: 1064/2021).

Data availability
The RNA-seq data of patients will be uploaded in GEO before publishing the manuscript.

Acknowledgments
Our gratitude goes to the patients who contributed to this study to advance our understanding of COVID-19. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.    This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 14, 2021.   This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 14, 2021. ; https://doi.org/10.1101/2021.05.11.21256862 doi: medRxiv preprint Figure 1 for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 14, 2021. ; https://doi.org/10.1101/2021.05.11.21256862 doi: medRxiv preprint Figure 2 for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 14, 2021. ; https://doi.org/10.1101/2021.05.11.21256862 doi: medRxiv preprint Figure 3 for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 14, 2021. ; https://doi.org/10.1101/2021.05.11.21256862 doi: medRxiv preprint Figure 4 for use under a CC0 license. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.