Single Dose of a VSV-Based Vaccine Rapidly Protects Macaques From Marburg Virus Disease

Marburg virus (MARV) is a member of the filovirus family that causes hemorrhagic disease with high case fatality rates. MARV is on the priority list of the World Health Organization for countermeasure development highlighting its potential impact on global public health. We developed a vesicular stomatitis virus (VSV)-based vaccine expressing the MARV glycoprotein (VSV-MARV) and previously demonstrated uniform protection of nonhuman primates (NHPs) with a single dose. Here, we investigated the fast-acting potential of this vaccine by challenging NHPs with MARV 14, 7 or 3 days after a single dose vaccination with VSV-MARV. We found that 100% of the animals survived when vaccinated 7 or 14 days and 75% of the animal survived when vaccinated 3 days prior to lethal MARV challenge. Transcriptional analysis of whole blood samples indicated activation of B cells and antiviral defense after VSV-MARV vaccination. In the day -14 and -7 groups, limited transcriptional changes after challenge were observed with the exception of day 9 post-challenge in the day -7 group where we detected gene expression profiles indicative of a recall response. In the day -3 group, transcriptional analysis of samples from surviving NHPs revealed strong innate immune activation. In contrast, the animal that succumbed to disease in this group lacked signatures of antiviral immunity. In summary, our data demonstrate that the VSV-MARV is a fast-acting vaccine suitable for the use in emergency situations like disease outbreaks in Africa.


INTRODUCTION
In 1967, a newly emerging pathogen caused hemorrhagic disease outbreaks with primary cases in Marburg, Germany (1) and Belgrade, Yugoslavia (2). Infected patients presented with fever, nausea, vomiting, conjunctivitis and later on petechiae as well as excessive bleeding suggesting clotting defects (3). In total, 32 people were infected, 7 of whom died (21.9% case fatality rate (CFR)) (1). Marburg virus (MARV) was identified as the causative agent and has since caused infrequent outbreaks in Africa with a CFR up to 90% (3). MARV is a member of the Filoviridae, has a 19 kb (-) single-stranded genome encoding 7 proteins, and forms filamentous particles that bud from the surface of infected cells (4). There is neither an approved treatment nor a vaccine against Marburg virus disease (MVD) and due to its high pathogenicity and effective human-to-human transmission MARV is classified as a select agent in the United States and the World Health Organization (WHO) added it to its list of priority pathogens (5).
Several vaccine approaches have been explored ranging from inactivated virus, subunit, and viral vector vaccines (6). One of the most promising vaccines is based on vesicular stomatitis virus (VSV) expressing the MARV glycoprotein (GP) as the viral antigen. A single dose of this vaccine is uniformly protective in nonhuman primates (NHPs) when administered 28 days prior to challenge (7,8). It also demonstrated a great potential for use in a post-exposure scenario in NHPs (9). This vaccine is very similar to the VSV-based Ebola virus (EBOV) vaccine VSV-EBOV (also known as rVSV-ZEBOV) which has been approved by the US Food and Drug Administration (FDA) (10) and the European Medicines Agency (EMA) under the name "Ervebo" for human use to prevent Ebola virus disease (EVD) (11). A single dose of VSV-EBOV has been shown to protect NHPs and humans within 7-10 days from lethal disease (12,13) and, more recently, it has been shown that a low dose of 10 PFU protects NHPs from disease (14). Because of certain similarities between MVD and EVD and their respective VSV-based vaccines, we wanted to explore the fast-acting potential of VSV-MARV in NHPs. We demonstrate that a single dose of VSV-MARV protects NHPs uniformly from lethal disease when administered up to 7 days before challenge. Notably, three out of four NHPs survived when the vaccine was administered 3 days prior to MARV infection with only one animal developing moderate disease. All control animals succumbed after developing classical signs of MVD as did one NHP in the day -3 vaccine group. In the day -3 group, whole blood (WB) transcriptional analysis revealed strong innate immune activation only in surviving NHPs. Collectively, these data demonstrate that the VSV-MARV is a fast-acting vaccine suitable for use in outbreak situations.

Ethics Statement
All the work involving infectious MARV was performed following standard operating procedures (SOPs) approved by the Rocky Mountain Laboratories (RML) Institutional Biosafety Committee (IBC) in the maximum containment laboratory at the RML, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health. Animal work was performed in strict accordance with the recommendations described in the Guide for the Care and Use of Laboratory Animals of the National Institute of Health, the Office of Animal Welfare and the United States Department of Agriculture and was approved by the RML Animal Care and Use Committee (ACUC). Procedures were conducted by trained personnel under the supervision of veterinary staff on animals anesthetized with ketamine. All efforts were made to ameliorate animal welfare and minimize animal suffering in accordance with the Weatherall report on the use of nonhuman primates in research (https://royalsociety.org/policy/publications/2006/ weatherall-report/). Animals were housed in adjoining individual primate cages that enabled social interactions, under controlled conditions of humidity, temperature, and light (12-h light:12-h dark cycles). Food and water were available ad libitum. Animals were monitored and fed commercial monkey chow, treats, and fruit at least twice a day by trained personnel. Environmental enrichment consisted of commercial toys, music, and video. Endpoint criteria based on clinical score parameters as specified and approved by the RML ACUC were used to determine when animals were humanely euthanized.

Animal Cohorts and Study Design
Sixteen male and female cynomolgus macaques (Macaca fascicularis) 5-7 years of age and 4.0-6.7 kg in weight were used for this study. Three groups of cynomolgus macaques (n=4 per group) were vaccinated with a single intramuscular (IM) injection of 1x10 7 plaque forming units (PFU) VSV-MARV at 14, 7 or 3 days before challenge. Control animals were IMvaccinated with 1x10 7 PFU VSV-EBOV (n=1 each day -14 and -7, n=2 day -3; n=4 total). Challenge was performed by IM injection of 1,000 PFU of MARV (confirmed by back-titration) as previously described (8). Clinical exams including a blood draw on anesthetized animals were conducted on days -14, -7, -3, 0, 3, 6, 9, 14, 21, 28, 35 and 42 post-challenge. The animals were observed at least twice daily for clinical signs of disease according to a RML ACUC-approved scoring sheet and humanely euthanized when they reached endpoint criteria. The study ended on day 42 post-challenge when all surviving animals were humanely euthanized.

Hematology and Serum Chemistry
The total white blood cell, neutrophil, lymphocyte, and platelet counts were determined from EDTA blood with the IDEXX ProCyte DX analyzer (IDEXX Laboratories, Westbrook, ME). Serum biochemistry including aspartate aminotransferase (AST), alkaline phosphatase (ALP), glucose, creatinine, and total bilirubin was analyzed using the Piccolo Xpress Chemistry Analyzer and Piccolo General Chemistry 13 Panel discs (Abaxis, Union City, CA).

Virus Loads
Viremia was determined from EDTA WB samples using Vero E6 cells (mycoplasma negative). To this end, cells were seeded in 48well plates the day before titration. At the day of titration, blood samples were thawed, and 10-fold serial dilutions were prepared in DMEM without supplements. Media was removed from cells and triplicate wells were inoculated with each dilution. After one hour, DMEM supplemented with 2% FBS, penicillin/streptomycin and L-glutamine was added and cells were incubated at 37°C. Cells were monitored for cytopathic effect (CPE) and 50% tissue culture infectious dose (TCID 50 ) was calculated for each sample employing the Reed and Muench method (18).

Humoral Immune Responses
Post-challenge NHP sera were inactivated by gamma irradiation (5 MRad) (19) and removed from the maximum containment laboratory according to RML standard operating procedures (SOP) approved by the RML IBC. Antibody titers were determined using ELISA kits based on recombinant soluble MARV-Angola GPdTM (Alpha Diagnostics, San Antonio, TX). Serum samples were diluted 1:200 for assessment in the ELISA. The ELISA was performed, and titers were calculated as per manufacturer's instructions. cDNA Library Preparation for RNA Sequencing RNA quality and concentration were assessed using an Agilent 2100 bioanalyzer. A total of 100 ng was used for library preparation. Total RNA samples were treated with 1 U of DNase I (New England BioLabs, Ipswich, MA) at 37°C for 10 min and cleaned with RNAClean XP beads (Beckman Coulter, Brea, CA). Ribosomal RNAs were depleted using rRNA removal beads (Illumina, San Diego, CA). The remaining RNA was fragmented at 94°C for 8 min to yield a median fragment size of 155 nucleotides (nt). Libraries were prepared using the TruSeq stranded total RNA library prep kit according to the instructions provided by the manufacturer (Illumina, San Diego, CA). Final library quality was assessed using an Agilent high-sensitivity DNA kit. Libraries with unique barcoded adaptors were pooled and sequenced in a single end read (1x100) on the Illumina NextSeq 500 platform.

Bioinformatics Analysis
Data analysis was performed using the RNA-Seq workflow module of the systemPipeR package (20). RNA-Seq reads were demultiplexed, quality filtered, and trimmed using Trimmomatic with an average phred score cutoff of 30 and minimum length of 50 bp. Quality reports were generated with the FastQC function. Trimmed reads were mapped to the Macaca fascicularis reference genome (Macaca_fascicularis.Macaca_fascicularis_ 5.0.dna.toplevel.fa) using HISAT2 and the corresponding gene annotation (Macaca_fascicularis.Macaca_fascicularis_5.0.94.gtf) from Ensembl. Uniquely mapped reads were counted using summarizeOverlaps in strand-specific mode. Normalization and statistical validation of differentially expressed genes (DEGs) was performed using EdgeR package's pair-wise function (21). 0 DPV data were used as the reference. DEGs were defined as protein coding genes with human homologues with at least a 2-fold change in expression, a multiple hypothesis Benjamini-Hochberg false discovery rate (FDR) corrected P value less than 0.05 and an average of at least 5 read per kilobase of transcript per million mapped reads (RPKM).
Temporal gene expression patterns and signatures that distinguish vaccine groups, survivors and non-surviving, and negative controls were analyzed using maSigPro, which is a twoway regression-based approach that finds a set of statistically significant DEGs for the entire time course (22). Only proteincoding genes with human homologs and an average of at least 5 read per kilobase of transcript per million mapped reads (RPKM) were included in this analysis.

Functional Enrichment and Data Visualization
DEGs were first mapped to human homologs using BioMart (Ensemble Gene 94). Only protein-coding genes with human homologs were included for further analysis. The functional enrichment of DEGs was assessed using Metascape (23).
Heatmaps, Venn diagrams, bar graph and volcano were generated using R packages VennDiagram, dplyr, and ggplot2. Line graphs were generated using GraphPad Prism V8 (San Diego, CA).

Statistical Analysis
Clinical data were examined for statistical significance using Prism version 9 (GraphPad, San Diego, CA). Survival curves were analyzed with Mantel-Cox test and values representing groups were analyzed by two-way ANOVA with Tukey's multiple comparisons. Statistical analysis of the maSigPro data was carried out using Prism version 8 (GraphPad, San Diego, CA). Significance was determined using a one-way ANOVA with a Dunnett's multiple-comparison test. Statistically significant differences are indicated as follows: p<0.0001 (****), p<0.001 (***), p<0.01 (**), and p<0.05 (*).

VSV-MARV Vaccination Protects NHPs Within 7 Days From Lethal Disease
The minimum time between vaccination and challenge for protection with a single IM dose of 1x 10 7 PFU VSV-MARV was determined by vaccinating groups of 4 NHPs at 14, 7 and 3 days prior to challenge. The control group consisted of 4 VSV-EBOV-vaccinated NHPs that were vaccinated with a single IM dose of 1x 10 7 PFU on day -14 (n=1), day -7 (n=1) and day -3 (n=2). On day 0, all NHPs were challenged with 1,000 PFU MARV by IM injection. The control animals developed signs of MVD and were humanely euthanized 6-and 8-days postchallenge (DPC) when they reached IACUC-approved endpoint criteria ( Figure 1A). One NHP in the day -3 vaccine group developed clinical signs of MVD as demonstrated by the increase in clinical score ( Figure 1B) and was euthanized 7 DPC ( Figure 1A). Of the remaining three animals in the day -3 group, one developed moderate signs of MVD, two developed very mild signs of MVD, and all three recovered ( Figure 1B). None of the NHPs in the day -14 or -7 vaccine groups developed signs of MVD ( Figure 1B). Only the 5 NHPs that succumbed to disease developed hallmarks of MVD including thrombocytopenia ( Figure 1C), high titer viremia ( Figure 1D), and increased levels of AST ( Figure 1E) and ALP ( Figure 1F). Other parameters examined in WB and serum after challenge demonstrated changes to abnormal levels for cell populations and metabolites for the 5 NHPs that succumbed to MVD ( Figure S1). Similarly, only these 5 NHPs developed cytokine levels indicative of the disease-associated cytokine storm ( Figure S2).

Surviving NHPs Developed MARV GP-Specific IgM and IgG Responses
VSV-based vaccines have been shown to primarily confer protection by humoral responses (15,24). Therefore, antigen- . Statistical significance among groups was determined using two-way ANOVA with Tukey's multiple comparisons. Statistical significance is indicated as **p < 0.01, and *p < 0.05. specific immune responses were determined in serum samples collected from all NHPs throughout the study. Levels of MARV GP-specific IgM peaked in the day -14 group 11 days postvaccination (DPV; -3 DPC) before dropping by DPC 0 (Figure 2A). MARV GP-specific IgM was detected in the day -7 group 3 DPC increasing slightly 6 DPC dropping 9 and 14 DPC (Figure 2A). In the day -3 group MARV-GP specific IgM was detected 6 DPC increasing slightly 9 and 14 DPC (Figure 2A). Levels of MARV GP-specific IgG were detected 11 DPV (day -14 group), 3 DPC (day -7 group) and 9 DPC (day -3 group) ( Figure 2B) coinciding with peak levels of MARV GPspecific IgM (Figure 2A). For the day -14 vaccination group the MARV challenge served as a boost as demonstrated by the increase in MARV GP-specific IgG 6 and 9 DPC ( Figure 2B).

VSV-MARV Vaccination Induces Transcriptional Changes Indicative of Antiviral Defense and B Cell Activation
Longitudinal transcriptional changes in WB in response to VSV-MARV and VSV-EBOV vaccination were assessed using RNA sequencing (RNA-Seq). To that end, transcriptomes of all available WB samples collected on 3, 7, 11, and 14 DPV were compared to those observed 0 DPV. Gene expression changes were evident as early as 3 DPV with VSV-MARV with 88, 73, 483 and 828 DEGs detected 3, 7, 11, and 14 DPV, respectively ( Figure 3A). The largest number of DEGs was detected 14 DPV ( Figure 3A). DEGs enriching to gene ontology (GP) terms associated with innate immunity to viral infection such as "defense response to virus" and "response to interferon alpha" were detected 3 DPV and remained differentially expressed throughout 14 DPV ( Figure 3B). Notable DEGs in this group include interferon stimulated genes (ISG; ISG15, BST2, IRF7, IFIH1) as well as genes that play a role in antigen presentation (HLA-B), cytokine production (MIF) and antiviral response (RSAD2, CARD9, DDX58, DDX60, and DHX58) ( Figure 3C). DEGs important for lymphocyte activation and adaptive immunity (e.g., GO terms "lymphocyte activation" and "lymphocyte mediated immunity") were detected 11 DPV and included genes important for T cell activation and effector functions (CD3D, CD8B, IL2RA, IL2RB1, GZMB and CD274). DEGs important for B cell activation and antibody secretion (MZB1, BATF, IGHM, IGHA1) were upregulated 14 DPV (Figures 3B, C).

Animals Vaccinated as Early as 7 Days Prior to Challenge Generate a Recall Adaptive Immune Response and Show Limited Gene Expression Changes
Next, we examined gene expression changes post MARV challenge in the day -14 and -7 groups. All comparisons were made relative to 0 DPV in a pairwise manner in each group. Limited transcriptional changes were detected with the exception of 9 DPC in the day -7 group ( Table 1). DEGs detected in the day -14 group at 3 and 6 DPC mapped to GO terms associated with viral life cycle, leukocyte migration, and lymphocyte activation ( Figure 4A). DEGs important for lymphocyte activation include genes involved in T cell activation (CD3E, decreased 9 DPC), and B cell activation and antibody secretion (IGLL1, IGHG2, IGHD, represents a Gene Ontology term; the size of each bubble represents number of DEGs enriching to that GO terms while color intensity represents the statistical significance (shown as -log 10 of the q-corrected p-value). (C, F) Heatmap of DEGs enriched to (C) "positive regulation of cytokine production", "defense response to virus" and "lymphocyte activation"; (F) "defense response to virus", "activation of innate immune response", "type I interferon signaling pathway", "lymphocyte proliferation", and "regulation of B cell proliferation." Each column represents the median normalized transcript counts (RPKM) for each gene at each time point. The range of colors is based on scaled and centered RPKM values of the entire set of genes, with red indicating highly expressed genes and blue indicating lowly expressed genes. and IGHA1; increased 9 DPC) ( Figure 4B). DEGs detected 14-42 DPC played a role in "oxidative phosphorylation" and "ATP metabolic process" including (ATP5F1D, ATP5F1E, COX2, NDUFA2) ( Figure 4A).
As described above, transcriptional changes detected in the day -7 group following MARV challenge enriched to GO processes associated with viral life cycle, myeloid cell differentiation, and apoptotic signaling pathway ( Figure 4C). DEG that enriched to GO term "viral life cycle" included genes that play a role in viral budding (VPS28, VPS37B, GRK2), and Interferon-stimulated genes and chemoattractant (ISG15, ISG20, CXCL8, CCL5). DEGs important for adaptive immunity are detected 9 DPC with DEGs enriching to "T cell activation", and "B cell activation" including DEGs that played a role in antigen presentation as well as T and B cell activation (CD74, HLA-DRB5, HLA-C, CD247, CD19, CD79A, IGHA1, and CD27) ( Figure 4D). Interestingly, DEGs enriching to GO terms associated with lymphocyte activation were down regulated 9 DPC ( Figure 4D). Finally, DEGs detected 14-42 DPC played a role in "myeloid cell differentiation" including ITGA2B, SLC4A1, STAT5B, H3C11, and EPB42 ( Figure 4C).
To further capture any transcriptional differences between the 3 surviving animals and the animal that succumbed to MARV challenge, we compared the WB transcriptional landscape between these animals at 0 DPV. This comparison identified 55 upregulated genes in the surviving relative to non-surviving animal on the day of vaccination some of which play a role in myeloid leukocyte activation (CD14, GRN, HEXIM1, and SIGLEC5); regulation of MAPK cascade (INPP5K, CEACAM1, and MINK1); and regulation of cell adhesion (ST3GAL4, ADAM8, and LAMA5) (Figures S3E, F).  (20) "lymphocyte activation", "T cell activation", and "B cell activation." Each column represents the median normalized transcript counts (RPKM) for each gene at each time point. The range of colors is based on scaled and centered RPKM values of the entire set of genes, with red indicating highly expressed genes and blue indicating lowly expressed genes. (E) maSigPro analysis identified four distinct and significant temporal expression clusters of genes. The median RPKM of each cluster over time is plotted. Statistical significance among time points was determined using one-way ANOVA with Dunnett's multiple-comparison test relative to 0 DPV. Statistical significance is indicated as **p < 0.01, and *p < 0.05. (F) Functional enrichment of the four clusters presented in (E). Horizontal bars represent the number of genes mapping to each GO term with color intensity representing the negative log of the FDR-adjusted p-value -log(q-value).

DISCUSSION
In a previous study we demonstrated that a single dose of this VSV-based MARV vaccine administered 35 days prior to challenge is highly efficacious protecting 100% of vaccinated NHPs from MVD (8). Here, we determined the minimum time between vaccination and challenge with a single dose of 1x 10 7 PFU. NHPs vaccinated 14 and 7 days prior to challenge with a lethal dose of MARV were uniformly protected from lethal MVD, while vaccination 3 days before challenge resulted in 75% survival. A previous study with VSV-MARV demonstrated the potential of this vaccine to be used as an effective postexposure treatment. When VSV-MARV was administered 20-30 min after lethal MARV challenge, all treated NHPs survived the challenge without signs of MVD (9). This result highlighted the fast-acting potential of this vaccine which was investigated in detail here. In addition, our observations are very similar to a study conducted with the FDA-approved VSV-EBOV vaccine against EBOV in NHPs where 67% of animals survived when animals were vaccinated 3 days before challenge (12), indicating that this platform can provide rapid protection and is suitable for deployment during an outbreak as has been described for VSV-EBOV (13).
Protection against filovirus infection and disease using VSVbased vaccines is mainly mediated by antigen-specific humoral immune responses (15,24). The dynamics of MARV GP-specific IgM and IgG responses demonstrated here were very similar to those reported for VSV-EBOV peaking around 10-14 days post vaccination (12). Although VSV-EBOV induced neutralizing IgM antibodies in humans (25), neutralizing antibody responses were not examined here as previous studies showed that they do not correlate with outcome after MARV infection (26). The development of the antibody response correlates with the induction of genes important for lymphocyte activation and function 11 and 14 DPV. Specifically, differential gene expression revealed a significant over-representation of genes that enrich to GO terms associated with B cell activation, humoral responses, T cell migration, and lymphocyte mediated immunity. Notable DEGs detected at these time points encoded heavy chains and factors important for memory B cell differentiation. These observations are in line with those recently reported for WB analysis following VSV-EBOV vaccination as well as PBMC analysis following VSV-MARV vaccination (8).
Since previous studies have shown limited contributions of T cells to vaccine-mediated protection against filoviruses (15, 27, 28), we did not measure the frequency of MARV GP-specific T cells in this study. Nevertheless, gene expression changes associated with T cell activation were evident 11-14 DPV. For instance, expression of genes important for TCR signaling and T cell survival were upregulated 11 DPV, and their expression remained increased at 14 DPV. Additionally, expression of genes important for antigen presentation, especially those encoding MHC molecules increased significantly 11 DPV. Overall, these findings support a role for T cells in mediating protection, most likely by providing help for affinity maturation and class switching of B cells as well as other potentially anti-viral functions as indicated by previous T cell depletion studies for VSV-EBOV (15,29).
VSV-MARV vaccination also induced a tightly regulated innate immune response. Several ISGs, Dead Box helicases and nucleic acid sensors were acutely induced 3 DPV. Interestingly, this induction was less robust than that observed for VSV-EBOV as indicated by data presented here and published elsewhere (29)(30)(31). In contrast, antiviral innate immune responses in control animals and those that succumb to infection are progressively amplified as disease progresses and are associated with adverse outcomes and high mortality both in NHP and humans (3,4,16,(32)(33)(34)(35). Although circulating levels of type I IFN do not significantly change with VSV-based vaccination, it is possible that the local production of IFNa/b in draining lymph nodes is sufficient to induce substantial expression of these ISGs (36). Moreover, control animals generated an excessive inflammatory response punctuated by substantial upregulation of genes encoding cytokines and chemokines and those important for neutrophil activation and function. This is in line with previous reports of neutrophilia and cytokine storm following lethal challenge with both EBOV and MARV (8,12,16,27,(37)(38)(39)(40). Finally, robust induction of genes that play a role in wound healing was also detected in control animals, indicative of the coagulopathy that accompanies EVD and MVD (37,39,40).
Innate antiviral responses have been shown to mediate the rapid protection using the VSV-EBOV vaccine in NHPs as determined by transcriptional analysis (28). Similarly, animals that were vaccinated 3 days before lethal challenge and survived exhibited hallmarks of robust innate immune responses that were attenuated by 14 DPC in line with resolution of the mild to moderate MVD symptoms they experienced. In contrast, the animal that succumbed to infection experienced sustained innate immune activation, and significant upregulation of genes that play a role in neutrophil activation. Interestingly, expression of genes important for myeloid cell activation and innate immunity was lower in blood samples obtained from this animal before vaccination. These observations, albeit preliminary, suggest that baseline expression of components of the innate immune branch may regulate individual responses to vaccination (and infection) and may modulate susceptibility to disease. This question merits additional investigation.
Another important observation in this study is that the robust innate immune responses induced by VSV-EBOV failed to protect animals from lethal MARV challenge when administered close to the time of challenge. In other words, the induction of ISGs by itself was not sufficient to "buy time" for antigen-specific responses to develop in animals vaccinated close to challenge. These findings suggest that MARV GP modulates the quality of the innate immune response and that antigenspecific responses induced by the vaccine rather than the virus under those circumstances are critical for shaping a protective immune response.
In summary, this study demonstrates that the VSV-MARV vaccine provides uniform protection within 7 days after administration of a single dose. In addition, 75% of the NHPs vaccinated only 3 days prior to challenge survived the lethal MARV infection. This rapid protection described here is comparable to that observed following administration of the approved VSV-EBOV vaccine (12) highlighting the fast-acting potential of this vaccine platform. VSV-EBOV was successfully used in a ring-vaccination approach in 2015 during the West African EBOV epidemic (13). In light of the first ever identified MARV case in Guinea in August 2021 (41) this study promotes the use of the VSV-MARV as an emergency vaccine suitable for ring-vaccination approaches.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI SRA BioProject, accession no: PRJNA763231.

ETHICS STATEMENT
The animal study was reviewed and approved by the Rocky Mountain Laboratories Animal Care and Use Committee, NIAID.

AUTHOR CONTRIBUTIONS
AM designed the study. AM and IM secured funding. AM, FF, and PH performed the animal studies. AM, AJ, ARM, JC, KO'D, and AP processed the samples and performed assays. AM, AJ, ARM, KO'D, AP, and IM analyzed the data. AM and IM wrote the manuscript. All authors contributed to the article and approved the submitted version. Resources and the National Center for Advancing Translational Sciences, NIH, through grant UL1 TR001414 awarded to IM.