key: cord-0428882-ht7ofusm authors: Aid, Malika; Vidal, Samuel J.; Piedra-Mora, Cesar; Ducat, Sarah; Chan, Chi N.; Bondoc, Stephen; Colarusso, Alessandro; Starke, Carly E.; Nekorchuk, Michael; Busman-Sahay, Kathleen; Estes, Jacob D.; Martinot, Amanda J.; Barouch, Dan H. title: Ad26.COV2.S Prevents SARS-CoV-2 Induced Pathways of Inflammation and Thrombosis in Hamsters date: 2021-09-30 journal: bioRxiv DOI: 10.1101/2021.09.30.462514 sha: a27077cff738678aae61b95cc40d613223892671 doc_id: 428882 cord_uid: ht7ofusm Syrian golden hamsters exhibit features of severe disease after SARS-CoV-2 challenge and are therefore useful models of COVID-19 pathogenesis and prevention with vaccines. Recent studies have shown that SARS-CoV-2 infection stimulates type I interferon, myeloid, and inflammatory signatures similar to human disease, and that weight loss can be prevented with vaccines. However, the impact of vaccination on transcriptional programs associated with COVID-19 pathogenesis and protective adaptive immune responses is unknown. Here we show that SARS-CoV-2 challenge in hamsters stimulates antiviral, myeloid, and inflammatory programs as well as signatures of complement and thrombosis associated with human COVID-19. Notably, single dose immunization with Ad26.COV2.S, an adenovirus serotype 26 vector (Ad26)-based vaccine expressing a stabilized SARS-CoV-2 spike protein, prevents the upregulation of these pathways such that the gene expression profiles of vaccinated hamsters are comparable to uninfected animals. Finally, we show that Ad26.COV2.S vaccination induces T and B cell signatures that correlate with binding and neutralizing antibody responses. These data provide further insights into the mechanisms of Ad26.COV2.S based protection against severe COVID-19 in hamsters. Author Summary In this study, we show that vaccination with Ad26.COV2.S protected SARS-CoV-2 challenged hamsters from developing severe COVID-19 disease by attenuating excessive proinflammatory responses, such as IL-6 and IL-1, macrophages and neutrophils signaling. Ad26 vaccination also prevented the upregulation of pathways associated with thrombosis such coagulation and clotting cascades associated with infection, and the transcriptomic profiles of vaccinated animals were largely comparable to control uninfected hamsters. In contrast, SARS-CoV-2 challenged unvaccinated hamsters showed significant increase of these proinflammatory and prothrombotic pathways and significant weight loss compared to vaccinated hamsters. The COVID-19 pandemic has sparked intense interest in the rapid development of vaccines as well as animal models to evaluate vaccine candidates and to define molecular and immunologic correlates of protection. We and others have reported that animal models such as rhesus macaques and hamsters can be infected with SARS-CoV-2 and show robust viral replication in the upper and lower respiratory tract, enabling studies of COVID-19 pathogenesis and prevention with vaccines (1) (2) (3) (4) (5) . Golden Syrian hamsters show productive viral replication, lung pathology, and mortality when challenged with SARS-CoV2 (3, (6) (7) (8) , making them pertinent for vaccine evaluation. Reports of transcriptomics and proteomics profiling of blood and lung tissues from hamsters infected with SARS-CoV-2 have shown significant upregulation of interferon and proinflammatory pathways, activation of the complement system, and recruitment of neutrophils and macrophages to the lung of infected hamsters that correlates with the presence of SARS-CoV-2 viral RNA (6, 9) , supporting the role of these pro-inflammatory responses in COVID-19 severity (10) . Therefore, it is important to test whether vaccines developed against COVID-19 modulate the host immune and transcriptional responses and protect from excessive proinflammatory responses induced by SARS-CoV-2. We recently demonstrated that a single immunization with Ad26.COV2.S, an adenovirus serotype 26 (Ad26) vector-based vaccine expressing a stabilized SARS-CoV-2 spike protein, elicited binding and neutralizing antibody (NAb) responses and protected hamsters against SARS-CoV-2-induced weight loss, pneumonia, and mortality (3) . These preclinical data stimulated clinical trials that demonstrated immunogenicity and efficacy of Ad26.COV2.S in humans (11, 12) . In this study we performed in depth analyses of bulk RNA-Seq transcriptomic profiling of lung tissues at day 4 post SARS-CoV-2 challenge from Ad26.COV2.S vaccinated and sham unvaccinated hamsters. To characterize Ad26 vaccine-mediated protection from severe COVID-19 in hamsters, we integrated the transcriptomics data with virological data as well as adaptive immune responses elicited by Ad26 at weeks 2 and 4. We show that a single immunization with Ad26.COV2.S (also termed Ad26.S.PP) attenuated the upregulation of proinflammatory pathways and prevented the upregulation of thrombosis associated pathways such as platelet aggregation, blood coagulation and the clotting cascade. We also find that Ad26.COV2.S vaccination upregulated signatures of CD4+, CD8+, and B cell responses that correlated with the magnitude of Ad26-elicited humoral immune responses weeks following immunization. Together, these results provide new insights into the molecular and the immunological mechanisms of Ad26.COV2.S protection from severe COVID-19. We recently reported a study in which recombinant, replication-incompetent Ad26 vectors expressing SARS-CoV-2 Spike constructs prevented clinical disease in Syrian golden hamsters after challenge (3) . We studied two SARS-CoV-2 Spike constructs administered at either 10 10 or 10 9 viral particles (vp) as well as a sham control (5 total groups, n=10 per group). One Spike construct encoded a deletion of the transmembrane region and cytoplasmic tail reflecting the soluble ectodomain with a foldon trimerization domain (S.dTM.PP), while the other encoded the full-length spike (S.PP; renamed Ad26.COV2.S for clinical development). Both constructs contained two stabilizing proline mutations at the furin cleavage site. All 50 hamsters were challenged at four weeks with 5.0 × 10 5 TCID50 of the USA-WA1/2020 strain. While both constructs prevented severe weight loss, we observed that the S.PP construct exhibited superior immunogenicity and protection (3) . To gain insights into the molecular mechanisms of severe COVID-19 and Ad26-mediated protection in Syrian golden hamsters, we euthanized a subset of these animals at 4 days post-infection (4 dpi) with SARS-CoV-2 and performed bulk RNA sequencing (RNA-seq) from lung tissues. Specifically, we studied three sham unvaccinated hamsters, four S.dTM.PP-vaccinated hamsters (two at 10 10 and two at 10 9 ), and five S.PP-vaccinated hamsters (two at 10 10 and three at 10 9 ). For comparison, we additionally obtained lung tissues from three naïve animals that were neither vaccinated with Ad26 nor challenged with SARS-CoV-2 (S1A Fig). Transcriptomic profiling using bulk RNA-Seq analysis showed that the transcriptomic profiles of animals within each group of hamsters were homogenous, with no significant differences between groups of animals that received different vaccine doses (S1B and S1C Fig). We mapped the RNA-Seq reads to the SARS-CoV-2 genome and observed a significant number of reads that mapped to SARS-CoV-2 transcripts, ranging from 300 to 315000 reads, in the sham unvaccinated hamsters. However, we observed lower number of reads were mapped to SARS-CoV-2 transcripts in the S.dTM.PP-vaccinated hamsters, ranging from 10 to 16000 reads ( Fig 1A and S1D Fig). In contrast, SARS-CoV-2 reads were undetectable in the Ad26.S.PP vaccinated hamsters with only 0-12 reads mapping SARS-CoV-2 transcripts (Fig 1A and S1D Fig) . Consistent with these observations, we previously showed that detection of SARS-CoV-2 Envelope gene subgenomic RNA (sgRNA) by quantitative polymerase chain reaction (qPCR) was markedly diminished by Ad26 vaccination (3), validating our RNA-seq data. Differential expression genes analysis (DEGs) showed differences in genes upregulated or downregulated in vaccinated and sham unvaccinated animals compared to naïve animals (S2A Fig) An excessive inflammatory response to SARS-CoV-2 is a major cause of disease severity and death in COVID-19. Pathways of type I and II interferon responses, and proinflammatory cytokines and chemokines were previously reported in blood and lung tissue of hamsters infected with SARS-CoV-2 and in COVID-19 patients (6, 13, 14) (S3A and S3B Fig) . We interrogated our transcriptomic data in vaccinated hamsters compared to naïve and sham groups for proinflammatory pathways that have been shown to be central for the pathogenesis of severe COVID-19 (10, (15) (16) (17) . Gene set enrichment analysis (GSEA) of DEGs at 4 dpi showed pathways of interferon signaling, inflammasome, and proinflammatory cytokines signaling such as interferon alpha, TNF, IL-1, and IL-6 signaling pathways, were significantly increased in (Fig 1B and 1C) . Further, the expression of major proinflammatory cytokines and chemokines that contribute to SARS-CoV-2 pathogenesis, such as Il-6, Il-1a, Il-1b, Stat1/2/3, Furin, IL-2rg, Ccl21, Cxcl10, Csf, Tnf, Ccr6 and Ccr7, were significantly increased in sham unvaccinated (shown in green) compared to naïve hamsters (Wald test, BH-adjust P<0.05) (Fig 1D and 1E) . Whereas, the expression levels of these proinflammatory markers is comparable between vaccinated and naïve animals (Fig 1C and 1E) . Reduction of interferon signaling in vaccinated hamsters was validated by immunohistochemistry (IHC) for the interferon-inducible protein Mx1 in the lung of Ad26.S.PP vaccinated hamsters compared to sham infected unvaccinated hamsters (p=0.02) (Fig 1F) . Furthermore, we observed that proinflammatory pathways of IL6 JAK-STAT3 and TNF alpha signaling correlated negatively with Ad26 elicited ELISA binding and NAbs at weeks 2-4 following immunization (Fig 1G and S4 Fig). In contrast, these pathways correlated positively with viral loads in the lung and nares of vaccinated hamsters (Fig 1G, S4 Fig) . Previous bulk and single cell transcriptomics and proteomics studies in COVID-19 patients, macaques and hamsters infected with SARS-CoV-2 showed upregulation of macrophage and neutrophils signature that correlate with disease severity (4, 6, 10, 14, (16) (17) (18) . We observed that pathways related to monocyte (NES=1.79, FDR=0.006), macrophage M1(NES=2.41, FDR<10 -6 ) and M2 (NES=1.86, FDR=0.008), and neutrophil (NES=1.58, FDR=0.04) signaling were significantly increased in sham unvaccinated hamsters compared to naïve animals (Fig 2A) . Vaccination with Ad26 prevented the upregulation of these pathways in the Ad26.S.PP group at 4 dpi (Fig 2A and 2B ), whereas these pathways were slighter upregulated in the Ad26.S.dTM.PP vaccinated group (Fig 2A and 2B) . pathways (Fig 2A) . Further, we observed that markers of proinflammatory macrophages including Ccl3, Cd163, Cd68, Csf1r and Marco (6), were significantly increased in sham unvaccinated hamsters compared to naïve animals, while in the Ad26.S.PP vaccinated hamsters the expression of these macrophages markers was comparable to naïve animals (Fig 2B) . Consistent with these observations, vaccinated hamsters had fewer macrophages (Iba-1+) and neutrophils (myeloperoxidase, MPO+) by immunohistochemistry in lung, 4 days following challenge by quantitative image analysis (Fig 2C and 2D) . We observed that signatures of monocytes, M0, M1, M2 macrophages and neutrophils correlated negatively with Ad26 induced ELISA and NAbs titers at weeks 2-4 but were positively correlated with viral loads in the lung and nares of vaccinated hamsters (Fig 2E and We next examined the effect of Ad26 vaccination on additional pathways known to play prominent roles in the pathogenesis of severe COVID-19 including activation of complement and coagulation cascades (14, 19, 20 (Fig 3A) . Markers of these SARS-CoV-2 induced pathways in sham hamsters were enriched in complement components C3, C7, C2; clotting and coagulation cascade and tissue factors F3, F5, Plau, Fga; and platelet activation and aggregation markers such as Clu, Timp1, Thbs1, Sh2b2, and Vav1 (Fig 3B) . These markers were less or not increased in the Ad26 vaccinated hamsters (Fig 3B) . When (Fig 3C) , suggesting Ad26 vaccination protected hamsters from developing thrombosis associated pathways. Moreover, we found that pathways of complement cascade, platelet activation and aggregation, coagulation cascade, and fibrinolysis correlated negatively with weeks 2-4 ELISA binding and neutralizing antibody titers (NAb) elicited by Ad26 vaccination, and positively correlated with 4 dpi viral loads in lung and nares of vaccinated hamsters (Fig 3D and S6 Fig) . We next investigated whether Ad26 vaccination induced molecular signatures consistent with protective adaptive immune mechanisms observed in SARS-CoV-2 infected humans, macaques and hamsters vaccinated with Ad26 (3, 11, 21, 22) . GSEA showed that signatures of (Fig 4A and 4B) . GO term analysis of CD4+ T cell pathway (Cd28, Cd3g, Cd3e, Cd3d, Lck, Fyn, Cd247, Ccr7, Dpp4, Icos), showed enrichment of markers of lymphocyte co-stimulation (FDR=2.65 x 10 -12 ), leukocyte activation (FDR=7.25 x10 -12 ), and T cell activation (FDR=1.18 x 10 -10 ). Moreover, key markers of cytotoxic CD8+ T cell adhesion, activation, and effector memory functions such as Prf1, Lef1, Lck, Tcf7 and Nkg7, were increased by Ad26 in vaccinated compared to naïve control animals (Fig 4B) , suggesting CD8+ T cell-induced by Ad26 are engaged in the recognition and elimination of infected cells. Similarly, Th1 cell signatures were significantly upregulated in Ad26-vaccinated compared to naïve control animals and include Th1 cytotoxic genes (NES=2.01, FDR<10 -6 ), genes increased in Th1 cells (NES=1.6, FDR=0.02) and STAT4 signaling, a major transcription factor involved in Th1 differentiation and proliferation (23) (NES=1.31, FDR=0.02) (Fig 4C) . We next correlated T cell signatures induced in vaccinated hamsters with humoral immune responses elicited by the Ad26 vaccine at weeks 2-4. We observed a significant positive correlation of markers of CD4+ T cell (Fig 4D) , signatures of CD8+ T cell and T regulatory cell signatures with ELISA and neutralizing antibody titers at weeks 2 and 4, whereas these signatures correlated negatively with viral loads in the lung and nares of vaccinated hamsters (Fig 4E) . Similarly, we found that STAT4 targets, a major transcription factor involved in Th1 activation (23, 24), were positively correlated with ELISA and neutralizing antibody titers at weeks 2 and 4, and negatively correlated with viral loads in the lung and nares of vaccinated hamsters (Fig 4E) . Given the critical role of humoral responses in immunity to SARS-CoV-2 infection (25), we next characterized B cell responses following Ad26 vaccination in hamsters. We observed a significant increase in signatures of B cell activation and differentiation and other markers regulating B cell fate and development in vaccinated hamsters compared to sham and naïve animals, as shown by the SLEA scores in each individual animal (Fig 5A) and by the upregulation of the top markers of B cell activation and differentiation and markers regulating B cell fate and development such as Cd79, Pou2af1, Bank1and Tcf4 (Fig 5B) . These pathways correlated positively with binding and neutralizing antibody titers elicited by Ad26 at weeks 2 and 4 (Fig 5C and S7 Fig) but were correlated negatively with viral loads in the lung and nares of vaccinated hamsters (Fig 5C and S7 Fig) . Additionally, pathways of B cell activation, differentiation, and development correlated positively with S-specific and RBD-specific antibody IgG, IgG2a, IgG3, IgM, Fc-receptors FcRγ2, FcRγ3, and FcRγ4 and antibody-dependent complement deposition (ADCD) responses in vaccinated animals at week 4 assessed by systems serology and previously published by our group (3) (Fig 5D) . We previously showed that vaccination with Ad26.COV2.S, an Ad26 vector encoding a full-length prefusion stabilized S immunogen (S.PP), protected against SARS-CoV-2 challenge in hamsters (3) and macaques (21) . Here we used bulk RNA-Seq of the lung, and showed that Ad26-based vaccines both attenuate hallmark pathologic pathways associated with COVID-19 and stimulate signatures of protective adaptive immune responses. These observations are consistent with a decrease of lung injury and tissue damage in vaccinated hamsters reported previously (3) (26) and highlight the protective role of Ad26 vaccination from SARS-CoV-2 pathology in both macaques and hamsters (3, 21) . In a Phase 1/2a randomized, clinical study, Ad26.COV2.S vaccination generated robust CD8+ T cell, CD4+ T cell, and antibody responses following vaccination (11, 12) . Consistent with these clinical data, in the Ad26 vaccinated hamsters we observed upregulation of signatures of effector and cytolytic CD8+ T cell, CD4+ T cell signatures. More importantly, we showed that these signatures correlated positively with binding and neutralizing antibody titers elicited by Ad26 at weeks 2 and 4 and negatively with viral loads in lungs and nares at 4 dpi. Alignment was performed using STAR version 2.7.3a (27) with the MesAur1.0 (GCF_000349665.1) assembly and annotation of the hamster downloaded from NCBI. Transcript abundance estimates was calculated internal to the STAR aligner using the algorithm of htseqcount as described previously (4) . DESeq2 was used for normalization, producing both a raw and normalized read count table. Differential expression at the gene level were performed by DESeq2 implemented in the DESeq2 R package. A corrected p-value cut-off of 0.05 was used to assess significant genes that were upregulated or down regulated by SARS-Cov2 at day 4 post challenge in sham and vaccinated animals compared to naïve controls and in vaccinated compared to sham animals using Benjamini-Hochberg (BH) method. Tissues were fixed in freshly prepared 4% paraformaldehyde for 24 hours, transferred to 70% ethanol, paraffin embedded within 7-10 days, and blocks sectioned at 5 µm. Slides were baked for 30-60 min at 65 o C then deparaffinized in xylene and rehydrated through a series of graded ethanol to distilled water. For Iba-1 IHC, heat induced epitope retrieval (HIER) was performed using a pressure cooker on steam setting for 25 minutes in citrate buffer (Thermo; AP-9003-500) followed by treatment with 3% hydrogen peroxide. Slides were then rinsed in distilled water and protein blocked (BioCare, BE965H) for 15 min followed by rinses in 1x phosphate buffered saline. Primary rabbit anti-Iba-1 antibody (Wako; 019-19741 at 1:500) was applied for 30 minutes followed by rabbit Mach-2 HRP-Polymer (BioCare; RHRP520L) for 30 minutes then counterstained with hematoxylin followed by bluing using 0.25% ammonia water. Labeling for Iba-1 was performed on a Biogenex i6000 Autostainer (v3.02). In some cases, Iba The canonical pathway and upstream regulator functions of IPA core expression analysis tool (Qiagen) were used to interrogate the lists of genes upregulated or down regulated by SARS-CoV-2 at 4 dpi in sham infected hamsters compared to macaques infected with SARS-CoV-2. Canonical pathways and upstream regulators were considered significant if pathway activation Z-Score ≥ 2 and pathway overlap corrected p-value < 0.05 (using the Benjamini-Hochberg method). Functional analysis of statistically significant gene and protein changes was performed using Ingenuity Pathways Knowledge Base (IPA; Ingenuity Systems). For all gene set enrichment analyses, a right-tailed Fisher's exact test was used to calculate P-values associated with each biological function and canonical pathway. The calculated z-score signifies whether a gene or protein expression changes for known targets of each regulator are consistent with what is expected from the literature (z > 2, regulator predicted to be activated, z < -2, regulator predicted to inhibited). Additional functional module analyses were performed using Functional module detection (https://hb.flatironinstitute.org/module/). Go term enrichment analysis was performed using the GeneMania database (https://genemania.org/). For GSEA analysis, all significant pathways and molecular signatures, up-or down-regulated in the different groups of hamsters, were selected using a false discovery rate < 20 and a nominal p-value < 0.05. In IPA Global Canonical Pathways (GCP), a multiple-testing corrected p-value was calculated using the Benjamini-Hochberg (BH) method. Quality of RNA-Seq raw reads was examined using FastQC Heatmap presenting the SLEA z-score of each of these pathways. SLEA represent the mean expression of all significant genes within each pathway for each animal. An SLEA z-score greater than 0 corresponds to a pathway for which member genes are up-regulated while an SLEA z-score inferior to 0 corresponds to a pathway with genes downregulated in that sample. Columns correspond to individual animals and rows correspond to individual pathways scaled across all animals using the z-score R function. Negative correlations were shown in cyan and positive correlations were shown in pink. Correlations were assessed using Spearman correlation. Scatter plot of the SLEA score of IL6_JAK_STAT3 and TNF alpha pathways in vaccinated animals at 4 dpi as a function of the ELISA binding titers, neutralizing antibody titers and viral load. The x-axis represents the samples level enrichment score of each pathway (SLEA) and the y-axis shows immune responses elicited by Ad26 in vaccinated hamsters. The average expression of the genes within each pathway was calculated using the SLEA z-score method. A linear regression model (blue or red line), and its 95% confidence interval (gray zone), was fit between SLEA z-score, viral load and the different antibody responses. A Spearman correlation and a t test were performed to assess the significance of the correlation between pathways SLEA scores and each response. Each red dot correspond to an individual animal. Positive correlations were shown in red and negative correlation were shown in blue. SARS-CoV-2 infection protects against rechallenge in rhesus macaques Respiratory disease in rhesus macaques inoculated with SARS-CoV-2 Ad26 vaccine protects against SARS-CoV-2 severe clinical disease in hamsters Vascular Disease and Thrombosis in SARS-CoV-2-Infected Rhesus Macaques Animal models for COVID-19 Longitudinal omics in Syrian hamsters integrated with human data unravel complexity of moderate immune responses to SARS-CoV-2. bioRxiv Pathogenic and transcriptomic differences of emerging SARS-CoV-2 variants in the Syrian golden hamster model Pathogenesis and transmission of SARS-CoV-2 in golden hamsters Syrian hamsters as a small animal model for SARS-CoV-2 infection and countermeasure development Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophages Interim Results of a Phase 1-2a Trial of Ad26.COV2.S Covid-19 Vaccine Safety and Efficacy of Single-Dose Ad26.COV2.S Vaccine against Covid-19 Dysregulation of type I interferon responses in COVID-19 Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19 COVID-19 cytokine storm: The anger of inflammation Impaired immune cell cytotoxicity in severe COVID-19 is IL-6 dependent Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19 Baricitinib treatment resolves lower airway inflammation and neutrophil recruitment in SARS-CoV-2-infected rhesus macaques Pulmonary Vascular Endothelialitis, Thrombosis, and Angiogenesis in Covid-19 Factors associated with COVID-19-related death using OpenSAFELY Single-shot Ad26 vaccine protects against SARS-CoV-2 in rhesus macaques Activated STAT4 has an essential role in Th1 differentiation and proliferation that is independent of its role in the maintenance of IL-12R beta 2 chain expression and signaling Signal transducer and activator of transcription 4 is required for the transcription factor T-bet to promote T helper 1 cell-fate determination Correlates of protection against SARS-CoV-2 in rhesus macaques. Nature. 2020. 26. van der Lubbe JEM Ad26.COV2.S protects Syrian hamsters against G614 spike variant SARS-CoV-2 and does not enhance respiratory disease We thank Kathryn L. Pellegrini, Gregory Tharp, Steven E. Bosinger, Hanneke Schuitemaker,