1. Introduction
The ongoing coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), presents a global threat to public health, and it is still taking its toll, with over 750 M cases worldwide, as well as over 6.8 M deaths (As of 12 April 2023). Along with the pandemic, variants emerged, including alpha, which first appeared in the United Kingdom (UK or B.1.1.7), beta, originally described in South Africa (SA or B.1.351), gamma, identified in Brazil (P.1), delta, which emerged in India (21A/S:478K or B.1.617.2), and omicron (B.1.1.529) (
https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html (accessed on 11 June 2023)). These variants show increased transmissibility (largely alpha, delta, and omicron), immune evasion (primarily beta and gamma, but also delta and omicron), or both (mostly delta and omicron) [
1,
2,
3,
4,
5]. The ongoing evolution of SARS-CoV-2 and the consequential emergence of various SARS-CoV-2 variants, as well as subvariants, such as “Omicron variants under monitoring”, as recently defined by the WHO (
https://www.who.int/activities/tracking-SARS-CoV-2-variants, updated on 18 January 2023), which vary in their sequence, transmissibility, the nature of the disease they cause, challenge the efficacy of currently available vaccines and therapeutic solutions.
VSV-ΔG-spike is a replication-competent vaccine candidate, based on the Vesicular Stomatitis Virus (VSV) platform, in which the VSV glycoprotein G (VSV-G) was replaced by a full-length human codon-optimized Spike glycoprotein (S) gene of SARS-CoV-2. [
6]. VSV-ΔG-spike propagation was performed by serial passaging in vero E6 and vero cells, a process that was accompanied by spontaneous accumulation of several mutations, until reaching a stable genetic version of the vaccine [
6,
7]. These abovementioned mutations include known spike mutations at N501, E484, Q493, and G685, sites of major importance to immune escape and virus transmissibility, as previously described [
8,
9].
Previously, we showed that sera derived from individuals participating in phase II clinical trial vaccinated with VSV-ΔG-spike efficiently neutralized the original virus, as well as alpha, beta, gamma, and delta [
7]. This observed preservation may be attributed to the unique sequence composed of spike mutations that correspond with SARS-CoV-2 variants’ mutations, which occurred naturally during growth and evolution of the VSV-ΔG-spike vaccine [
7].
In the current work, we aimed to explore the robustness of the VSV-ΔG-spike vaccine candidate against various SARS-CoV-2 variants. To that end, VSV-ΔG-spike efficacy was evaluated against authentic SARS-CoV-2 variants: alpha, beta, gamma, and delta, in comparison to the original virus, in K18-hACE2 mice. K18-hACE2s are transgenic mice designed to express the human ACE2 receptor under the keratin 18 promoter (K18) in epithelial cells in various tissues, including airway epithelia, namely, the route of SARS-CoV-2 infection, and they are widely used as SARS-CoV-2 infection models [
10]. It was previously reported that SARS-CoV-2 infection of K18-hACE2 mice targets the brain as a major site, leading to severe neurological disease, neuroinflammation, and neuronal death, contributing to infection-associated morbidity and mortality [
11,
12].
Here, we tested the efficiency of the VSV-ΔG-spike candidate vaccine against authentic SARS-CoV-2 viruses. We show that VSV-ΔG-spike vaccination induces neutralizing antibodies against all tested variants (original virus, alpha, beta, gamma, and delta), and it provides protection against lethal doses of each variant. Moreover, we show a significant prevention of viral load in organs of vaccinated mice, as well as a reduction in lung and brain tissue damage, compared to unvaccinated mice infected with each of the variants. RNAseq analysis of infected brains revealed that, whereas infection with any of the tested SARS-CoV-2 variants led to induction of cytokine storm and extensive neuroinflammation response, vaccination prior to infection completely prevented these deleterious phenotypes. We previously showed the efficacy of the VSV-ΔG-spike vaccine in serum samples derived from vaccinees participating in phase II clinical trial of VSV-ΔG-spike vaccine against variants alpha, beta, gamma, delta, and omicron [
7]. Thus, though not tested in this current work, we suggest that VSV-ΔG-spike may prove efficacious against the currently predominant variant omicron and its subvariants.
Taken together, we show that disease outcome and harmful phenotype caused by lethal dose of the original virus or variants alpha, beta, gamma, or delta can be prevented by VSV-ΔG-spike vaccination, irrespective of variant identity.
2. Methods
2.1. Cells
Vero E6 cells (ATCC CRL-1586™) were grown in DMEM containing 10% fetal bovine serum (FBS), MEM nonessential amino acids (NEAA), 2 mM L-glutamine, 100 Units/mL penicillin, 0.1 mg/mL streptomycin, and 12.5 Units/mL nystatin (P/S/N). Calu3 cells (ATCC HTB-55) were grown in RPMI, supplemented with 10% FBS, NEAA, 2 mM L-glutamine, P/S/N, and 1% Na-pyruvate. All reagents were from Biological Industries, Beit-Haemek, Israel. Cells were cultured at 37 °C, being in conditions of 5% CO2 with 95% humidity.
2.2. Viruses
Virus stocks of the SARS-CoV-2 original virus (GISAID accession EPI_ISL_406862) were propagated (four passages) in Vero E6 cells. SARS-CoV-2 variants were provided by the Central Virology Lab of the Israel Ministry of Health [
13]. Alpha (B.1.1.7, GISAID accession EPI_ISL_4169857) was passaged once in Vero E6, followed by two passages in Calu3 cells. Beta (B.1.351, GISAID accession EPI_ISL_4169885), gamma (P.1, GISAID accession EPI_ISL_4169886), and delta (B.1.617.2, GISAID accession EPI_ISL_4169986) variants were propagated in Vero E6 cells. Variants were verified by whole-genome sequencing (WGS) [
7]. All virus stocks were titered on Vero E6 cells, as previously described [
6]. Handling and working with the SARS-CoV-2 virus were conducted in a BSL3 facility, in accordance with the biosafety guidelines of the Israel Institute for Biological Research (IIBR).
2.3. VSV-∆G-Spike Preparation
VSV-∆G-spike was developed and generated, as previously described [
6], followed by additional passaging in Vero cells (WHO, Vero RCB 10–87), which were widely accepted by the WHO and regulatory agencies for the manufacturing of human viral vaccines, and they were further purified, formulated, and titered, as previously described [
6].
2.4. Animal Experiments
All animal experiments, involving SARS-CoV-2, were conducted in a BSL3 facility. Animal studies were approved by the Israel Institute for Biological Research Institutional Animal Care and Use Committee (IACUC). Female K18-hACE2 transgenic mice (The Jackson Laboratory, Bar Harbor, ME, USA), being 6–8 weeks old, were maintained at 20–22 °C and a relative humidity of 50 ± 10% on a 12 h light/dark cycle. Animals were fed with commercial rodent chow (Koffolk Inc., Neot Hovav, Israel) and provided with tap water ad libitum. Mice were randomly assigned to experimental groups.
K18-hACE2 mice were vaccinated intramuscularly (i.m.) with 107 pfu/mouse of VSV-ΔG-spike at a prime-boost (P + B) regimen during a three-week interval.
SARS-CoV-2 original virus, as well as alpha, beta, gamma, or delta variants, were diluted in phosphate-buffered saline (PBS), supplemented with 2% FBS (Biological Industries, Beit-Haemek, Israel). Approximately 4 weeks following boost vaccination, anesthetized animals (75 mg/kg ketamine, 7.5 mg/kg xylazine in PBS) were infected by (20 µL i.n. instillation) with a pre-determined lethal dose as follow: 40 LD50 of the original virus, or 1000 LD50 variants alpha, 30 LD50 of beta, 100 LD50 of gamma, or 100 LD50 of delta. Animals’ body weight was monitored daily. Animals were sacrificed at 2 and 5 dpi for the following analyses: (1) viral load in lungs (2 and 5 dpi), nasal turbinates (2 dpi), and brains (5 dpi); (2) histopathological analysis in lungs and brains (5 dpi); and (3) RNAseq analysis of brains (5 dpi). Viral load, histopathological analysis, and RNAseq analysis are described below. The numbers of animals per experiment are indicated for each experiment.
2.5. Plaque Reduction Neutralization Test
Vero E6 cells were seeded in 12-well plates (4 × 105 cells/well) and grown overnight in growth medium. Sera from mice following prime and boost vaccination were heat inactivated (HI) (56 °C for 30 min), and then they were diluted in twofold serial dilutions (between 1:40 and 1:5120) in 300 µL of infection medium (MEM containing 2% FBS with NEAA, glutamine, and P/S/N), mixed with 300 µL of 300 pfu/mL of original SARS-CoV-2, or each variant, and incubated at 37 °C, 5% CO2 for 1 h. Then, Vero E6 cells were infected with virus–serum mixtures (200 µL/well) and incubated at 37 °C and 5% CO2 for 1 h. 2 mL per well overlay (MEM containing 2% FBS and 0.4% tragacanth (Merck, Jerusalem, Israel)) were added to each well, and plates were incubated at 37 °C 5% CO2 for 72 h for original virus, beta, gamma, and delta variants or 96 h for alpha variant. Following incubation, the overlay was aspirated, and cells were fixed and stained with 1 mL/well of crystal violet solution. The number of plaques in each well was determined, and the serum dilution that neutralizes 50% of the virions (NT50) was calculated using Prism software (GraphPad Software Inc.).
2.6. Viral Load Determination in Organs
Viral load was determined for nasal turbinates (2 dpi), lungs (2 and 5 dpi), and brains (5 dpi). Organs were harvested, and they were processed, and infectious virus quantitation was performed by plaque assay, as previously described [
6]. Briefly, serial dilutions of extracted organs were prepared in infection medium (MEM containing 2% FCS) and used to infect Vero E6 monolayers in duplicates (200 µL/well). Plates were incubated for 1 h at 37 °C to allow viral adsorption. Then, 2 mL/well of overlay (MEM containing 2% FBS and 0.4% tragacanth (Merck, Jerusalem, Israel) was added to each well, and plates were incubated at 37 °C, 5% CO
2, for 72 h (for original, beta, gamma, and delta variants) or 96 h (for alpha variant). The medium was then aspirated, and the cells were fixed and stained with 1 mL/well of crystal violet solution (Biological Industries, Beit-Haemek, Israel). The number of plaques in each well was determined. Viral load, as well as LOD, were calculated based on volume of cell infection (200 µL/well), dilution factor, and tissue processing volume (1 mL for nasal turbinates, 1.5 mL for lungs and brains), and they are presented as pfu/organ.
2.7. Histopathology and SARS-CoV-2 Immunohistochemical Staining
Virus detection was performed by immunohistochemistry (IHC), and mice (n = 3 per group) were anesthetized, and then they were perfused transcardially with PBS. Lungs and brains were isolated and fixed in 4% paraformaldehyde at room temperature (RT) for 7 days, and then they were transferred to 70% ethanol for further fixation. Lungs were spread horizontally in the cassette, and brains were dissected into four coronal sections, showing the striatum, hippocampus, midbrain, cerebellum, and brainstem. Serial sections, 4 µm thick, were performed using the Leica BOND-MAX system (Leica Biosystems Newcastle Ltd., Newcastle upon Tyne, UK). For 3, 30-diaminobenzidine (DAB) and H&E staining, pictures were taken using the Olympus microscope (BX60, Serial No. 7D04032, Tokyo, Japan) equipped with the microscope’s camera (Olympus DP73, Serial No. OH05504) at objective magnifications of 1.25× and 10×.
Paraffin sections (4 microns thick) were cut, put on glass slides, and stained for viral detection by IHC. Immunohistochemical staining was performed on 4 µm paraffin sections using the Leica Bond max system (Leica Biosystems Newcastle Ltd., Newcastle upon Tyne, UK). Sections were dewaxed and pretreated with epitope-retrieval solution for 10 min (ER1, Leica Biosystems Newcastle Ltd., Newcastle upon Tyne, UK), followed by 30 min incubation with a rabbit anti SARS-CoV-2 antibody (in-house preparation of rabbit polyclonal anti-RBD, diluted 1:800). Pictures were taken using an Olympus microscope (BX60, serial NO. 7D04032), equipped with a microscope camera (Olympus DP73, serial NO. OH05504) at objective magnification of ×1.25, ×10.
2.8. Gene Expression and Functional Analysis
RNA was isolated from mice brains at 5 dpi using Qiagen RNeasy mini kits (Qiagen, Valencia, CA, USA) with an on-column DNase step (Qiagen, Valencia, CA, USA), according to the manufacturers’ instructions. RNA quantification was carried out in a Qubit fluorometer using the Qubit RNA HS assay kit (Invitrogen, Carlsbad, CA, USA). Quality control analysis of RNA integrity was performed using High Sensitivity RNA ScreenTape and the TapeStation Analysis software (Agilent Technologies, Santa Clara, CA, USA). All the samples reached a RIN (RNA Integrity Number) score of >7. RNA-seq was performed at the JP Sulzberger Columbia Genome Center (New York, NY, USA). Libraries were generated using the Illumina TruSeq stranded mRNA kit, according to the manufacturer’s instructions. Polyadenylated RNA enrichment was performed. Sequencing of the 100-bp paired-end reads was performed on the Illumina NovaSeq 6000 system. DESeq2 files were constructed by a comparison of the data of the virally infected mice with the corresponding carrier buffer control of the same time point for log2 fold change and p-values. The log2 fold change was converted to fold change values in some cases for simplicity.
Functional analysis was performed by Ingenuity Pathway Analysis (IPA®; Qiagen Redwood City, CA, USA) to evaluate the influence of the viral infection on biological canonical pathways and networks. IPA cutoff criteria for input list of differentially expressed genes was set to log2 fold-change > ±1.3 and adjusted p-value < 0.05. Graphical representation of networks with their scores were computed by IPA. Within networks, nodes represent genes and lines indicate biological relationships (direct or indirect) with other genes, based on the published literature within the IPA software. The network score is based on hypergeometric distribution and is calculated with the right-tailed Fisher’s exact test.
2.9. Statistical Analysis
Data were analyzed with GraphPad Prism 6 software. Exact p values are provided for each analysis. Statistical significance for neutralization of SARS-CoV-2 variants was determined by one-way ANOVA, as well as the Kruskal-Wallis test with Dunns’ multiple comparisons test. Statistical analysis for mortality was performed using the Log-Rank (Mantel-Cox) test. Significance for viral load analysis was determined by the Mann-Whitney nonparametric test, per organ (vaccinated vs. unvaccinated). For IPA, the network score was based on hypergeometric distribution and was calculated with the right-tailed Fisher’s exact test.
4. Discussion
Since the beginning of SARS-CoV-2 pandemic, emergence of new variants has been challenging the effectiveness of currently available vaccines, due to variants’ increased transmissibility’ and immune evasion [
3,
4,
14]. Here, we explore the efficacy of VSV-ΔG-spike vaccine against the original SARS-CoV-2, as well as against several variants: alpha, beta, gamma, and delta in the K18-hACE2 mice model. We show complete protection of K18-hACE2 mice infected with SARS-CoV-2 variants, decrease in viral load and prevention of tissue damage by VSV-ΔG-spike vaccination. We also show that vaccination prevents brain-related cytokine storm exemplified by the robust transcriptomic disease profile, irrespective of variant identity.
Preclinical studies with our VSV-ΔG-spike vaccine candidate in four animal species (mice, rabbits, hamsters, and pigs) demonstrated its excellent safety profile [
15,
16]. Soon to be published data from the clinical trial in Israel, aiming to evaluate the safety, immunogenicity, and potential efficacy of VSV-ΔG-spike vaccine candidate (BriLife
®) further indicates no significant adverse effects in humans.
Many have explored the ability of sera obtained from convalescent patients or human vaccinees to neutralize SARS-CoV-2 variants. Significant reduction in neutralization of the beta variant is often reported [
3,
4,
8,
9,
17], whereas only a minor reduction in neutralization of alpha variant is described [
2,
4]. Neutralization of the delta variant is reduced by ~2–11 fold, which is more significant than the alpha, but less than the beta variant neutralization [
4,
13,
18,
19,
20]. Here, we show significant neutralization of alpha, beta, gamma, and delta variants by vaccinated K18-hAEC2 mice sera with no significant difference in neutralization of the original virus. These data correspond with our previously reported results of sera samples from individuals receiving two-doses of VSV-∆G-spike as part of phase II clinical trial showing maintenance of the neutralization potential against several SARS-CoV-2 variants: alpha, beta, gamma, delta, and omicron, which were not tested in the current work [
7]. During its development, rVSV-ΔG-spike vaccine spontaneously acquired mutations at specific sites, including, but not exclusive to, the receptor-binding domain (RBD). Some of them are identical to or correspond to key mutations or areas that can be found in several SARS-CoV-2 variants: N501, E484, Q493, and G685. These mutations are known to affect SARS-CoV-2 variants’ transmissibility and immune evasion. We suggest that these acquired mutations enable vaccine efficacy against the variants tested here, and potentially towards current and future SARS-CoV-2 variants.
With omicron and its subvariants being the predominant variants, we also aimed at establishing a robust infection model of K18-hACE2 with omicron variant. Infection of K18-hACE2 mice with omicron did not result in significant weight loss, nor led to lethality, even at high infection doses. This is in accordance with others reopting an attenuated disease in rodents, including K18-hACE2 mice, showing limited weight loss, if any, and a low viral burden in lower and upper respiratory tracts in K18-hACE2 mice [
21], thus making the K18-hACE2 mouse model of low value in evaluating vaccine efficacy against Omicron. Hence, our work focuses on variants alpha, beta, gamma, and delta, presenting a lethal disease following infection.
Though the respiratory system is considered the major target of SARS-CoV-2, COVID-19 is a multi-organ disease, accompanied by central nervous system symptoms, such as memory and attention deficits, and neurological manifestations associated with severe prognosis [
22]. ACE2 is expressed in human brains, mostly in neurons, but also in non-neuronal cells, such as astrocytes, oligodendrocytes, and endothelial cells, thus making the brain susceptible to SARS-CoV-2 infection [
23,
24]. Additionally, autopsies of post-mortem patients’ brain tissues revealed profound neuroinflammation phenotype, with substantial immune activation in the CNS, neuropathology, as evident by astrocytosis, axonal damage, blood–brain barrier leakage, and viral antigen presence in ACE2-receptor-positive cells enriched in the vascular compartment [
25,
26]. Recapitulating COVID-19 brain pathology in animal models is of major importance for assessing SARS-CoV-2 derived damage, as well as vaccine efficacy.
K18-hACE2 mice develop a robust respiratory disease resembling COVID-19 upon SARS-CoV-2 infection, with an emphasis on the brain as a major SARS-CoV-2 target. A recent examination of the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with several SARS-CoV-2 variants showed that the original virus, as well as alpha, beta, and delta infection were restricted to neurons and appeared to spread from the olfactory bulb mainly in basally oriented regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage, or demyelination [
27]. The pathological damage to the lungs and brains of infected mice were recently found to be similar to the clinical symptoms of patients with severe COVID-19 disease, further suggesting that the K18-hACE2 model can recapitulate SARS-CoV-2 infection in humans [
10]. In this work, we show that K18-hACE2 succumbed to infection with lethal doses of the original virus, as well as alpha, beta, gamma, and delta. Additionally, the SARS-CoV-2 tissue distribution depicted herein, for unvaccinated animals, is in accordance with other reports on K18-hACE2 mice model. It was previously reported that K18-hACE2 mice infected with 10
5 pfu SARS-CoV-2 display the highest virus titers in lungs at three dpi, declining on days five and six dpi [
11], and that disease severity, as well as death in K18-hACE2 mice, are correlated with peak SARS-CoV-2 viral load in the brain at five to six dpi [
11]. Here, we show that VSV-ΔG-spike vaccination prevents all deleterious changes and effects, regardless of variant identity, in a similar manner. It was observed that there occurred prevention of weight loss and mortality, as well as elimination of virus presence, in lungs at two and five dpi and in brains at five dpi.
Viral load analysis showed high viral titers in unvaccinated mice brains, as well as histopathological evaluation, showing positive staining for SARS-CoV-2 in several brain areas in all brains. The robust VSV-ΔG-spike vaccination efficacy is further highlighted in these two analyses, where viral load is significantly reduced (below LOD), and viral antigen is not detected, irrespective of SARS-CoV-2 variant. In line with these findings, we also show a significant disease transcriptomic profile in the unvaccinated brains, representing disease state, as opposed to the vaccinated mice brains. Whether the lack of brain viral load, lack of viral antigen, and significant transcriptomic changes in vaccinated brain are derived from prevention of virus entry to the brain, or rather by its elimination, will require additional research, beyond the scope of this work. Nevertheless, thorough transcriptomic analysis of unvaccinated mice brains shed light on the extent and nature of brain disease caused by the different SARS-CoV-2 variants that are prevented by vaccination.
It was previously shown that early intranasal host response in the nasal epithelia, namely, site of infection, may underlie and precede disease severity [
28]. This is mediated by interferon-stimulated genes (ISG). In mild and moderate COVID-19, there is an expression of anti-viral/interferon-responsive genes, which mediate the protective response, whereas severe disease is accompanied by impaired anti-viral gene expression on one hand, and recruitment of highly inflammatory macrophages, as the source of proinflammatory cytokines, on the other hand [
28]. The significant reduction in virus titers we show in the nasal turbinates at two dpi in vaccinated mice may induce/correlate with ISG expression, leading to increased viral control, thus facilitating the lack of virus presence in the brain, altogether contributing to reduced disease severity and increased protection of K18-hACE2 mice. It is noteworthy that, whereas vaccination nearly eliminated virus presence in lungs and brains, in the nasal turbinates, at two dpi, there is a significant reduction in virus presence, but not complete exclusion. Hence, it would be interesting to explore intranasal administration of VSV-ΔG-spike, which may abolish viral presence, already, at the site of infection and reduce the risk of virus spread from infected individuals, thus providing an additional level of protection against SARS-CoV-2.
Due to the nature of the COVID-19 disease, and its impact on the respiratory system, many have explored gene expression patterns in lungs following SARS-CoV-2 infection. Winkler et al. [
29] thoroughly analyzed lung inflammation and damage following SARS-CoV-2 infection with the original virus at several time points, showing high levels of viral RNA and viral titers at several time points following infection, decline in pulmonary function, as well as extensive upregulation of the innate immune response, including hundreds of DEGs, many of which associated with IFN signaling, nuclear factor-κB-dependent cytokine responses, or leukocyte activation [
29]. Other work focused on RNAseq and bioinformatic data comparing transcriptional responses of mouse lungs (including K18-hACE2 mice), human lungs/lung tissues following SARS-CoV-2 infection showed differences between human and mice in single copy orthologues DEGs (scoDEGs), but high agreement in pathway analyses between species, specifically in immune signatures and inflammation pathways [
30].
Here, we also aimed to characterize the transcriptomic landscape of brains infected with SARS-CoV-2 variants, representing diseased brains. We compared the diseased brain following infection with the original virus to infection with alpha, beta, gamma, or delta variants, and moreover, we tested the effect of VSV-∆G-spike vaccination prior to infection on the brain transcriptome. We show that diseased brains at five dpi upregulate hundreds of genes, whereas there are hardly any DEGs in vaccinated brains. Gamma and delta infection led to the highest number of upregulated DEGs in unvaccinated brains. There is a high proportion of upregulated DEGs that are shared by all variants. The top genes upregulated by the original SARS-CoV-2 infection are also among the 10–20 top genes upregulated by the other variants tested in this work: alpha, beta, gamma, and delta. The same applies for the top canonical pathways expressed by all variants in set 1, in which infection with any of the variants leads to significant activation of the same specific pathways: the first two pathways are shared by all tested variants: “Role of Hypercytokinemia/hyperchemokinemia in the Pathogenesis of Influenza” and “Role of Pattern Recognition Receptors in Recognition of Bacteria and Viruses”, followed by “Agranulocyte adhesion and diapedesis” and “Granulocyte adhesion and diapedesis”, pathways concerning with migration of leukocytes and white blood cells to areas of infection, and inflammatory response to pathogens, which are also highly ranked. Hypercytokinemia, or cytokine storm, and hyperchemokinemia accompanied by an uncontrolled proinflammatory response, play important roles in COVID-19 pathogenesis, and induction of pro-inflammatory cytokines was also correlated with disease severity in animal models and humans [
12,
31]. This is represented by the top expressed and activated pathway in our analysis. Though this pathway’s title highlights “…Pathogenesis of Influenza”, the molecules involved in this pathway are very much relevant for SARS-CoV-2, as well, including many cytokines and innate immune molecules, namely, IFNβ, IL1, IL6, IL10, Cxcl10, Cxcl2, IRF3, IRF7, IRF9, etc. SARS-CoV-2 infection is a potent inducer of proinflammatory cytokines, and as such, strong upregulation and activation of this pathway upon SARS-CoV-2 infection by any of the tested variants is expected [
12,
31].
The second highest-enriched canonical pathway by all variants in “unvaccinated set 1” is “Role of Pattern Recognition Receptors in Recognition of Bacteria and Viruses” (
Figure S2). Pattern Recognition Receptors (PRRs) are important components of innate immunity, they recognize pathogen-associated molecular patterns (PAMPs), and this recognition leads to the expression of cytokines, chemokines, and co-stimulatory molecules, with the aim of eliminating infection and activating adaptive immunity [
32]. SARS-CoV-2 was previously shown to be recognized by several distinct classes of PRRs: Toll-Like Receptors (TLRs), Retinoic-acid Inducible Gene-I Receptors (RLRs), and C-type Lectin Receptors (CLRs) [
33]. Signaling through these PRRs leads to production of inflammatory cytokines and chemokines, as is also evident by the high activation of hypercytokinemia/hyperchemokiemia, as well as interferon signaling pathways.
In contrast, vaccination prior to infection with any of the tested variants leads to a relatively small amount of DEGs, mostly downregulated. For each variant in the vaccinated set, most DEGs (either upregulated or downregulated) were exclusive, with several genes overlapping between a few of the variants, and only two genes are shared by all variants. Such minimal overlap of DEGs between all vaccinated and infected mice with any of SARS-CoV-2 variants may imply that the vaccine prevents disease state by initiating a robust host response to the infection, such that the shared pathways and genetic signatures that were induced in the diseased animals are no longer activated. Additionally, canonical pathway expression and activation are less significant upon vaccination, as expressed by the relatively low p-value and z-scores. Thus, transcriptomic analysis highlights the ability of VSV-∆G-spike vaccination to prevent brain disease response caused by each of the tested variants in a similar manner.
Phagocytosis is an essential innate immune defense mechanism for the elimination of microbial pathogens. During viral infection, one important immune response is autophagy, which destroys the virus-infected cells, and viruses are often found to modulate it [
34]. Here, “phagosome formation” is the second most activated canonical pathway shared by all tested variants (
Figure S3). Formation is the first step of the three characteristic steps of phagocytosis: formation, maturation, and resolution [
35]. The strong activation of this pathway indicates an immediate defense response against viral infection that activates immune and inflammatory response against each one of the tested variants. Such activation, as part of disease state in the brain, has deleterious consequences.
An additional interesting and highly activated pathway in brains of all unvaccinated and infected mice is the “Neuroinflammation signaling pathway” (
Figure S4). Neuroinflammatory signaling plays a key role in maintaining CNS homeostasis, including destruction and removal of damaging agents and clearing injured neuronal tissue. However, uncontrolled inflammatory response can result in destruction of normal brain tissue.
When comparing the top 50 canonical pathways shared by all variants in the disease set to the top 50 canonical pathways common to all variants in the vaccinated and infected set, we show a pattern in which, in the “disease state”, all but one of the top 30 pathways are highly upregulated, and the trend is shared by all variants. “PPAR signaling” is the first pathway to be inactivated by all variants. The Peroxisome Proliferator-Activated Receptor (PPAR) family acts as a ligand-activated transcriptional regulator. PPARs are known to have an established immune modulatory role, thus their down regulation in disease state may abolish this activity, leading to exacerbated inflammatory damage. PPAR agonists were previously reported as effective adjuvants for COVID-19 vaccines by modifying immunogenetics [
36].
On the contrary, analysis of vaccinated and infected canonical pathways yielded mostly downregulated pathways among the top 50, if any. Original virus and delta-infected mice were the most affected groups in “vaccinated and infected set 2”, with the largest number of canonical pathways. For all variants following vaccination, −log(p-value) and z-score values were significantly lower, with a maximum of 6.2 −log(p-value) and z-score ranging between −3 to 1, representing mostly inactivation or no effect on canonical pathways. Taken together, canonical pathway analysis shows that vaccinated mice brains, following SARS-CoV-2 vaccination, are very much similar to naïve mice brains, indicating prevention of disease as a result of vaccination.
Interestingly, despite the massive gene expression and activation of many canonical pathways, including pathways characterized by an excessive immune response in unvaccinated infected mice, these animals ultimately succumbed to the disease caused by any of the tested variants. Moreover, these deleterious disease manifestations were not at all evident at brains of vaccinated mice, further strengthening vaccine efficacy in disease prevention. Whether the immune response prevented virus entry to the brain awaits further analyses.
As previously reported, VSV-∆G-spike vaccine development was accompanied by spontaneous acquisition of mutations that are similar or identical to some of the key naturally occurring mutations in SARS-CoV-2 variants [
7], which may be beneficial in VSV-∆G-spike efficacy against a variety of SARS-CoV-2 variants, and this may account for its robustness. The striking similarity in the disease transcriptomic landscape between the original virus in relation to alpha, beta, gamma, and delta variants, despite their differences in virulence and infectivity, may also account for variants that may emerge in the future. In this work, VSV-ΔG-spikes’ efficacy was evaluated against five different viruses and was able to prevent the deleterious infection outcomes, irrespective of variant identity, thus showing promising potential for VSV-ΔG-spike vaccine ability to also protect against omicron and its subvariants, as well as against future variants.
In light of the predominance of mRNA vaccines for SARS-CoV-2, we wish to point out several advantages of the VSV vaccination platform utilized by VSV-ΔG-spike. First, VSV-ΔG-spike has two mechanisms of action: the antigen (SARS-CoV-2 Spike) is presented on the surface of the virus particle, similarly to its presentation on the surface of SARS-CoV-2, and it is also expressed in vaccine-infected cells that further secrete it and/or present it to the immune system. In contrast, mRNA vaccines use only the latter. Second, the VSV-ΔG-spike may be applicable at various vaccination routes. The VSV-ΔG-spike vaccine is a live virus and the spike antigens presented on its surface mediate infection of target cells expressing hACE2. This enables utilization of this vaccine by variety of vaccination routes, such as i.m., i.n., or orally (p.o), to induce preferred local and systemic immune response.
Taken together, we suggest that VSV-∆G-spike provides robust protection against a variety of SARS-CoV-2 variants in a similar manner. Maintenance of neutralization potential against each variant, complete protection from morbidity and mortality, significant viral load reduction for some organs, and complete elimination for others, as well as “naïve-like” transcriptomic profile, regardless of variant identity, highlight the robustness of the VSV-∆G-vaccine, as well as its potential against future SARS-CoV-2 variants.